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ABSTRACT 


In  this  thesis,  the  separation  bubble  behavior  and  its  effect  on  the  steady 
and  dynamic  stall  characteristics  of  a  thin  airfoil  in  a  compressible  flow  at  a 
transitional  Reynolds  number  was  studied.  For  such  flows,  laminar  separation 
occurs  near  the  airfoil  leading  edge,  but  turbulent  reattachment  occurs  within  a 
short  distance  downstream,  forming  a  separation  bubble  in  the  underlying  region. 
Two  experimental  techniques,  point  diffraction  interferometry  (PDI)  and  laser 
doppler  velocimetry  (LDV),  were  used  to  acquire  detailed  flowfield  information 
that  showed  the  development  of  the  leading-edge  separation  bubble  and  its 
subsequent  bursting  at  higher  angles  of  attack.  The  initiation  of  the  stall  process 
from  the  leading-edge  separation  bubble  as  opposed  to  trailing-edge  flow  reversal 
pointed  to  the  need  for  transitional  flow  analysis.  Both  in  the  boundary  layer  and 
Reynolds-averaged,  Navier-Stokes  (N-S)  analysis  methods,  transition  models  were 
incorporated  to  determine  the  location  and  extent  of  the  transition  zone  that  best 
modeled  the  measured  separation  bubble  behavior.  Computed  results  for  steady 
flow  gave  remarkable  agreement  with  the  measurements.  The  computations 
compared  favorably  with  the  measurements  for  an  airfoil  oscillating  in  pitch  about 
the  quarter-chord  point  during  the  airfoil  upstroke.  However,  the  computations 
did  not  predict  the  light  stall  and  vorticity-shedding  process  that  was  measured 
during  the  airfoil  downstroke. 
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I.         INTRODUCTION 

It  is  the  objective  of  this  thesis  investigation  to  contribute  to  the 
understanding  and  prediction  of  the  dynamic  stall  phenomenon,  which  is  observed 
to  occur  when  an  airfoil  is  pitched  rapidly  past  the  dynamic  stall  angle.  Dynamic 
stall  may  be  encountered  on  helicopter  and  propeller  blades,  on  fighter  aircraft 
wings  and  on  turbomachinery  blades.  Usually,  dynamic  stall  is  an  event  one  seeks 
to  avoid  because  it  may  generate  unacceptably  high  stresses,  which  may  cause 
blade  failures.  However,  more  recently  it  has  also  been  recognized  that  the  brief 
period  of  increased  lift  that  occurs  during  the  dynamic  stall  process  might  be 
usable  on  fighter  aircraft  to  improve  maneuverability.  Hence,  there  is  a  continuing 
need  and  interest  in  a  number  of  important  technical  fields  to  improve  our 
knowledge  of  the  dynamic  stall  phenomenon. 

In  general,  one  distinguishes  between  two  types  of  dynamic  stall.  Deep 
dynamic  stall  occurs  when  the  airfoil  is  pitched  to  angles  significantly  above  the 
static  stall  angle,  and  light  dynamic  stall  occurs  when  the  static  stall  angle  is 
exceeded  only  slightly.  In  both  cases,  a  so-called  dynamic  stall  vortex  is  observed 
to  form  near  the  leading  edge  of  the  airfoil,  which  grows  in  size  while  moving 
toward  the  trailing  edge.  During  this  phase  of  the  dynamic  stall  process,  the  lift 
continues  to  increase  quite  substantially  above  the  maximum  static  lift  value 
achievable  for  a  given  airfoil.  However,  as  soon  as  the  dynamic  stall  vortex 
reaches  the  trailing  edge,  the  lift  starts  to  break  down  quite  abruptly.  This  is 
accompanied  by  a  sharp  spike  in  the  pitching  moment.  This  general  sequence  of 
events  is  illustrated  in  Figure  1.1.  Naturally,  the  precise  physics  of  the  dynamic 
stall  phenomenon  is  much  more  complicated  and  depends  on  many  important 
parameters,  such  as  airfoil  shape,  Reynolds  number,  Mach  number,  etc.  To  review 
the  dynamic  stall  physics  in  more  detail,  refer  to  the  recent  surveys  of  Carr  [Ref.  1] 
and  Carr  and  Chandrasekhara  [Ref.  2]. 

In  the  past,  most  experimental  investigations  of  dynamic  stall  were 
restricted  to  incompressible  flow  studies  using  either  water  tunnels  because  of 
their  superior  flow  visualization  capability  or  relatively  inexpensive  low-speed 
wind  tunnels.  However,  it  became  increasingly  apparent  that  the  compressibility 


of  the  air  may  have  an  important  effect  on  the  dynamic  stall  process  even  though 
the  free-stream  Mach  number  may  be  less  than  0.3  (normally  cited  as  the  threshold 
beyond  which  compressibility  has  to  be  considered  in  steady-state  flow  analysis). 
Therefore,  a  special  dynamic  stall  wind  tunnel  was  built  in  the  Fluid  Mechanics 
Laboratory  of  the  National  Aeronautics  and  Space  Administration  (NASA)  Ames 
Research  Center,  Moffett  Field,  Calif,  which  made  it  possible  to  raise  the  free- 
stream  Mach  number  to  0.45.  Furthermore,  the  airfoil  was  mounted  such  that  it 
permitted  an  unobstructed  view  of  the  dynamic  stall  flowfield.  This  made  it 
possible  to  visualize  and  measure  the  precise  flow  details  and  events  that  lead  to 
the  formation  of  the  dynamic  stall  vortex. 
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Figure  1.1  Dynamic  Stall  Events  on  the  NACA  0012  Airfoil  [Ref.  1], 


In  this  thesis,  a  National  Advisory  Committee  for  Aeronautics  (NACA) 
0012  airfoil  is  studied  in  a  flow  of  Mach  number  0.3  and  Reynolds  number 
540,000.  This  choice  of  Mach  and  Reynolds  number  allows  the  study  of 
compressibility  effects  in  a  Reynolds  number  range  where  it  is  known  that  a 
separation  bubble  forms  near  the  leading  edge  of  the  airfoil.  The  central  question 
addressed  in  this  thesis  investigation  is  the  study  of  separation  bubble  behavior 
near  static  and  dynamic  stall  for  the  purpose  of  understanding  the  flow  physics, 
which  leads  to  the  initiation  of  static  and  dynamic  stall.  To  this  end,  a  new  flow 
visualization  technique,  point  diffraction  interferometry  (PDI)  was  used  to 
measure  the  density  gradients  in  the  flow,  and  laser-doppler  velocimetry  (LDV) 
was  applied  to  measure  the  velocity  field  near  the  leading  edge.  The  second 
question  addressed  is  the  assessment  of  currently  available  computational  methods 
concerning  their  adequacy  or  inadequacy  to  predict  the  type  of  flow  described. 

Therefore,  the  results  of  the  investigation  are  documented  as  follows: 
Chapter  II  provides  a  review  of  airfoil  flows  with  separation  bubbles,  Chapter  III 
contains  a  description  of  the  experimental  setup  and  measurement  techniques  used 
in  this  investigation,  Chapter  IV  provides  a  characterization  of  the  computational 
methods,  Chapter  V  summarizes  the  experimental  results,  and  Chapter  VI 
documents  the  computational  results,  and  finally,  Chapter  VII  gives  data 
comparisons  between  the  measurements  and  the  predictions.  Chapter  VIII 
presents  the  conclusions  and  recommendations. 


EL         REVIEW  OF  LOW-SPEED  AIRFOIL  FLOWS 
WITH  LEADING-EDGE  SEPARATION  BUBBLES 

A.  INTRODUCTION 

Information  about  the  stalling  characteristics  of  airfoils  is  extremely 
important  for  the  design  of  aircraft  wings,  helicopter  blades,  propellers,  turbo- 
machinery  blades,  windmills,  etc.  In  the  early  1930s,  Jones  [Ref.  3]  suggested  to 
distinguish  between  trailing-edge  stall  and  two  types  of  stall  caused  by  flow 
separation  near  the  airfoil  leading  edge.  In  the  late  1940s,  McCullough  and  Gault 
[Ref.  4]  performed  detailed  measurements  and  confirmed  the  existence  of  three 
types  of  stall,  namely  trailing-edge  stall  characterized  by  the  movement  of  the 
separation  point  of  the  turbulent  boundary  layer  forward  from  the  trailing  edge 
with  increasing  incidence;  leading-edge  stall  caused  by  an  abrupt  separation  of 
the  laminar  boundary  layer  near  the  leading  edge  without  subsequent 
reattachment;  and  thin-airfoil  stall  where  laminar-flow  separation  occurs  near  the 
leading  edge,  followed  by  turbulent  flow  reattachment  that  moves  progressively 
downstream  to  the  trailing  edge  with  increasing  incidence.  As  a  result,  a 
recirculatory  flow  region  forms  between  the  separation  and  reattachment  points 
that  is  generally  referred  to  as  a  separation  bubble. 

It  is  the  objective  of  the  present  investigation  to  make  a  contribution  to  the 
understanding  of  the  physics  of  the  separation  bubble.  Therefore,  the  review  of 
previous  work  is  limited  to  studies  of  airfoil  flow  at  Mach  numbers  equal  to  or  less 
than  0.3  where  leading-edge  separation  bubbles  were  found  to  occur. 

B.  FUNDAMENTAL  FLOW  PHYSICS  OF  LEADING-EDGE 
SEPARATION  BUBBLES 

The  type  of  bubble  that  forms  depends  on  the  airfoil  shape.  In  a 
comprehensive  review  of  separation  bubbles,  Tani  [Ref.  5]  drew  attention  to  the 
need  to  further  distinguish  between  long  and  short  bubbles.  The  measurements 
of  McCullough  and  Gault  revealed  the  presence  of  a  separation  bubble  near  the 
leading  edge  of  the  NACA  63-009  airfoil  at  an  incidence  of  4°  and  a  Reynolds 
number  of  5.8  x  106.    Separation  of  the  laminar  boundary  layer  took  place  before 


stall  near  the  leading  edge  of  an  airfoil  at  low  incompressible  flow  speeds  to  form 
the  separation  bubble  even  at  a  Reynolds  number  greater  than  1  million.  An 
increase  in  incidence  caused  the  bubble  to  move  forward  and  contract  in 
streamwise  extent,  which  was  on  the  order  of  1  %  airfoil  chord.  A  small  region  of 
constant  pressure  existed  within  the  bubble.  At  higher  angles  of  incidence,  the 
flow  no  longer  reattached  to  the  surface  downstream  of  the  laminar  separation,  or 
in  effect,  experienced  short  bubble  burst.  This  occurred  at  an  incidence  of  9  °  and 
caused  a  complete  redistribution  of  the  pressure  and  an  abrupt  loss  of  lift.  In 
contrast,  a  separation  bubble  first  occurred  at  low  incompressible  speeds  and  the 
same  Reynolds  number  of  5.8  x  106  on  the  NACA  64A006  airfoil  at  3  °  incidence. 
As  the  incidence  was  increased  to  5°,  the  bubble  broke  down,  causing  a  slight 
drop  in  lift.  However,  this  bubble-burst  did  not  lead  to  complete  flow  separation. 
Instead,  the  flow  separation  persisted  over  some  distance  and  reattached  farther 
downstream.  As  the  incidence  was  increased,  the  reattachment  point  moved 
progressively  downstream  until  it  reached  the  trailing  edge  at  about  9°  incidence 
angle.  This  is  referred  to  as  the  long  bubble. 

Tani  suggested  to  distinguish  short  from  long  bubbles  by  their  effect  on 
the  airfoil  pressure  distribution  rather  than  by  their  different  lengths.  The 
presence  of  separation  bubbles  makes  the  pressure  distribution  markedly  different 
from  that  in  inviscid  flow.  In  particular,  for  a  long  bubble  a  sharp  suction  peak 
near  the  leading  edge  is  not  realized.  Instead,  a  pressure  plateau  of  reduced  level 
extends  over  the  length  of  the  separation  bubble.  On  the  other  hand,  the  suction 
peak  continues  to  increase  followed  by  a  small  pressure  plateau  up  to  the  stall 
angle  of  attack  for  a  short  bubble  in  a  high  Reynolds  number  flow.  For  the  flow 
of  consideration  in  this  thesis  at  a  Reynolds  number  540,000  and  M^  =0.3,  the 
leading-edge  short  bubble  is  shown  to  have  a  more  extensive  effect  on  the 
pressure  distribution. 

The  formation  of  a  separation  bubble  clearly  depends  on  the  existence  of 
an  adverse  pressure  gradient  that  is  steep  enough  to  cause  laminar  separation. 
However,  no  bubble  will  form  if  the  Reynolds  number  is  sufficiently  high  and  the 
flow  transitions  from  laminar  to  turbulent  ahead  of  the  point  where  flow 
separation  would  have  occurred  if  the  boundary  layer  had  remained  laminar.  On 
the  other  hand,  at  a  sufficiently  low  Reynolds  number  the  separated  flow  will  be 


unable  to  reattach  and,  again,  no  bubble  will  be  formed.  Therefore,  a  Reynolds 
number  range  exists  within  which  separation  bubbles  form.  Below  the  critical 
Reynolds  number  for  bubble  formation  the  lift  is  low  and  the  drag  is  high.  As 
shown  by  McMasters  and  Henderson  [Ref.  6],  the  maximum  L/D  of  an  airfoil 
improves  dramatically  in  the  range  of  Reynolds  numbers  from  104  to  106,  which  is 
shown  in  Figure  2.1  [Ref.  7].  For  Re  >  106,  typical  of  large  aircraft,  boundary 
layer  transition  occurs  upstream  of  the  theoretical  laminar  separation  point.  As  a 
consequence,  the  turbulent  boundary  layer  is  much  more  resistant  to  separation, 
and  the  flow  remains  attached  to  relatively  high  incidence  angles.  Eventually,  the 
stall  will  occur  as  a  result  of  the  forward  movement  of  the  separation  point  from 
the  trailing  edge. 

In  the  range  of  Reynolds  numbers  from  104  to  106,  many  complicated 
phenomena  take  place  within  the  boundary  layer.  Separation,  transition,  and 
reattachment  can  all  occur  within  a  short  distance.  As  shown  in  Figure  2.1, 
surface  roughness  has  an  important  effect.  Similarly,  the  turbulence  level  of  the 
free  stream  is  known  to  affect  the  aerodynamic  characteristics. 
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Figure  2.1.  Airfoil  Performance  as  a  Function  of  Chord  Reynolds  Number 
[Ref.  7]. 


The  basic  flow  physics  of  bubble  formation  appears  to  have  been  first 
explained  by  von  Doenhoff  in  NACA  TN  639  of  March  1938  [Ref.  8].  He 
suggested  that  the  onset  of  turbulence  in  the  separated  flow  is  responsible  for  the 
reattachment  because  the  turbulent  entrainment  caused  the  flow  to  return  to  the 
airfoil  surface.  The  typical  structure  of  a  laminar  separation  bubble  was  described 
by  Horton  in  1968  [Ref.  9].  Figure  2.2  shows  the  flowfield  in  the  vicinity  of  a 
laminar  separation  bubble.  At  the  separation  point,  the  dividing  streamline  is 
observed  to  rise  slowly  away  from  the  wall.  After  the  boundary  layer  undergoes 
transition  to  turbulent  flow,  the  streamline  quickly  reattached  to  the  wall.  The 
center  of  the  reversed  flow  vortex  lies  near  the  reattachment  point  and  causes  the 
fluid  to  move  faster  in  this  region.  Within  the  separated  region,  the  wall  pressure 
remains  constant  near  the  separation  point,  while  a  steep  pressure  gradient 
accompanies  the  high-velocity  region  near  the  reattachment  point.  Tani  noted 
that  the  size  of  the  separation  bubble  was  reduced  as  the  adverse  pressure 
gradient  increased  [Ref.  5].  This  was  explained  by  the  more  rapid  boundary  layer 
transition. 
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Figure  2.2.  Structure  of  Flow  With  Laminar  Separation  Bubble  [Ref.  7]. 


Tani  postulated  that  the  Reynolds  number  based  on  the  boundary  layer 
thickness  at  the  separation  point  provides  a  better  measure  than  the  chord-based 
Reynolds  number.     He  suggested,  based  on  available  experiments,  that  the 


Reynolds  number  based  on  the  displacement  thickness  at  separation  must  exceed 
780  if  reattachment  is  to  occur.  Later,  Owen  and  Klanfer  [Ref.  10]  put  forward  a 
similar  criterion.  This  Tani,  Owen,  Klanfer  criterion  can  be  interpreted  as  an 
assurance  that  transition  will  occur  in  the  separated  flow. 

Crabtree  [Ref.  1 1]  sought  to  develop  a  criterion  for  bubble  breakdown.  He 
made  the  assumption  that  the  breakdown  occurs  because  a  maximum  possible 
pressure  value  exists  that  can  be  recovered  in  the  turbulent  entrainment  process, 
which  causes  the  flow  reattachment.  He  found  from  the  analysis  of  the  available 
experimental  results  that  the  value  of  the  pressure  recovery  coefficient  increases 
with  either  an  increase  in  incidence  at  constant  Reynolds  number  or  with  a 
decrease  in  Reynolds  number  at  constant  incidence.  The  latest  measurements  by 
Malkiel  and  Mayle  [Ref.  12]  showed  that  the  intermittency  development  in 
separation  bubbles  is  well  modeled  with  the  turbulent  spot  theory  used  in 
attached  boundary  layers. 

The  foregoing  discussion  of  the  main  features  of  separation  bubbles  makes 
it  clear  that  the  transition  to  turbulent  flow  has  a  controlling  effect  on  the  extent 
of  laminar  separation  bubbles.  At  chord  Reynolds  numbers  less  than  106,  the 
local  Reynolds  number  based  on  the  boundary  layer  development  is  often  too 
low  for  the  Tollmien-Schlichting  instabilities  to  promote  natural  transition  before 
the  laminar  boundary  layer  reaches  separation.  Once  laminar  separation  has 
occurred,  the  resulting  inflectional  velocity  profiles  promote  a  more  rapid 
amplification  of  boundary  layer  fluctuations,  which  eventually  reach  transition 
levels.  Laminar  separation  is  also  accompanied  by  a  significant  increase  in  the 
boundary  layer  thickness.  After  transition,  increased  mixing  promotes  growth  of 
the  turbulent  separated  shear  layer,  which  eventually  reattaches  to  the  airfoil 
surface.  A  decrease  in  Reynolds  number  tends  to  delay  transition  within  the 
separation  bubble.  On  the  other  hand,  the  addition  of  external  disturbances,  such 
as  free  stream  turbulence,  acoustic  disturbances,  or  surface  roughness  promote 
transition  within  the  separation  bubble.  Accordingly,  the  transition  onset 
location  must  be  considered  a  critical  parameter  and  its  prediction  becomes  a 
crucial  element  in  the  development  of  prediction  methods  for  airfoil  flows  with 
leading-edge  separation  bubbles. 


C.  TRANSITION  ONSET  PREDICTION 

A  widely  used  empirical  criterion  for  the  onset  of  boundary  layer  transition 
in  high  Reynolds  number  airfoil  flows,  Re  >  106,  has  been  given  by  Michel  [Ref. 
13].  Early  attempts  at  transition  prediction  in  the  Reynolds  number  range  from 
104  to  106  focused  on  correlation  of  local  flow  parameters  and  local  Reynolds 
numbers,  proposed  by  Gault  [Ref.  14],  Gaster  [Ref.  15],  Horton  [Ref.  9],  O'Meara 
and  Mueller  [Ref.  16].  More  recently,  linear  stability  calculations  have  been 
introduced.  Such  calculations  assume  an  en  transition  criterion  based  on  the 
stability  of  Falkner-Skan  reverse-flow  profile  solutions.  This  approach  assumes 
that  the  onset  of  transition  can  be  correlated  with  the  location  where  the 
amplification  of  the  Tollmien-Schlichting  instabilities  has  reached  a  certain  level, 
expressed  by  the  value  of  (n)  in  en.  This  value  is  usually  assumed  to  be  in  the 
range  from  8  to  9.  Additional  details  about  this  method  can  be  found,  for 
example,  in  Cebeci  and  Bradshaw  [Ref.  17].  Another  promising  procedure  for  the 
prediction  of  transition  has  been  recently  developed  by  Herbert  and  Bertolotti 
[Ref.  18]  and  is  known  as  the  parabolized  stability  equation  (PSE)  method. 

D.  TRANSITION  LENGTH  PREDICTION 

The  key  variable  in  describing  a  boundary  layer  during  transition  from 
laminar  to  turbulent  flow  is  the  intermittency,  ytr(x),  defined  as  the  fraction  of 

time  that  the  flow  is  turbulent.  Therefore,  the  transition  zone  may  be  defined  as 
beginning  where  ytr(x)  has  just  departed  from  zero  and  extending  to  where  it  is 

unity.  In  constant  pressure  boundary  layers,  a  large  amount  of  experimental  data 
has  been  gathered  and  a  satisfactory  model  based  on  Emmons'  theory  of  spots 
and  Narasimha  and  Dhawan's  [Ref.  19]  hypothesis  of  concentrated  breakdown 
have  been  developed.  Chen  and  Thyson  [Ref.  20]  extended  this  model  to 
account  for  the  effect  of  pressure  gradient  by  making  the  hypothesis  that  the 
turbulent  spot  propagation  velocity  at  any  point  x  is  proportional  to  the  local 
free-stream  velocity,  but  the  transition  growth  of  the  spot  across  streamlines  is 
assumed  to  be  unaffected.  This  Chen-Thyson  transition  model  is  widely  used  for 
attached  high  Reynolds  number  airfoil  flows,  as  documented  by  Cebeci  and 
Bradshaw  [Ref.  17].  The  effect  of  pressure  gradient  is  to  lengthen  the  transition 
zone  in  favorable  gradient  pressure  flows  and  to  shorten  it  in  adverse  pressure 
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gradient  flows.  This  effect  is  incorporated  in  the  Chen-Thyson  model  by  the 
choice  of  the  transition  length  value.  A  lower  transition  length  value  forces 
transition  to  take  place  over  a  shorter  distance.  Recently,  Cebeci  [Ref.  21]  has 
proposed  a  correlation  formula  to  include  the  effect  of  Reynolds  number  so  that 
the  Chen-Thyson  model  can  be  applied  in  the  Reynolds  number  range  from  104 
to  106.  Another  transition  length  model  that  attempts  to  account  for  the  effects  of 
free-stream  turbulence  and  pressure  gradients  has  been  suggested  by  Abu- 
Ghannam  and  Shaw  [Ref.  22].  It  assumes  the  existence  of  similarity  in  the  spot- 
generation  rate.  Very  recently,  Gostelow  et  al.  [Ref.  23]  developed  a  new  model 
to  calculate  the  transition  length  as  a  function  of  pressure  gradient  and  free- 
stream  turbulence  level,  which  is  based  on  extensive  new  measurements  of  the 
spot-generation  rate  in  attached  boundary  layers. 

It  is  important  to  recall  that  no  transition  measurements  exist  as  yet  in  free 
shear  layers.  Hence,  the  use  of  the  Chen-Thyson,  Abu-Ghannam  and  Shaw,  or 
Gostelow  and  Walker  transition  length  models  for  the  modeling  of  the  transition 
process  in  separation  bubbles  is  tentative  at  best.  For  this  reason,  in  the  present 
study  the  sensitivity  of  leading-edge  separation  bubble  calculations  to  the  choice 
of  both  the  transition  onset  location  and  the  transition  length  is  examined. 

E.        EFFECT  OF  AIRFOIL  OSCILLATION  ON  LEADING-EDGE 
SEPARATION  BUBBLES 

It  is  well  known  that  the  dynamic  stall  phenomena  caused  by  oscillating  an 
airfoil  about  a  mean  angle  of  attack  that  is  close  to  the  static  stall  angle  are  quite 
different  from  the  static  stall  phenomena.  For  a  review  of  the  extensive  literature 
on  dynamic  airfoil  stall,  the  author  refers  the  reader  to  the  review  papers  by 
McCroskey  [Ref.  24],  Carr  [Ref.  1],  Carr  and  McCroskey  [Ref.  25],  and  Carr  and 
Chandrasekhara  [Ref.  26]  that  discuss  the  available  experimental  and 
computational  results.  It  is  apparent  from  this  literature  survey  that  little  work  has 
been  devoted  to  the  study  of  the  so-called  "light  dynamic  stall"  problem  in  the 
Reynolds  number  regime  from  104  to  106  in  order  to  understand  the  effect  of  small 
amplitude  oscillations  on  the  separation  bubble  behavior.  Therefore,  it  was  a 
major  objective  of  the  present  work  to  investigate  this  effect.  To  this  end,  a 
NACA  0012  airfoil  was  oscillated  in  pitch  about  the  quarter-chord  in  a  free-stream 
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flow  at  a  Mach  number  of  0.3  and  a  Reynolds  number  of  540,000.  Two 
experimental  techniques  were  used  to  acquire  detailed  flowfield  information, 
namely  point  diffraction  interferometry  (PDI)  and  laser  doppler  velocimetry 
(LDV).  These  experimental  methods  and  results  are  documented  in  Chapters  III 
and  V,  respectively.  The  computational  analysis  methods  and  results  are 
described  in  Chapters  IV  and  VI.  The  comparisons  between  the  computations 
and  experiments  are  given  in  Chapter  VII. 

F.        ANALYSIS  METHODS  FOR  AIRFOIL  FLOWS  WITH 
SEPARATION  BUBBLES 

Over  the  past  30  years,  so-called  viscous-inviscid  interaction  methods  have 
been  developed  to  overcome  the  limitation  of  classical  boundary-layer  methods  in 
the  analysis  of  attached  airfoil  flows  due  to  the  breakdown  of  the  equations  at 
the  point  of  flow  separation.  Instead  of  computing  the  pressure  distributions  from 
an  inviscid  flow  computation  and  using  it  as  input  into  the  boundary  layer 
computations,  the  coupling  between  the  viscous  and  inviscid  parts  of  the 
flowfield  is  achieved  by  an  inverse  boundary  layer  calculation.  Various  such 
coupling  methods  have  been  developed  for  the  calculation  of  steady  airfoil  flows. 
The  two  most  widely  used  viscous-inviscid  interaction  codes  were  developed  by 
Cebeci  [Ref.  27]  and  the  airfoil  design-analysis  code  developed  by  Drela  and 
Giles  [Ref.  28]  at  Massachusetts  Institute  of  Technology  (MIT)  at  Cambridge, 
Mass.  Example  calculations  of  steady  separation  bubbles  at  Reynolds  numbers 
105-106  have  been  given  by  Evangelista  et  al  [Ref.  29],  using  the  Drela  and  Giles 
code  and  by  Walker  et  al.  [Ref.  30],  using  the  Cebeci  code.  The  latter  authors 
found  that  the  separation  bubble  occurring  on  a  NACA  65-213  airfoil  at  a 
Reynolds  number  of  240,000  could  be  predicted  quite  well,  provided  the 
transition  length  constant  was  reduced  from  1 200  to  between  20  and  40,  thereby 
reducing  the  transition  length  in  the  Chen-Thyson  model. 

Relatively  little  work  has  been  done  to  extend  the  viscous-inviscid 
interaction  methods  to  the  analysis  of  the  flow  over  oscillating  airfoils  at  low 
Reynolds  numbers.  The  numerical  solution  of  the  Reynolds-averaged  Navier- 
Stokes  (N-S)  equations  is  becoming  an  increasingly  attractive  alternative  to  the 
viscous-inviscid  interaction  approaches  because  of  the  increasing  availability  and 
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speed  of  modern  workstations.  This  approach  offers  the  advantage  of  solving 
compressible,  steady  and  unsteady  attached,  mildly  separated  and  fully  separated 
flows  with  one  and  the  same  computational  method.  Viscous-inviscid  interaction 
approaches,  on  the  other  hand,  are  limited  to  flows  with  relatively  small  regions  of 
separation.  Therefore,  in  this  thesis,  the  N-S  approach  was  chosen  as  the  main 
method  for  the  analysis  of  the  flow  over  the  steady  or  oscillating  NACA  0012 
airfoil.  This  approach  is  supplemented  by  an  unsteady  inviscid,  incompressible 
panel  code  and  a  direct  boundary  layer  code.  These  three  codes  are  described  in 
Chapter  IV. 
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IE.      EXPERIMENTAL  SETUP  AND  MEASUREMENT  TECHNIQUES 

A.  INTRODUCTION 

A  novel  facility  known  as  the  Compressible  Dynamic  Stall  Facility  (CDSF) 
has  been  designed  at  the  NASA  Ames  Research  Center  Fluid  Mechanics 
Laboratory  (FML)  to  quantitatively  document  the  compressibility  effects  on 
dynamic  stall  of  rapidly  pitching  airfoils.  The  wind-tunnel  facility  is  in  use  for 
experimental  research  of  oscillating  airfoil  dynamic  stall  that  applies  to  helicopter 
rotors  [Ref.  31].  The  research  also  has  application  for  maneuvering  missile  fins 
and  aircraft  surfaces  in  addition  to  those  discussed  in  the  Introduction.  For  control 
or  eventual  beneficial  use  of  dynamic  stall,  a  clearer  understanding  of  the  flow 
phenomena  that  affect  the  inception  and  development  of  the  dynamic-stall  process 
is  required.  The  CDSF  test  section  uses  optical  glass  windows  to  support  an  airfoil 
with  pins  that  are  smaller  than  the  airfoil  thickness  to  give  an  unobstructed  view 
through  the  tunnel  of  the  flow  over  a  NACA  0012  test  airfoil  surface.  The  flow 
can  be  studied  near  the  airfoil  leading  edge  where  the  separation  bubble  forms  and 
the  dynamic  stall  process  occurs,  which  is  the  area  of  focus  for  this  work,  as  well 
as  away  from  the  airfoil  surface.  Point  diffraction  interferometry  (PDI)  and  laser 
doppler  velocimetry  (LDV)  measurements  are  used  to  characterize  (1)  flows  over 
steady  airfoils,  (2)  high-amplitude,  deep-dynamic  stall  airfoil  flows,  and  (3)  low- 
amplitude,  lightly  stalled  compressible  flows.  For  the  light  dynamic  stall,  the 
measurements  showed  that  the  viscous  effects  normal  to  the  airfoil  upper  surface 
were  on  the  order  of  the  airfoil  thickness,  giving  ample  data  to  evaluate  the 
capability  of  various  numerical  methods  to  predict  the  flow  characteristics  over  the 
leading  edge  of  the  airfoil.  A  major  part  of  this  work  investigated  the  effect  of 
small  amplitude  oscillations  on  the  separation  bubble  behavior. 

B.  EXPERIMENTAL  TEST  FACILITY 

The  CDSF  is  operated  as  one  of  four  in-draft  tunnels  at  the  FML  complex, 
driven  by  a  6-MW,  108-m3/s,  continuous-running  evacuation  compressor, 
described  in  detail  by  Carr  and  Chandrasekhara  [Ref.  32].  The  flow  in  the  tunnel 
is  controlled  by  a  choked,  variable-area  downstream  throat  capable  of  producing 
Mach  numbers  from  0  to  0.5,  allowing  the  study  of  a  wide  range  of  free-stream 
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conditions  from  incompressible  to  moderately  compressible.  The  evacuation 
compressor  drives  the  system  and  automatically  maintains  sufficient  vacuum 
downstream  of  each  test  facility  in  the  tunnel  complex  to  cause  sonic  velocity  in 
the  smallest  cross  section  of  each  test  chamber.  This  choked-throat  is  used  to 
isolate  each  test  bay  and  maintain  clean,  constant  flow  conditions  for  a  wide  range 
of  operating  conditions  during  simultaneous  operation  of  multiple  tunnels. 


The  CDSF  test  section  is  35-cm  high,  25-cm  wide,  and  100-cm  long.  The 
15.2-cm-diameter,  2.54-cm-thick,  D-shaped  interferogram-quality  glass  (BK-7) 
windows  are  supported  in  magnesium  frames  that  are  mounted  on  bearing  races, 
which  permit  rotation  of  the  window.  A  seal  prevents  leakage  of  air  into  the  test 
section.  Detailed  mechanical  design  information  on  the  drive  mechanism  can  be 
found  in  Sticht  [Ref.  33].  The  window/airfoil  combination  is  driven  in  sinusoidal 
oscillation  by  a  four-bar,  push-rod-flywheel  system,  about  the  25%  chord  position. 
A  variable-speed  AC-drive  motor  with  a  controller  is  used  to  maintain  the  airfoil 
oscillation  frequency  to  within  1%  of  the  desired  value.  The  drive  mechanism  is 
capable  of  oscillating  the  airfoil  about  mean  incidence  angles  of  0  to  15°, 
oscillation  amplitudes  of  2  to  10  °,  and  frequencies  up  to  100  Hz.  The  phase  angle 
in  the  cycle  is  determined  by  an  optical  encoder  keyed  to  the  flywheel.  Figure  3.1 
shows  a  schematic  of  the  drive  system  and  includes  specific  details  about  the  drive 


system. 
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Figure  3.1.  Side  View  of  the  Compressible  Dynamic  Stall  Facility  [Ref.  32]. 
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The  primary  instrumentation  for  this  study  was  non-intrusive,  including  a 
two-component  LDV  and  a  stroboscopic  PDI  system.  Three  high-resolution 
encoders  provided  accurate  position  and  frequency  information.  The  mean  angle 
of  attack  was  indicated  by  an  absolute  position  encoder  with  3600  counts  per 
revolution  (0.1  deg/count).  The  instantaneous  angle  of  attack  was  provided  by  an 
incremental  encoder  mounted  on  the  window.  The  frequency  and  phase  angle 
were  provided  by  another  incremental  encoder  with  800  counts/cycle  of  oscillation 
and  was  keyed  to  the  flywheel  drive  shaft.  All  encoder  outputs  are  digital  and 
were  read  into  a  Micro  VAX  II  computer,  which  was  used  to  control  the 
experiment.  [Ref.  34] 

C.        UNSTEADY  LDV  MEASUREMENT  TECHNIQUE 

A  TSI,  Inc.,  two-color,  two-component,  four-beam,  frequency-shifted, 
argon-ion  laser-based,  off-axis,  forward-scatter  system  was  used  for  velocity 
measurements.  Traversing  was  accomplished  by  directing  the  four  beams  onto  a 
352-mm  focal-length  lens  mounted  on  a  computer-controlled  traverse.  Figure  3.2 
presents  a  schematic  of  the  CDSF  instrumentation.  The  scattered  light  was 
collected  at  15  degrees  off-axis  in  the  forward  scatter  mode  by  the  500-mm  focal- 
length  receiving  optics.  This  arrangement  gave  a  reduced  probe  volume  length 
and  improved  the  measurement  resolution.  The  receiving  optics  were  mounted  on 
a  separate  traverse  mechanism,  but  were  slaved  to  the  traverse  on  the  transmitting 
side.  The  individual  photo-multiplier  signals  were  processed  by  TSI  1990 
counters. 

Special  phase-locking  circuitry  enabled  handling  of  the  random  LDV  data 
and  the  unsteady-position  data.  There  was  a  need  to  read  the  phase  angles  the 
instant  each  LDV  data  sample  was  validated,  which  were  simultaneous  Doppler 
signals  in  both  the  U-  and  V-components  of  velocity.  The  LDV  data  were 
acquired  in  the  coincidence  mode  with  the  acceptance  window  width  arbitrarily  set 
to  50  |is.  The  coincidence  pulse  was  used  to  trigger  data  acquisition  and  freeze  the 
rapidly  changing  encoder  values  until  data  transfer  to  the  computer  could  be 
completed  as  shown  in  Figure  3.3.  The  data  from  the  system  were  acquired  by  the 
computer  in  the  DMA  mode  using  an  eight-channel  multiplexer  to  check  for 
coincidence  of  the  data  on  two  channels.  The  encoder  data  were  also  input  to  the 
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multiplexer  through  a  system  that  records  the  encoder  value  only  when  coincident 
LDV  data  were  obtained.  The  LDV  data-acquisition  scheme  was  described  in 
detail  by  Chandrasekhara  and  Ahmed  [Ref.  34]. 
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Figure  3.2.  Schematic  of  the  LDV  System  in  Use  in  the  CDSF  [Ref.  32], 
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Figure  3.3.  Unsteady-Flow  LDV  Data- Acquisition-Method  Schematic  [Ref.  34], 
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The  LDV  data  rate  was  dependent  on  the  seeding  particles,  which  were  1- 
(im-diameter  polystyrene  latex  particles  (PSL)  suspended  in  alcohol  and  injected 
into  the  tunnel  inlet  by  a  TSI  9306  six-jet  atomizer.  The  particles  traveled  3 
meters  before  passing  through  the  probe  volume  in  the  test  section,  and  by  that 
time  the  alcohol  had  evaporated,  leaving  only  the  PSL. 

The  flow  Mach  number  for  the  tests  was  0.3  and  the  Reynolds  number  was 
540,000.  The  airfoil  oscillation  frequency  was  21.6  Hz  or  a  reduced  frequency  of 
0.1  based  on  the  full  chord.  The  flow  over  the  NACA  0012  airfoil  was 
investigated  at  steady  angles  of  attack  and  for  oscillations  about  the  quarter  chord 
point,  a  =  10°  - 10°  sinco  t  and  a  =  10°  -  2°  sinco  t  . 

Data  were  acquired  in  a  rectangular  x-y  grid.  The  origin  was  determined  at 
the  leading  edge  of  the  airfoil  when  the  angle  of  attack  was  0  degree.  The  LDV 
probe  volume  was  traversed  in  the  range  -0.167  <  x/c  <  0.167,  and  0.083  <  y/c  < 
0.167.  The  rectangular  measurement  grid  embedded  in  the  computational  grid  is 
shown  in  Figure  3.4.  The  resolution  in  both  the  x  and  y  direction  was  1.27  mm. 
The  measurement  grid  remained  stationary  as  the  airfoil  rotated  about  the  quarter - 
chord  point. 


Figure    3.4    Measurement  Locations  (6x21)  Overlaid  on  Computational  Grid 

(275x81).   a  =  8°,  NACA  0012  Airfoil. 
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D.        STROBOSCOPIC  PDI  MEASUREMENT  TECHNIQUE 


The  PDI  technique  permits  the  real  time  recording  of  flow  interf erograms. 
The  PDI  technique  as  used  in  the  present  test  is  based  on  a  variation  of  the  work 
by  Smartt  [Ref.  35  and  36].  The  PDI  technique,  described  by  Can  et  al.  [Ref.  37], 
used  a  "pin  hole"  point  discontinuity  in  an  exposed  photographic  plate  to  diffract  a 
portion  of  the  incoming  light.  The  plate  was  placed  at  the  focal  plane  of  the  point 
source,  producing  a  reference  spherical  wave  front  from  the  "pin  hole"  by 
diffracting  a  portion  of  the  incident  light.  The  light  passing  through  the 
photographic  plate  surrounding  the  pinhole  interacted  with  the  reference  wave 
front,  giving  the  interference  fringes  as  seen  in  Figure  3.5.  The  optical  path 
difference  caused  by  density  variations  in  the  flowfield  produced  the  interference 
fringes. 


Wavefront 
under  test 


image 


Figure  3.5.  Principle  of  Point  Diffraction  Interf erometry  [Ref.  38]. 

1.         Implementation 

The  interferogram  images  were  obtained  by  triggering  a  laser  light  source  at 
defined  phase  angles  of  the  airfoil  through  the  oscillation  cycle  to  give  the  density 
gradients  in  the  flow  field.  The  light  source  was  a  pulsed  Nd:YAG  laser  with  a 
frequency  doubler,  operating  at  a  wavelength  of  532  nm.  Figure  3.6  shows  a 
schematic  of  the  arrangement,  having  a  long  focal  length  that  permits  use  of  the 
standard  Z-type  schlieren  arrangement  with  minimum  astigmatism.  The  laser  light 
was  expanded  through  a  negative  lens  to  partially  fill  one  schlieren  mirror  and 
directed  through  the  test  section,  and  then  refocused  by  the  other  schlieren  mirror 
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before  passing  through  the  exposed  photographic  plate.  The  laser  replaced  the 
conventional  spark  as  the  light  source  and  was  also  pulsed  with  enough  energy  to 
burn  off  the  emulsion  located  at  the  focal  plane  of  the  second  mirror.  The  spot  or 
"pin  hole"  was  created  by  passing  light  through  the  test  section  at  a  no-flow 
condition.  [Ref.  37] 


Schiaran  Mirror 


Row  Direction 
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Figure  3.6.  Schematic  of  the  PDI  Optical  Arrangement  [Ref.  39]. 

The  tunnel  was  turned  on  and  the  real-time  interference  fringes  were 
recorded  on  Polaroid  film  (ASA  3000)  for  immediate  viewing.  This  capability  was 
of  great  value  for  analysis  of  the  dynamic  flow  field  under  investigation,  since  it 
permitted  rapid  review  of  the  progress  of  the  dynamic  stall  process  and  on-line 
retaking  of  any  photographs  that  were  not  of  the  required  quality.  This  also 
permitted  an  on-line  study  of  the  flow  field  as  it  developed,  since  any  phase  angle 
in  the  cycle  could  be  accessed  directly.  This  ability  to  actively  search  for  the  onset 
of  dynamic  stall  using  the  PDI  technique  is  in  strong  contrast  to  the  more 
conventional  holographic  interferometers,  which  require  a  major  post-processing 
effort  before  the  interferogram  can  even  be  verified,  much  less  analyzed. 

The  phase  angle  for  each  interferogram  was  chosen  by  setting  switches  on  a 
specially  designed  comparator  circuit  box.  Since  the  laser  took  a  finite  time  to 
actually  emit  the  light,  additional  electronics  were  incorporated  in  the  circuitry  to 
display  the  phase  angle  when  the  laser  actually  fired.  This  was  accomplished  by 
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sensing  the  light  using  a  photo  diode;  when  the  laser  fired,  the  output  of  the 
photodiode  was  used  to  freeze  the  angle-of-attack  display  so  that  the  precise  angle 
at  which  the  image  was  taken  could  be  recorded. 

2.         Fringe  Counting  To  Determine  Mach  Number 

The  PDI  interferometer  produces  images  that  are  equivalent  to  those 
produced  by  Mach-Zehnder  interferograms.  It  is  a  common  path  interferogram 
that  is  less  sensitive  to  vibration  and  uses  the  mirrors  of  an  existing  schlieren  setup 
to  direct  the  laser  beam.  The  path  length  difference  attributable  to  density  (or 
phase)  changes  can  be  related  to  the  fringe  number  e  for  a  two-dimensional  flow. 
The  equation  is 

e  =  (n-nr)i  (3.1) 

For  bright  fringes,  e  is  zero  or  an  integer.  The  dark  fringes  represent  half 
integers.  The  refractive  indexes,  n  and  nr,  are  for  the  signal  beam  and  the 
reference  beam,  respectively.  Using  the  perfect  gas  and  the  Gladstone-Dale 
equation  [Ref.  40],  the  above  expression  can  be  reduced  to 

p-p'=7^(T+A£ 

(n0-l)  V   L   J 

The  constant  A  can  easily  be  determined  for  the  present  experiment  because 
the  laser  wavelength,  X0  =  532  nm,  the  wind  tunnel  width,  L  =  25  cm,  the 
refractive  index  of  air,  (n0-l)  =  2.733  x  10-4  [Ref.  40],  and  the  standard 
atmospheric  density,  p0  =  1.21  kg/m3.  Thus,  by  dividing  the  above  equation  by  p0 
and  rearranging 

p      pr      0.009421  e  ,00. 

J—  =  ^-L  + (3.2) 

Po       Po  Po 

For  isentropic  flow,  the  free  stream  Mach  number  gives  the  density  ratio 
pr/p0  where  pr  =pTC  corrected  for  wind-tunnel  screen  losses.  Thus,  by  knowing 
the  fringe  number,  p/p0  can  be  determined.   A  negative  fringe  number  represents 
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an  accelerated  flow  and  vice  versa.  Since  the  Mach  number  distribution  is  known 
from  the  fringes,  the  pressure  distribution  can  be  calculated  in  the  flow  field  and 
on  the  airfoil  surface. 


C    = 


'p.; 


p    ymL^Poo    ) 

To  quantify  the  data  in  the  interferograms,  the  PDI  images  were  scanned, 
creating  enhanced  gray  scale  images  (300  dots  per  inch)  that  were  processed 
further  on  a  workstation.  The  three  registration  markers,  one  each  above  and 
below  the  airfoil  quarter  chord  point  and  one  below  the  leading  edge,  were  placed 
on  one  of  the  rotating  windows  to  give  proper  airfoil  orientation  and  scaling  of  the 
fringe  patterns,  as  shown  in  Chapter  V.  The  fringe  program  allowed  the  user  to 
pick  points  on  the  image  which  were  the  centers  of  the  dark  fringes  and  the  points 
where  the  fringes  intersected  with  the  airfoil  surface.  Starting  at  the  stagnation 
contour,  the  user  is  prompted  for  the  first  fringe  value  and  then  all  fringes  are 
mapped  over  the  lower  and  upper  surfaces  of  the  airfoil.  A  cursor  arrow  was  used 
to  select  the  fringe  ends  that  intersected  the  airfoil  surface,  and  by  using  the  middle 
mouse  button  the  point  was  mapped  to  the  nearest  pixel  on  the  airfoil  surface. 
Fringe  incrementing  could  be  cycled  as  desired  from  selected  values  of  1 ,  0  or  - 1 . 
Zoom,  back  up,  manual  incrementing,  and  curve  fitting  options  were  available 
during  data  processing  of  each  density  contour  or  fringe.  Each  fringe  was  only  2 
to  3  pixels  wide  near  the  leading  edge,  and  it  was  readily  apparent  after  curve 
fitting  (in  green  on  screen)  each  contour  if  the  center  of  the  fringe  was  properly 
mapped.  The  mapping  process  was  immediately  repeated  if  required  by  using  the 
back-up  mode.  When  the  upper  surface  mapping  was  completed,  the  right  mouse 
was  used  to  save  the  entire  image  data  in  a  QPLOT  data  format.  Several  PDI 
images  and  resultant  QPLOT  plots  are  part  of  the  presentation  and  discussion  in 
Chapters  V  and  VII. 

E.        MEASUREMENT  UNCERTAINTIES 

LDV  measurements  were  made  at  the  mid-span  plane  in  the  center  of  the 
tunnel  over  a  NACA  0012  airfoil  that  was  machined  to  0.001 -inch  accuracy. 
Relative  uncertainties  in  the  data  were  (1)  ±0.005  for  Mach  number,  (2)  ±0.05  for 
angle  of  attack,  and  (3)  ±0.005  for  reduced  frequency.   The  LDV  velocity  data 
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uncertainties  are  estimated  to  be  <5%  in  attached  flows.  The  z-component  of 
vorticity  was  obtained  by  interpolating  the  velocity  component  data  with  a  spline 
curve  fit  and  then  calculating  the  velocity  gradients  from  the  fitted  curve. 
Therefore,  the  uncertainty  of  the  vorticity  data  was  estimated  to  be  rather  high  at 
20%,  but  is  common  to  all  the  distributions  and  a  discussion  of  comparisons  was 
still  considered  valid.  The  fringe  uncertainty  in  the  PDI  images  was  estimated  to 
be  one  in  general,  but  increasing  to  three  near  the  airfoil  leading  edge  where 
fringes  converged  in  the  region  of  the  suction  peak.  Each  fringe  represented  a 
density  variation  of  Ap/p0  =0.0085,  and  the  deduced  uncertainty  in  Cp  from  the 
fringes  was  then  estimated  at  +0.2. 
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Previous  experimental  studies  in  the  facility  have  clearly  established  the 
two-dimensional  nature  of  flow  up  to  the  onset  of  dynamic  stall  and  even  in  the 
attached  region  of  partially  separated  flows  [Ref.  41  and  42].  In  addition,  previous 
tests  in  the  tunnel  attest  to  the  flow  quality  as  the  entrance  section  was  previously 
used  for  aeroacoustic  studies  [Ref.  43].  Also,  excellent  flow  uniformity  of  0.25% 
at  58  m/s  was  reported  in  the  tunnel,  in  addition  to  low  turbulence  intensity 
(0.083%)  with  a  bandwidth  of  50  to  50,000  Hz  [Ref.  43  and  44]. 

F.        DATA  VALIDATION  (LDV  SYSTEM) 

The  data-acquisition  and  processing  software  incorporated  (1)  checks  for 
detecting  oscillating  drive  frequency  variations  beyond  a  pre-set  tolerance  and  (2) 
phase  averaging  by  binning  the  data  appropriately  and  plotting  histograms  and 
phase  distributions  of  the  velocity  components.  A  minimum  of  10,000  samples 
was  collected  per  channel  at  each  measurement  point  for  the  oscillatory  cases 
along  with  encorder  information  for  each  sample.  For  the  limited  number  of 
steady  angle-of- attack  test  cases,  5,000  samples  were  collected  at  each  point. 
Anytime  the  data  were  widely  scattered,  the  entire  data  set  was  rejected  and  a  new 
set  of  data  was  acquired.  The  electronic  components  were  carefully  set  and 
maintained  for  the  large  dynamic  flow  range,  and  at  some  locations  30  minutes 
were  required  to  complete  data  acquisition. 

For  the  low-amplitude  test  case,  the  data  were  sorted  into  36  bins, 
corresponding  to  different  phase  angles  in  a  range  of  ±5°.  Also,  during  the  data- 

24 


analysis  stage,  the  raw  data  files  for  the  high-amplitude  test  case  were  processed 
further  by  sorting  the  data  into  120  bins  at  a  resolution  of  +1.5  °.  Whenever  the 
stipulated  number  of  coincident  samples  (50  in  this  case)  was  not  present  in  any 
bin,  that  bin  was  said  to  contain  a  "hole."  A  spline  curve  fit  was  used  to  fill  the 
"hole"  data  between  validated  data  bins,  which  worked  successfully  over  most 
phase  angles  [Ref.  34]. 

Histograms  of  raw  LDV  indicate  the  quality  of  the  recorded  data,  the 
distribution  of  recorded  velocities,  the  mean  velocity  value  normalized  by  UM,  and 
the  root  mean  square  (RMS)  intensity  or  standard  deviation  of  the  recorded 
samples.  The  initial  histogram  is  for  5,000  samples  at  a  =  8°  and  x/c  =  0.000  and 
y/c  =  0.083  in  steady  flow  for  the  U-component  velocity.  The  measured  height 
above  the  airfoil  surface  is  4.82%  chord.  Over  50%  of  the  data  samples  were 
recorded  at  1.1  U^  with  a  mean  velocity  of  1.07  UM  for  the  entire  distribution,  as 
shown  in  Figure  3.7.  The  standard  deviation  of  1 1%  is  high  for  the  attached  flow 
condition  and  measurements  outside  the  boundary  layer,  but  is  attributed  to  low- 
level  noise  in  the  system  as  noted  by  the  recorded  samples  around  0.0  UM  and 
0.511^,  which  are  well  outside  of  the  2  a  range.  In  such  cases,  the  data  were 
reprocessed,  excluding  these  samples  and  decreasing  the  intensity  considerably.  In 
this  way,  only  mean  values  were  considered. 
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Figure  3.7.  Histogram  of  Normalized  U-component  Velocity  at  x/c  =  0.000; 
y/c  =  0.083.  M00  =  0.3,  a  =  8°,  k  =  0.0,  Re  =  540,000. 
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For  the  second  example,  histograms  of  both  raw  LDV  velocity  components 
are  given  for  the  high-amplitude  oscillation  case  (a  =  10°  -10°  sin  cot)  forward  of 
the  airfoil  at  x/c  =  -0.133  and  y/c  =  0.083  in  Figure  3.8.  At  the  phase  angle  of  180 
degrees  ( a  =  10°  degrees  up),  the  number  of  coincident  data  samples  recorded  was 
125,  or  higher  than  the  83  samples,  if  evenly  distributed.  Near  50%  of  all  U- 
component  velocities  were  recorded  at  the  mean  velocity  of  1.31  U^.  Both 
velocity  component  data  distributions  give  appearance  of  the  desired  bell  curve. 
For  the  U-component,  data  points  at  0.9  U^  and  1.05  UM  are  outside  2  a  and  are 
considered  as  system  noise  that  contributes  to  the  8.0%  RMS  value.  The  V- 
component  samples  were  recorded  with  a  wider  data  distribution  around  the  mean 
velocity  of  0.38  U^ .  Velocities  at  0.15  U^  and  0.23  U^  can  also  be  considered  part 
of  system  noise. 
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Figure  3.8.   Histograms  of  Normalized  U-  and  V-component  Velocities  at  x/c  = 
-0.133,  y/c  =  0.083.  M.,  =  0.3,  a  =  10°  - 10°  sinco  t  =  10°(up),  k  =  0.1 . 

Normalized  velocities  of  both  the  U-  and  V-components  were  recorded 
throughout  the  oscillation  cycle  at  the  same  location  forward  of  the  airfoil,  as 
shown  in  Figure  3.9.  The  figure  legend  indicates  that  if  the  number  of  coincident 
particles  validated  in  both  components  was  less  than  20,  the  data  points  were 
denoted  by  a  "hole."  U-component  velocities  follow  the  motion  of  the  airfoil  up  to 
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the  dynamic  stall  angle  of  attack  at  a  =  15.9°.  The  U  velocity  increased  rapidly 
from  1.13  U^  to  1.32  UM  while  the  angle  of  attack  decreased  from  a  =  16°  to 
a  =  5°  down,  during  flow  reattachment.  The  V-component  velocities  increased 
very  rapidly  during  airfoil  upward  motion  before  reaching  the  dynamic  stall  angle 
of  attack.  Holes  in  the  data  were  shown  from  16  °  angle  of  attack  on  the  upstroke 
to  the  same  angle  of  attack  on  the  downstroke  because  the  airfoil  blocked  the 
recording  of  data  by  the  off-axis  forward  scatter  arrangement.  Data  "holes"  were 
interpolated  from  surrounding  validated  samples,  or  were  not  used  in  the 
subsequent  analysis.  When  holes  were  present  over  a  wide  phase  angle  range,  the 
data  were  not  used  at  all  for  the  results  given  in  Chapter  V.  Even  though  data 
"holes"  are  shown  in  this  Figure  for  <20  validated  samples,  the  data  were  not  used 
if  <50  samples  were  recorded  in  each  bin. 
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Figure  3.9.  Phase  Distribution  of  Normalized  U-  and  V-component  Velocities  at 
x/c  =  -0.133,  y/c  =  0.083.  M00=  0.3,  a  =  10°  -10°  sin  cot,  k  =  0.1. 

For  the  low-amplitude  forced-oscillations  (a  =  10°  -2°sincot),  example 
LDV  raw  data  histograms  in  Figure  3.10  show  a  wide  range  of  instantaneously 
recorded  velocities  at  the  closest  measurement  location  of  1.93%  chord  above  the 
airfoil  surface  at  x/c  =  0.033.  The  number  of  coincident  data  samples  (368)  was 
large.   The  mean  U-component  velocity  was  0.95  UM  at  10  °  angle  of  attack  (up), 
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and  instantaneous  velocities  ranged  from  0.56  <  U/U„  <  1.44  with  a  maximum  of 
7%  of  the  samples  recorded  at  0.85  U^.  Instantaneous  V-component  velocities 
ranged  from  slightly  less  than  -0.02  U^  to  0.6  UM  with  a  mean  average  velocity  of 
0.20  U^,  but  favoring  a  distribution  towards  0.3  UM.  The  range  of  measured 
velocities  may  be  explained  by  the  unsteady  nature  of  the  flow  within  the 
separation  bubble  as  indicated  in  the  bell  curve  distributions. 


8.00 


0.00 


(a)  U-component. 


(b)  V-component. 


Figure  3.10.   Histograms  of  Normalized  U-  and  V-component  Velocities  at  x/c  = 
0.033;  y/c  =  0.083.  MM=0.3,  cc  =  10°  -2°  sin  cot  =  10°  (up). 

The  maximum  and  minimum  U-component  velocities  in  Figure  3.11 
occurred  at  totally  different  angles  of  attack  when  compared  with  the  distribution 
given  at  the  airfoil  leading  edge  (not  shown).  U-component  velocities  increased 
rapidly  after  the  light  dynamic  stall  at  1 1  °  angle  of  attack  (down)  to  a  maximum 
value  of  1.27  UM  at  the  bottom  of  the  oscillation  cycle.  The  U  velocities  decreased 
from  8  to  1 1  °  angle  of  attack  (up),  which  was  the  opposite  of  what  was  expected 
because  of  airfoil  rotation,  and  is  due  to  the  presence  of  the  separation  bubble  at 
the  measurement  point.  Validated  U  and  V-component  velocities  were  sparse 
between  1 1.5  °  angle  of  attack  (up)  and  1 1.6°  angle  of  attack  (down)  because  of 
the  measurement  location  in  the  well  developed  separation  bubble  region.  V- 
component  velocities  were  low  at  0.2  U^  and  appeared  to  be  fairly  constant 
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throughout  the  oscillation  cycle,  except  near  the  top  of  the  oscillation  cycle.  The 
velocities  were  50%  less  than  those  at  the  airfoil  leading  edge.  The  V-component 
data  distribution  with  a  reduced  validated  sample  criteria  (5)  shows  that  "holes" 
were  still  present  for  angles  of  attack  approaching  the  top  of  the  oscillation  cycle 
and  shortly  thereafter  during  the  rapid  occurrence  of  light  dynamic  stall.  In  this 
case,  the  data  from  11.5  °  (up)  to  1 1.6°  angle  of  attack  on  the  downstroke  (region 
of  partial  airfoil  blockage)  were  curve  fitted  from  validated  bins  (>50)  of  adjoining 
phase  angles. 
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Figure  3.11.  Phase  Distribution  of  Normalized  U-  and  V-component  Velocities  at 
x/c  =  0.033,  y/c  =  0.083.  Moo=0.3,  cc  =  10°  -2°  sin  cot,  k  =  0.1. 
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IV.       DESCRIPTION  OF  COMPUTATIONAL  ANALYSIS  METHODS 

Potential  flow,  boundary  layer,  and  N-S  analysis  methods  were  used  to 
assess  the  prediction  capabilities  of  these  methods  by  comparing  the  computational 
results  with  the  measured  data.  Potential  flow  analysis  is  generally  regarded  as 
sufficient  to  predict  an  airfoil's  low  incidence  lift  and  moment  characteristics  at 
reasonably  large  Reynolds  numbers  because  the  flow  is  attached  and  the  boundary 
layers  on  the  airfoil  upper  and  lower  surfaces  are  thin.  The  skin  friction  and  drag 
characteristics  can  be  obtained  by  coupling  the  potential  flow  analysis  with  a 
solution  of  the  boundary  layer  equations.  In  the  classical  approach,  the  computed 
potential  flow  pressure  distribution  is  used  as  input  into  the  boundary  layer 
equations.  This  approach  limits  the  analysis  to  attached  flow  problems  only 
because  the  boundary  layer  equations  become  singular  at  the  point  of  flow 
separation.  In  recent  years,  this  limitation  could  be  overcome  by  starting  with  a 
reasonable  estimate  for  the  boundary  layer  displacement  thickness  distribution, 
computing  the  pressure  distribution  from  the  potential  flow  analysis,  and 
continuing  this  viscous-inviscid  interaction  process  to  convergence.  These 
viscous-inviscid  interaction  methods  permit  the  computation  of  separation  bubbles 
and  of  airfoil  flows  with  significant  regions  of  trailing-edge  stall.  However,  a 
more  general  and  powerful  approach  has  become  available  with  the  introduction  of 
large-memory  high-speed  computers,  which  permit  the  direct  numerical  solution 
of  the  Reynolds-averaged  N-S  equations  and  enable  the  computation  of  steady  and 
unsteady  attached  and  separated  airfoil  flows  with  one  and  the  same  method. 
Naturally,  this  greater  capability  comes  at  the  price  of  significantly  larger 
computation  times.  In  this  thesis,  the  author  therefore  uses  only  the  potential  flow, 
the  classical  boundary  layer,  and  the  N-S  methods  for  the  analysis  of  the  NACA 
0012  airfoil  experiments  described  in  Chapter  V. 

A.        UNSTEADY  POTENTIAL  FLOW  CODE  (UPOT) 

1.         Introduction 

Potential  flow  solutions  for  steady  or  unsteady  inviscid  and  incompressible 
airfoil  flow  are  provided  by  the  UPOT  code.  The  conservation  equations  in  this 
case  reduce  to  the  Laplace  equation,  which  can  be  solved  by  distributing  sources 
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and  vortices  on  the  airfoil  surface.  Whenever  the  airfoil  experiences  an  angle  of 
attack  change,  a  vortex  is  shed  from  the  trailing  edge.  The  wake  vortices,  in  turn, 
influence  the  flow  over  the  airfoil,  At  each  time  step,  airfoil  pressure  distribution, 
force  and  moment  coefficients  are  calculated. 

2.         Theory  and  Numerical  Method 

The  UPOT  code  is  capable  of  computing  the  unsteady  flow  over  a  two- 
dimensional  (2-D)  airfoil  based  upon  the  panel  method  developed  by  Hess  and 
Smith  [Ref.  45].  In  this  approach,  the  flow  is  modeled  by  using  uniform  source 
distributions  (qp,  which  are  allowed  to  vary  from  panel  to  panel,  and  uniform 
vorticity  distributions  (y),  which  are  held  constant  on  all  the  panels.  These 
distributions  automatically  satisfy  the  Laplace  equation.  Use  of  the  superposition 
principle  allows  the  buildup  of  complicated  flow  fields  by  using  simple  source  and 
vortex  flows,  if  an  appropriate  boundary  condition  can  be  satisfied  at  the  airfoil 
surface  boundary.  Thus,  the  total  velocity  potential  is  constructed  by 
superimposing  the  velocity  potential  for  the  free-stream  flow,  the  source  flows, 
and  the  vortex  flows.  Detailed  descriptions  of  all  equations  are  given  in  Reference 
46,  and  only  primary  equations  are  listed  here.  The  total  velocity  potential  is  given 
by 

<D  =  V00(xcosa  +  ysina)+I  jpanelj    -^lnr--^-e   ds  (4.1) 


j=i 


^2tt  2je    j 


Once  the  velocity  potential  is  calculated  for  the  steady  flow  problem,  the 
velocity  is  evaluated  by  taking  the  grad(  <E> ),  and  the  pressure  is  determined  from 
the  steady  Bernoulli  equation.  The  unsteady  problem  is  discussed  later. 

The  boundary  conditions  consist  of  the  airfoil  surface-flow  tangency 
condition  and  the  Kutta  condition,  which  requires  that  the  pressures  be  equal  at  the 
upper  and  lower  trailing  edge  panels  of  the  airfoil.  This  leads  to  a  system  of  linear 
simultaneous  equations  for  the  unknown  source  and  vortex  strengths  on  the  panels 
and  to  the  influence  coefficient  concept.  An  influence  coefficient  is  defined  as  the 
velocity  induced  at  any  point  in  the  flow  field  (field  point)  by  a  unit  strength 
singularity  (source  or  vortex)  placed  anywhere  in  the  flow  field.  A  detailed 
derivation  of  influence  coefficients  is  found  in  Reference  46.     Finally,  the 
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equations  are  written  in  matrix  form  and  solved  by  the  Gaussian  elimination 
method. 


Li,i 


_an+l,l 


al,n+l 

"qi" 

'  b,  ' 

• 

q2 

b2 
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qn 

bn 

an+l,n+l_ 

_y  . 

.bn+l_ 

(4.2) 


The  unsteady  problem  is  formulated  by  describing  the  vortex-shedding 
process,  shown  in  Figure  4.1,  The  vortex-shedding  process  creates  non-linearities, 
because  the  wake  vortices  influence  the  flow  over  the  airfoil,  which  in  turn  alters 
the  vortex-shedding  process  as  the  airfoil  maneuvers  in  time.  The  disturbance 
potential  is  approximated  by  using  backward  finite  differences  to  integrate  the 
velocity  field  in  two  stages  from  upstream  free-stream  conditions  to  the  leading 
edge  of  the  airfoil  and  along  the  airfoil  surface  to  each  panel  control  point.  The 
rate  of  change  of  the  disturbance  potential  must  be  evaluated  at  each  time  step  for 
all  panel  control  points  and  now  includes  the  contributions  from  the  shed  vorticity 
panel  and  the  wake  core  vortices. 


The  total  circulation  in  a  potential  flow  must  be  preserved  according  to  the 
Helmholtz  vorticity  conservation  theorem.  The  vortex-shedding  process  requires 
that  any  change  in  circulation  around  the  airfoil  must  be  countered  by  a  change  in 
vorticity  in  the  wake  of  opposite  sign,  but  equal  in  magnitude.  Thus,  the  vorticity 
in  the  wake  element  (shed  vorticity  panel)  of  length  Ak  is  equal  to  the  negative 

change  in  circulation  around  the  airfoil  for  the  previous  time  step.  The  overall 
circulation  Tk  at  time-step  tk  is  simply  yk  multiplied  by  the  airfoil  perimeter,  £. 

The  shed  vorticity  panel  provides  the  desired  communication  to  carry  the  solution 
from  one  time  step  to  the  next,  where  the  vorticity  conservation  condition  has  to  be 
satisfied  at  every  time  step. 


MYw)k  =  rk_,-rk=  ^(Yk-i-Yk) 


(4.3) 


The  airfoil  trailing-edge  wake  panel  with  unknown  shed  vorticity  strength, 
unknown  length,  and  unknown  trailing-edge  angle  introduces  three  additional 
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unknowns  to  the  system.  In  addition,  N  unknown  source  strengths  and  one 
unknown  vorticity  strength  on  the  airfoil  are  all  time  dependent  and  must  be 
solved  by  iteration  techniques  to  a  final  solution. 


Vortex  Shedding  at  Time  Step,  tk 


Helmholtz's  Theorem: 

Ak(Tw)k  +  rk  =  rk_, 

Panel  j 

Source  Distribution  (qj  )k 

Vorticity  Distribution  yk 


(Tk-3-rk_2) 

^  ft 

(rk_4_rk_3) 

A, 


CVJk 


(UJ.         tan0k  = 


Ak  =(tk  -  tk_, )  ^[(Uw)k]2+[(Vw)k[ 
(Vw)k 


(Uw) 


w^lc 


Figure  4.1 .      Extension  of  Panel  Methods  Representation  for  Unsteady  Flow 
[Ref.45]. 

The  Kutta  condition  must  include  the  rate  of  change  of  the  potential  at  the 
trailing-edge  panel,  which  is  directly  related  to  the  rate  of  change  of  the  total 
circulation. 
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Two  additional  equations  are  required  to  solve  for  the  (N  +  4)  unknowns  in 
the  unsteady  problem,  and  assumptions  are  required  for  the  geometry  of  the  wake 
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panel  as  suggested  by  Basu  and  Hancock  [Ref.  47].  First,  the  wake  panel  is 
aligned  with  the  direction  of  the  local  resultant  velocity  at  its  midpoint  in  the 
moving  airfoil  reference  frame.  In  addition,  the  wake  panel  length  is  assumed  to 
be  proportional  to  the  magnitude  of  the  local  resultant  velocity  at  the  panel 
midpoint  and  the  size  of  the  time  step.  These  boundary  conditions  are  written  as 
follows: 

tan0k=^^  (4.5) 

(Uw)k 

Ak=(tk-tk_1)^[(Uw)kf  +  [(Vw)k]2  (4.6) 

The  unsteady  flow  Bernoulli  equation  for  the  pressure  coefficients  must  be 
written  with  respect  to  the  airfoil-fixed  system.  Giesing  [Ref.  48]  showed  the 
equation  to  be  written  as 

f_  "q°      _      *  stream       _      *  total        _       2      Oty 

(l/2)pVoc2  VB  I    VM  J        V^2  9t 


C     —  —         stream        _      T  total        _        ■*-      vyr  /a  n\ 


For  this  work,  2-D  aerodynamic  coefficients  of  lift,  drag,  and  pitching 

moment  are  computed  for  a  NACA  0012  airfoil  at  steady  angles  of  attack  and 

during  sinusoidal  oscillations  about  the  quarter-chord  point  by  integrating  pressure 
distributions,  assuming  a  constant  Cp  on  each  panel. 

B.        STEADY  BOUNDARY  LAYER  FLOW  CODE 

1.         Introduction 

The  boundary-layer  equations  are  incorporated  into  a  program  (BL2D) 
developed  by  Cebeci  et  al.  [Ref.  18].  The  equations  and  analysis  are  restricted  to 
2-D  steady,  incompressible,  and  viscous  flow.  Inputs  to  the  program  include  the 
Reynolds  number  and  estimated  or  prescribed  transition  locations,  as  well  as  panel 
coordinate  and  velocity  information  computed  by  a  panel  method.  Even  though 
the  flow  accelerates  over  the  airfoil  upper  surface  to  Mach  numbers  much  greater 
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than  the  free-stream  value  of  0.3  and  is  considered  compressible,  the 
incompressible  BL2D  code  with  a  transition  model  incorporated  is  used  to 
evaluate  the  extent  of  its  capability  to  predict  the  flow  physics  evident  in  the 
measurements.  Output  is  generated  for  several  boundary-layer  parameters, 
including  velocity  profiles,  skin  friction  coefficients,  and  displacement  thickness. 

2.         Theory  and  Models 

The  2-D,  incompressible,  steady  boundary-layer  equations  for  both  laminar 
and  turbulent  flows  are 

dx     dy 


3u        3u         1  dp 
u^  +  v  — -  = -  +  u 


fd2u]      d(uV) 


dx        dy        p  dx        \dy2 
dp 


(4.9) 


3y 


3y 
=  0  (4.10) 


where  u  =  (i/p  is  the  kinematic  viscosity,  [i  is  the  coefficient  of  viscosity,  and  u', 
v'  represent  the  turbulent  velocity  fluctuations.  These  equations  must  be  solved 
subject  to  the  boundary  condition  on  the  airfoil  surface  ( u  =  0  and  v  =  0),  and  at 
the  outer  edge  of  the  boundary  layer,  u  =  ue(x). 

To  use  Equation  4.9,  an  expression  must  be  found  for  the  Reynolds  shear 
stress  term.  Since  it  is  not  feasible  to  attempt  calculating  the  actual  value, 
empirical  models  must  be  used.  One  such  model  is  the  eddy-viscosity  concept: 

9u 


puV  =  pem^  (4.11) 

3y 

where  em  is  an  empirical  term  called  the  turbulent  eddy  viscosity.  A  refinement  of 
the  eddy-viscosity  concept  called  the  Cebeci-Smith  (CS)  [Ref.  13]  eddy-viscosity 
model  is  used  in  the  boundary  layer  code  for  this  work.  The  CS  model  divides  the 
viscous  region  into  an  inner  layer  and  an  outer  layer,  as  described  in  Reference  17. 
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The  Chen-Thyson  [Ref.  20]  transition  length  model  is  incorporated  into  the 
boundary  layer  code  for  this  study,  giving  the  intermittency  [  ytr(x)  ]  as 


Ytr(x)  =  l-exp 


vjY(r  u  ue  j 


(4.12) 


and  Falkner-Skan  transformation  variables,  £  and  r)  are  used  for  the  solutions. 
The  intermittency  factor,  ytr(x),  models  the  length  of  the  transition  region,  as 

explained  further  in  the  N-S  code  methodology  description.  Using  Equation  4.1 1, 
the  momentum  equation  for  a  turbulent  flow  (Equation  4.9)  can  be  rewritten  in  the 
same  form  as  for  a  laminar  flow: 


du       du        1  3p      d 
u  —  +  v —  = -  + 


'&? 


3x        dy        p  3x     3y  I     dy 


(4.13) 


where  b  =  \)  +  em.   Thus,  the  only  computational  difference  between  laminar  and 
turbulent  boundary  layers  is  the  addition  of  the  turbulent  eddy  viscosity. 

The  determination  of  the  transition  location  from  laminar  to  turbulent  flow 
is  sometimes  one  of  the  most  critical  factors  in  the  success  of  computational 
efforts  to  predict  or  reproduce  physical  phenomena.  Assuming  the  entire  flow  to 
be  turbulent  may  be  a  reasonable  approximation  for  some  applications,  especially 
when  the  Reynolds  number  is  high.  In  contrast,  there  are  important  flows  that 
require  more  accurate  predictions,  especially  flows  in  the  transitional  Reynolds 
number  regime. 

Until  the  transition  mechanism  and  the  many  factors  that  can  affect  it  are 
understood  more  fully,  we  must  rely  on  the  traditional  engineering  approach  of 
modeling  the  start  and  range  of  the  transitional  flow  region.  In  addition, 
experimental  values  have  an  expected  range  of  accuracy,  and  empirically 
calculated  values  may  deviate  from  the  actual  onset  of  transition.  Therefore,  to 
provide  an  initial  estimate  for  the  transition  location  when  no  other  method  of 
determination  was  available,  Michel's  criterion  [Ref.  1 3]  was  incorporated  into  the 
code: 
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Re    =1.174 


' \     22,400^ 

1  + 

Rex 

V  Atr     / 


Rex 
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where  Re  =ue0/v  is  the  Reynolds  number  based  on  mom  atum  thickness  at 
transition  and  Rx    is  the  Reynolds  number  based  on  the  transition  onset  location. 

Otherwise,  experimental  data  or  input  values  are  used  to  specify  the  onset  of 
transition. 

Then,  the  transitional  flow  region  is  calculated  by  using  the  Chen-Thyson 
model  [Ref.  20],  shown  as  ytr(x)  in  Equation  4.12.   The  flow  transition  range  is 

neither  fully  laminar  nor  fully  turbulent.  It  is  a  region  of  intermittency,  in  which 
turbulent  spots  gradually  appear  with  progression  in  the  streamwise  direction.  Use 
of  the  transition  model  in  this  work  will  be  shown  to  be  essential  for  obtaining 
results  in  agreement  with  the  measurements  for  low  Reynolds  number  flows  and 
will  be  discussed  more  fully  in  the  section  describing  the  N-S  code  methodology. 

3.         Numerical  Implementation 

The  x  /  c  and  y  /  c  coordinates  of  the  airfoil  are  supplied  as  part  of  the  input 
to  the  program  called  xc  and  yc.  The  program  redefines  these  coordinates  into  a 
single  parameter  corresponding  to  an  airfoil  surface  distance  from  the  stagnation 
point: 


xi  =  xi-i  +  V(xci  -  xci-i  )2 + (yci  -  yci-i  )2 


The  variable  x}  used  internally  by  the  program  in  the  boundary  layer  equations  is 
the  surface  coordinate. 

The  grid  variables  are  transformed  with  x  simply  scaled  by  a  reference 
length,  the  chord  of  the  airfoil,  so  that  £  =  x/c .  Since  the  surface  distance  x  is 
used,  these  will  not  be  the  same  as  the  x/c  coordinates  that  are  input.  The 
dimensionless  stream  function  f(x,T|)  is  defined  by 

f(X,lD  =  -JL-r 

Vue"ux 
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For  the  y  variable,  the  well  known  Falkner-Skan  transformation  is  used: 


V\)x 


where  y  is  the  normal  coordinate  along  which  the  thickness  of  the  boundary  layer 
is  measured.  The  dimensionless  similarity  variable  r\  eliminates  the  growth  of  the 
boundary  layer  in  laminar  flow  and  reduces  it  in  turbulent  flow.  This  enables 
larger  steps  in  the  stream  wise  direction  and  improves  computational  efficiency. 

Using  the  new  variables,  Equation  4.13  and  the  boundary  conditions  may 
be  rewritten  as 


4r—-f"— 


(bfT  +  ^-^ff"  +  m[l-(f')2]  =  £ 


(4.14) 


r,  =  0  f'  =  0  ffeO)  =  fw(£)  =  -      ^      J^Ld^ 


n  =  r|e 


f  =  1 


where  r|e  is  the  transformed  boundary  layer  thickness  corresponding  to  5,  and  RL 
is  the  Reynolds  number  based  on  the  reference  velocity  ( uM)  and  the  reference 
length,  c,  the  chord  of  the  airfoil.  The  prime  denotes  differentiation  with  respect 
to  tj  .  The  dimensionless  pressure  gradient  parameter  m  is  defined  by 

ue/Uoc       d£ 
The  velocity  components  u  and  v  are  related  to  the  dimensionless  stream  function 


by 


u  =  uef 


V  =  -^/ueUX 


f      d    i df     ~,dr\ 

T=^— >ex+— +r— ^ 

-^/UgX  dx  dx        dx 


Equation  4.14,  a  second-order  partial  differential  equation  can  be  solved  by 
various  numerical  methods,  such  as  the  Crank-Nicholson  or  Keller  Box  methods. 
The  Keller  Box  method  [Ref.  18]  is  used  in  this  work  and  proves  to  be  efficient  for 
the  boundary  layer  calculations,  but  first  requires  reformulating  the  higher-order 
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equations  into  a  set  of  first-order  equations.  At  each  rectangular  grid  section,  these 
equations  are  approximated  by  using  centered-difference  derivatives,  averaging 
values  at  the  center  of  the  "box,"  as  shown  in  Figure  4.2. 

The  higher-order  equations  are  formulated  into  a  set  of  first-order  equations 
with  a  second-order  truncation  error  by  using  the  following  definitions. 


f'  =  u 


u   =  V 


Equation  4.14  may  now  be  expressed  as  a  first-order  system: 

.,      ,     m  +  1  r        21     J    du        3f> 

(bv)  +_fv+m[1_u2]=^u__v_ 

T|  =  0  U  =  0  f  =  fw(x),  T|  =  *ne  U  =  l 


(4.15) 


(4.16) 


Note  that  the  u  and  v  in  Equations  4.15  and  4.16  are  not  the  velocity  components. 
They  are  two  new,  arbitrarily  selected  variable  names  for  the  expression  of  the 
first-order  system.  The  program  was  coded  using  these  particular  variable  names, 
and  the  present  numerical  discussion  uses  them  for  consistency. 


*n-i    ^n-'/,^n 


Figure  4.2.  Net  Rectangle  for  Difference  Approximations  [Ref.  18]. 

40 


Referring  to  Figure  4.2,  the  grid  points  may  be  described  as 
£n=£n-1  +  kn,        n  =  l,2,....,N  £,°  =  0 

Tijsiij.j+hj,       j  =  l,2,....,J  r|0=0      Tij=Tie 

where  the  capital  letters  N  and  J  are  the  maximum  number  of  grid  points  used  in 
the  streamwise  and  normal  directions,  respectively.  The  superscript  n  is  not  an 
exponent,  but  a  counter  just  like  the  subscript  j .  This  upper  and  lower  notation 
allows  using  both  counters  on  the  same  variable  when  needed. 

Considering  one  box  of  the  grid,  the  finite  difference  approximations  may 
be  written  for  the  midpoint  of  the  right  side  segment,  P2P4,  using  centered 
difference  derivatives: 


fn_fn  nn  4-nn 

-1J—  — ~2  U^ 


(4.17) 


Equation  4.16  may  be  approximated  in  the  same  manner  for  the  midpoint  of  the 
box  by  centering  first  in  one  direction  then  the  other.  The  resulting  finite 
difference  equation  is 

i-(bjvj  -b%.lvaH)  +  al(fv)l.y2  -a2(u2)?_1/2  +an(v^2f5L1/2  - f J-V?-,/2)  =  RjLl/2 

(4.18) 


where 


Tn-1     _ 

Lj-V2  - 


Rp/2  =  -L-J/2  +  an[(fv)j42  - (u2)^] -  mn 
±(b^^by_1vJ.1)  +  ^(fv)j.1/2  +  m(l-(u2)H/2) 


n-l 


an=i 


n-1/2 


k" 


«1  = 


mn+l 


+  a' 


a-,  =mn  +an 


The  boundary  conditions  are  rewritten  at  £,  =  %n  as 

41 


fg  =  fw  uS=0  Ujn=l 

The  resulting  implicit,  nonlinear  difference  equations  are  linearized  by 
Newton's  method  and  solved  by  a  block  elimination  method.  Equations  4.18 
comprise  a  set  of  3 J  +  3  equations  and  3 J  +  3  unknowns  (f  ",  u",  v",  where  j  =  0, 

1,  2,  ...,  J),  with  f-1-1,  u"-1,  and  v"-1  known.  The  method  assumes  that  an 
approximate  solution  is  known,  either  from  the  preceding  iteration  cycle  or  from 
the  previous  streamwise  station.  Then  small  unknown  quantities  are  added  to  the 
approximate  solution.  Using  the  arbitrary  iteration  variable  i  (superscript  n 
omitted  for  clarity),  the  approximate  equations  become 

fj+1  =fj  +  6fj  u*+1  =uj  +  du)  v)+]  =  v)+bv) 

with  i  =  0  corresponding  to  known  values  at  the  previous  streamwise  station, 
(£n_1),  these  expressions  may  be  substituted  into  Equations  4.17  and  4.18  for  the 
unknowns.  After  dropping  higher-order  terms  of  6,  a  linear  system  of  equations 
results: 

5fi-6fj_1-^-(5uj+8ui_1)  =  (r])j 

(s1)j5vi  +  (s2)j5v;_1+(s3)j8fj+(s4)j5fi_1+(s5)j5u]+(s6)j5u]_1=(r2)j 

5u5-8«S.1-^(6v5+8vi_1)  =  (r3)j 


where  the  right-hand  sides  are 


(ri)j  =  fj-i-fj  +  Vj-i/2 


n-l 


(r2)j  -Rj-l/2  ~ 


-J-(bjvj  -  bj_, vj_,)  +  a, (fv);_1/2  -  a2(u2);_1/2  +  an( v^2  fj_1/2 -f^2vj_1/2) 
j 


(r3)J  =  u;_1-u;  +  hjv;_1/2 
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and  the  coefficients  are 

(S2)J=l!:bJ-|+ffH-ffH'2 

Cs4)j-yij-i-— Vj_1/2 

(s5)j  =-a2u- 

(s6)j=-a2uj-i 
The  boundary  conditions  become 

5fj)  =  0  8^=0  8uj=0 

The  equations  may  be  arranged  into  a  block  tri-diagonal  structure  in  matrix - 
vector  form.  The  efficient  block  elimination  method  is  used  to  solve  for  the  small 
8  quantities.  Iterations  of  Newton's  method  continue  until  quantities  are  small 
enough  to  be  neglected. 

Velocity  profiles  are  output  by  the  boundary  layer  program.  At  each  station 
along  the  airfoil  surface,  the  program  calculates  the  u  velocities  for  each  value  of 
r\  in  the  normal  direction  of  the  grid.  To  retrieve  the  physical  distance  normal  to 
the  airfoil  surface,  an  inverse  transformation  is  required: 

y  =  -f4=  (4.19) 

>e/ux 

The  kinematic  viscosity  appears  only  indirectly  in  the  non-dimensional 

form  of  the  Reynolds  number.  Therefore,  the  following  equation  is  used  to  give 

the  actual  y  value: 

y  =  TlVx/RLue 
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where  RL  =  puoec/u.,  or  equivalently,  u^c/d,  and  x  and  ue  are  used  as  the  non- 
dimensional  forms  of  x/c  and  Ue/u^.  The  value  of  x  is  the  airfoil  surface 
distance. 

Plotting  the  shape  of  the  actual  velocity  profile  at  a  given  station  requires 
the  station  number,  the  corresponding  x  value,  the  u  velocities,  and  the 
corresponding  y  values.  All  values  are  provided  by  the  original  program  or 
Equation  4.19.  To  visualize  the  growth  of  the  boundary  layer,  the  height  of  the 
boundary  layer  is  determined  by  finding  where  the  u  velocity  has  reached  the  local 
external  velocity.  The  edge  of  the  boundary  layer  is  noted  when  u  reaches  0.995 
of  the  local  external  velocity  ( ue ). 

This  steady  boundary  layer  solution  can  also  used  for  the  computation  of 
the  boundary  layer  characteristics  on  the  harmonically  pitching  NACA  0012  airfoil 
because  the  pitch  rates  used  in  the  experiment  were  small  and  the  boundary  layer 
lag  effects  are  known  to  be  small  for  small  pitch  rates. 

C        NAVIER-STOKES  (N-S)  CODE  METHODOLOGY 

1.         Introduction 

A  numerical  solution  based  on  the  compressible  thin-layer  N-S 
approximation,  incorporating  an  empirical  transition  model,  is  used  to  study 
compressible  steady  and  unsteady  flows  over  an  airfoil  at  transitional  Reynolds 
numbers.  N-S  methods  require  computational  meshes  with  sufficient  grid 
clustering  near  the  wall  for  accurate  predictions  of  the  large-flow  gradients  in  the 
viscous  layers.  The  grid  outer  boundaries  must  be  placed  far  enough  (15  to  20 
chord  lengths)  from  the  body  of  interest  to  apply  the  inflow  and  outflow  boundary 
conditions.  Unsteady  airfoil  flow  solutions  were  computed  in  the  past  by  Rumsey 
and  Anderson  [Ref.  49],  Visbal  [Ref.  50],  and  Ekaterinaris  [Ref.  51]  with  single- 
grid  approaches,  but  none  accounted  for  transitional  flow  effects.  The 
development  of  leading-edge  separated  flows  are  investigated  by  discretizing  the 
equations  with  an  upwind-biased,  factorized,  iterative  scheme.  Solutions  obtained 
for  two-  and  three-dimensional  complex  flow  fields  by  Rai  [Ref.  52]  have  proven 
the  accuracy  of  upwind  methods. 
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2.         Numerical  Method 

For  this  study,  the  strong  conservation-law  form  of  the  governing  thin- 
layer,  compressible,  Reynolds-averaged,  N-S  equations  for  a  body-fitted 
coordinate  system  (£,  rj)  along  the  axial  and  normal  direction,  respectively,  is  as 
follows: 


-l 


dtq  +  d^F  +  dT]G=  Re-'a^S 


(4.20) 


The  Euler  equations  are  obtained  when  the  right-hand  side  is  set  equal  to 
zero.  In  Equation  4.20,  q  is  the  conservative  variable  vector,  F  and  G  are  the 
inviscid  flux  vectors,  and  S  represents  the  viscous  terms  in  the  normal  direction. 
The  conservative  variable  vector  q  is 


q  = 


p 

pu 

pw 

e 


N 


where  J  =  1  /  (xtz^  -  x^z* )  =  ^xi\z  -  ^zr|x  is  the  Jacobian  of  the  transformation.    F 
and  G  are  the  inviscid  fluxes  given  by 


F  =  i 

J 


pU 
PuU  +  £xp 

pwU  +  £zp 
(e  +  p)U-£tp 


and 


1 
G  =  - 

J 


pW 

pllW  +  T|xp 
pwW  +  Tlzp 

(e  +  p)W-n.tp 


The  viscous  term,  when  the  thin-layer  approximation  is  applied  and  also 

2li 
assuming  that  X  =  — -,  is  given  by 


J 


0 


lim]^^  +  ([i/3)m2r\z 
unijnvj  +  (|i/3)m2m4 
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where 


m]  =t\x2+t\z2, 


m2=rlxUTi+TlzWTr 


m3  =  —  (u2  +  w2)/2  +  KPr" 
arj 


^a2^ 


vaTly 


m4=(rixu  +  r|zw), 


a  is  the  local  speed  of  sound,  k  is  the  thermal  conductivity,  Pr  is  the  Prandtl 
number,  and  U,  W  are  the  contra- variant  velocity  components  given  by 


U  =  u£x+w^z  +  £t, 


W  =  ur|x  +  wtiz  +  r|t 


In  Equation  4.20,  all  geometrical  dimensions  are  normalized  with  the  airfoil 
root-chord  length,  c ;  the  density  is  normalized  with  the  free-stream  density,  p^ ; 

the  velocity  components  u  and  w  are  normalized  with  the  free-stream  speed  of 
sound,  aTO;  the  total  energy  per  unit  volume  is  normalized  with  p^ai;  and  the 
pressure  is  related  to  density  and  total  energy  through  the  equation  of  state  for  an 

t  i   -  2  .  ...2> 


ideal  gas,  p  =  (y  -  l)[e  - 1  -  |p(u   +  w  )]. 


3.         Numerical  Implementation 

The  numerical  integration  used  to  compute  the  mean  flow  is  performed  by 
using  a  third-order,  upwind-biased,  factorized,  iterative,  implicit  numerical  scheme 
given  by 


[i  +  h^( V^Atk  +  4A"k )]P  x  [i  +  n^V^k  +  Vf,B-k  -  Re-^M^ )]F 


x 


(4-21) 


(QE-QEk)  =  -[(Qf,k-Qu)+^(FP+1/2jk-FP.1/2>k)+ 

hri^i,k+l/2  _^i,k-l/2)~^e    hr|(Si,k  +  l/2      ^i,k— 1/2 )1 

where   h*  =  At/A£,  etc.,  A1  =  (df/dQ),  etc.,  are  the  flux  Jacobian  matrices,  A,  V  , 

and  8  are  forward,  backward,  and  central  difference  operators,  respectively.  The 
quantities  Fi+1/2k,  Gik+1/2,and  Sik+1/2  are  numerical  fluxes. 
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Time  accuracy  of  the  implicit  numerical  solution  is  obtained  by  performing 
Newton  iterations  to  convergence  within  each  time  step.  The  approximation  to 
Qn+I  at  each  subiteration  is  the  quantity  Qp.  When  p>2,  during  a  given 
subiteration,  Qp  =  Qn+1 ,  but  when  p  =  1  and  no  subiterations  are  performed,  then 
qp  =  Qn  and  Qp+1  =  Qn+1 .  By  subiterating  to  convergence,  linearization  and 
factorization  errors  are  minimized,  because  the  left-hand  side  of  Equation  4.21  can 
be  driven  to  zero  at  each  time  step.  The  inviscid  fluxes  F  and  G  are  evaluated 
using  Osher's  [Ref.  53]  upwinding  scheme. 


The  numerical  fluxes  for  a  third-order  accurate  upwind-biased  scheme  are 
given  by 


Fi+l/2.k  -Fi+l/2,k  +V6 


AF+        +2AF+ 

i-l/2.k  l  +  l/2.k 


-1/6 


AF"         +  2AF" 

i-?/2.k  1  +  1/2, k 


F(Qi.k^Qi+u)  +  V6[AF+(Qi+u,Qi>k)  +  2AF+(Qu,Qi+u)]- 
l/6[AF-(Qu,Qi+u)  +  2AF-(Qi+Kk,Qi,k)]; 


here  F  is  the  first-order  accurate  numerical  flux  for  the  Osher's  scheme  given  by 


Fi+1/2,k  =:  V2  Fiik  +  Fi+U  -  |§'*i  {F;  -  F"}dQ 


where  Fq  =  Fq  +  Fq ,  F*  =  (3F/8Q)1 ,  and  AF1  are  the  corrections  to  obtain  high- 
order  accuracy.  The  Osher  scheme  evaluates  the  flux,  assuming  a  shock  tube 
solution  where  Fq  is  piece-wise  continuous,  and  yields  good  predictions  of  the 

flux.  For  the  linearization  of  the  left-hand  side  of  Equation  4.21,  the  flux  Jacobian 
matrices  A,  B  are  evaluated  by  the  Steger- Warming  [Ref.  54]  flux- vector  splitting 
method.  The  linearization  errors  are  reduced  by  subiteration  to  convergence. 
Typically,  two  to  three  subiterations  are  sufficient  to  drop  the  residuals  two  orders 
of  magnitude  during  the  Newton  iteration  process.  Accurate  steady-state  solutions 
can  be  obtained  even  without  subiteration.  Two  subiterations  are  used  for  the 
unsteady  solutions  presented  in  this  thesis. 
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The  viscous  fluxes  Ss k+]/2  are  computed  with  central  differences  as 
follows: 

Sj,k+l/2  ~  S  |Qi.k+l/2'(Qr|)i.k+l/2'rli,k+l/2j' 

Qi.k+i/2  =  V2(Qi.k  +  Qi+u).  (QuWi/2  =  Qu+i  -  Qu- 

The  experimental  Reynolds  number  based  on  the  chord  length  for  the  test  case 
examined  is  540,000,  where  it  has  been  shown  that  transitional  flow  effects  should 
be  considered  in  the  region  near  the  airfoil  leading  edge. 

4.         Transition  Modeling 

For  flow  over  airfoils  at  Reynolds  numbers  based  on  the  airfoil  chord  length 
of  Rec  =500,000,  it  is  known  that  a  transitional  flow  region  exists  at  the  leading 

edge.  As  shown  by  Jang  et  al.  [Ref.  55],  the  incorporation  of  transition  modeling 
is  essential  for  the  successful  prediction  of  the  boundary  layers  on  airfoils. 

Empirical  formulations  are  used  in  engineering  methods  for  the  prediction 
of  flow  transition  in  boundary  layers.  Such  methods  assess  flow  transition  criteria 
by  (1)  using  momentum  thickness  as  a  parameter,  (2)  using  turbulence  models 
tuned  to  change  from  laminar  to  turbulent  flow  calculations  with  averaged 
transport  equations,  or  (3)  accounting  for  the  linear  stability  characteristics  of  the 
boundary  layer  by  calculating  amplitude  ratios  for  exponentially  growing 
eigenmodes.  These  formulations  have  been  used  with  success  for  incompressible 
and  subsonic  flows,  but  prediction  of  transition  location  with  these  methods  is 
uncertain  for  pressure-gradient-driven  compressible  flows. 

One  such  method  is  based  on  the  small  amplitude,  locally  parallel  flow 
assumptions,  which  enable  normal  mode  decomposition.  Computation  of  normal 
modes  and  their  growth  rates  give  a  criterion  to  predict  transition  onset  locations 
by  the  en  method.  Various  shortcomings  exist  with  this  method.  Optional 
strategies  for  calculating  (n)  factors  are  not  consistent  with  the  physical  constraints 
on  the  eigenmodes  under  the  assumptions  of  linear  stability  theory,  as  shown  by 
Mack  [Ref.  56]  and  Stuckert  et  al.  [Ref.  57].  The  method  is  typically  applied  for 
attached  flows  at  high  Reynolds  numbers. 
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The  method  used  in  this  work  was  previously  used  by  Walker  et  al.  [Ref.  30], 
showing  that  the  separation  bubble  occurring  on  a  NACA  65-213  airfoil  for 
incompressible  flow  and  Rec  =240,000  could  be  computed  successfully  by  using 

the  Chen-Thyson  [Ref.  20]  transition  length  model  in  combination  with  a  viscous - 
inviscid  interaction  method.  The  method  coupled  a  finite  difference  boundary- 
layer  code  with  a  potential  flow  code  for  the  outer  flow  field.  In  Reference  30,  to 
obtain  computed  results  that  were  in  agreement  with  the  measurements,  the 
transition  value,  Gy    had  to  be  chosen  in  a  range  between  20  and  40  rather  than 

the  recommended  value  of  1200.  Recent  experiments  by  Gostelow  et  al.  [Ref.  23] 
on  the  effects  of  free-stream  turbulence  and  adverse  pressure  gradient  on 
boundary-layer  transition  show  that  a  change  in  pressure  gradient  from  zero  to 
even  a  modest  adverse  level  is  accompanied  by  a  severe  reduction  in  transition 
length.  Larger  values  of  Gy    are  appropriate  for  attached  flows  at  high  Reynolds 

numbers.  Physically,  lower  transition  values  force  transition  to  take  place  over  a 
shorter  distance.  Recently,  Ekaterinaris  et  al.  [Ref.  58]  used  Michel's  empirical 
transition  onset  criterion  [Ref.  13]  and  Chen-Thy son's  transition  length  model  with 
some  success  to  compute  the  separation  bubble  on  the  NACA  0012  airfoil.  The 
objective  of  the  work  in  this  thesis  is  to  study  the  effect  of  systematic  variations  of 
the  transition  onset  location  and  transition  length  parameters  in  order  to  determine 
the  location  and  extent  of  the  transition  zone  that  best  models  the  measured 
separation  bubble  behavior.  Onset  of  transition  is  varied  from  the  initial  prediction 
given  by  Michel's  criterion,  which  is  based  on  high  Reynolds  number  data, 
Rec  >  106.  Michels'  criterion  predicts  that  boundary  layer  transition  has  occurred 
if  the  Reynolds  number  based  on  the  boundary  layer  momentum  thickness  (Ree)  is 
greater  than  the  quantity  { 1.174  (  1.0  +  22,400/Rex)Rex146},  where  Rex  is  based 
on  the  airfoil  surface  (  x )  location  from  the  stagnation  point. 

The  Chen-Thyson  intermittency  equation  (4.12)  was  incorporated  into  the 
N-S  code  as  was  done  for  the  boundary  layer  analysis.  The  equation  is  repeated 
here  for  completeness. 


ytr(x)  =  l-exp 


(-^)Re-(x-Xtr)i;^ 

'  tr  c  / 
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For  x  <  xtr,  ytr(x)  is  zero  and  downstream  from  transition  onset  increases 
exponentially  to  a  value  of  one,  which  corresponds  to  a  fully  turbulent  flow  region. 
The  GYtr  value  alters  the  slope  of  the  transition  curve  and  determines  the  length  of 
the  transition  region.   The  transition  value,  Gy  r ,  can  be  input  manually  and  was 

normally  varied  from  the  beginning  estimate  given  by  the  Cebeci  [Ref.  21] 
correlation  formula 

GYtr  =  213  [log(ReXtr )  -  4.7323]/3  (4.22) 

For  the  flow  transition  region,  an  effective  eddy  viscosity  is  obtained  by  scaling 
the  turbulent  eddy  viscosity  that  is  computed  from  the  mean  flow  by  ytr(x) . 
Therefore,  utrans  =  ytr(x)uturb  in  the  transitional  flow  region. 

The  computation  starts  with  a  laminar  solution  from  the  stagnation  point 
until  the  onset  of  transition  is  calculated  by  Michels'  criterion  or  transition  onset  is 
specified.  Chen-Thy son's  intermittency  equation  (4.12)  is  used  in  the  transition 
zone.  Aft  of  the  transitional  flow  region,  a  turbulence  model  is  used  for  flow 
analysis.  For  the  present  calculations,  the  two-layer  Baldwin-Lomax  [Ref.  59] 
turbulence  model  or  the  one  equation  Baldwin-Barth  [Ref.  60]  turbulence  model 
can  be  used  in  the  analysis.  The  Baldwin-Barth  model  is  used  for  all  the  results 
given  in  this  thesis.  The  effectiveness  of  various  turbulence  models  was 
investigated  by  Clarkson  et  al.  [Ref.  61]  and  Ekaterinaris  and  Menter  [Ref.  62]. 
The  conclusions  from  the  investigations  were  that  for  mildly  separated  flows,  such 
as  of  concern  for  this  work,  the  effect  of  the  turbulence  model  on  the  prediction  of 
aerodynamic  loads  was  not  very  significant. 

From  an  implementation  point  of  view,  the  transition  model  is  easily 
incorporated  into  numerical  analysis  codes,  is  simple  to  use,  and  can  be  combined 
with  any  turbulence  model.  The  model  accounts  for  local  flow  characteristics,  and 
other  effects,  such  as  curvature,  upstream  influence,  and  pressure  gradients  are 
included  only  through  the  computed  mean  flow.  The  model  is  also  independent  of 
free-stream  turbulence  level  and  other  flow  disturbances.  The  model  requires  an 
evaluation  of  boundary-layer  quantities  and  is  not  recommended  for  massively 
separated  flow  cases.    The  intermittency  function  has  streamwise  dependency, 
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while  the  effects  normal  to  the  wall  enter  through  the  turbulent  eddy  viscosity 
obtained  by  the  turbulence  model. 

5.         Boundary  Conditions  and  Gridding 

The  non-slip  condition  is  applied  for  velocities  on  the  airfoil  surface.  For 
the  unsteady  solutions,  the  surface  velocity  is  set  equal  to  that  obtained  by  the 
prescribed  airfoil  motion  as  follows: 


u  =  y(iUz-£tr|z), 


w  =  y(^Jlx-^x) 


Unsteady  solutions  for  the  oscillating  airfoil  are  obtained  by  rotating  the 
grid  with  the  airfoil.  For  subsonic  flow  conditions  at  the  inflow  and  outflow 
boundaries,  the  variables  are  evaluated  using  one-dimensional  Riemann  invariant 
extrapolation.  At  the  inflow  boundary,  there  are  three  incoming  characteristics 
and  one  outgoing  characteristic.  The  density  (p),  normal  velocity  (w),  and  the 
pressure  (p)  are  the  three  variables  specified  at  the  inflow  boundary.  The  fourth 
variable,  the  axial  velocity  (u),  is  extrapolated  from  the  interior  flow  region. 
Inflow  boundary  conditions  are  given  by 


/    2  \0/Y-D 


Pi  = 


Usi 


Sl  = 


(     \ 


=  flzfi(Rf-i;)f 


Ul=l/2(R?-  +  RJ), 


W]  =weo, 


Pi  = 


y  ) 


where  Rf  and  R2  are  the  incoming  and  outgoing  Riemann  invariants  given  by 


Rf  =  u00  +  2a00/(y-l), 


R2=u2-2a2/(y-l). 


At  the  outflow  boundary,  there  is  one  incoming  characteristic  and  three 
outgoing  characteristics.  Only  the  pressure  is  specified  at  the  outflow  boundary, 
while  the  other  variables  are  extrapolated  from  the  interior  flow  region.  For  the 
density  and  normal  velocity,  simple  first-order  extrapolation  is  used,  and  the  axial 
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outflow  velocity  is  obtained  from  the  zero-order  outgoing  Riemann  invariant. 
Outflow  boundary  conditions  are  given  by 


Pi=P2  u1=Rf-2a2/(y-l),  ai=VWPi. 


*i  =  VTP1/P1 '  w,=w2,  Pi=Poo- 

For  the  unsteady  analysis,  the  oscillating  airfoil  along  with  the 
computational  grid  rotates  about  the  one-quarter  chord  point  (x/c  =  0.25,  y/c  = 
0.0).  On  the  other  hand,  the  rectangular  LDV  measurement  grid  remains 
stationary  throughout  all  the  tests.  Therefore,  an  interpolation  scheme  was  used  to 
make  direct  comparisons  of  computed  and  measured  velocity  profiles  at  the 
stationary  nodes  of  the  measurement  grid  as  shown  in  Chapter  VII.  The 
computation  (275x81)  and  LDV  measurement  (6x21)  grids  near  the  leading  edge 
of  the  airfoil  were  shown  in  Figure  3.4. 

D.        COMPUTATIONAL  UNCERTAINTIES 

The  computations  were  done  with  grid  densities  ranging  from  275x81 
(baseline  grid)  to  601  x  161  with  refined  points  in  the  transition  region  to  establish 
grid  independence.  The  baseline  grid  had  16  surface  points  in  the  vicinity  of  the 
separation  bubble,  while  the  dense  grid  had  86  points  in  the  same  region.  The 
boundary  layer  had  17  and  38  points  in  the  respective  grids.  The  oscillatory  airfoil 
solutions  were  started  with  a  converged  steady  solution  at  8  degrees  angle  of  attack. 
The  density  residual  for  the  steady  angle-of- attack  solutions  was  on  the  order  of  10"6 
after  4000  time  steps,  with  each  time  step  being  0.005.  The  unsteady  analysis  used 
the  same  time-step  size  with  density  residuals  on  the  order  of  10-4,  which  was  not 
surprising  because  of  the  continual  flow  changes  that  were  imposed  by  the  oscillating 
airfoil.  Only  minor  differences  were  noted  in  the  computed  solutions.  For  example, 
at  8  degrees  angle  of  attack,  the  maximum  differences  in  pressure  and  skin-friction 
coefficient  were  less  than  1.5%,  occurring  around  the  suction-pressure  peak  (-3.4 
versus  -3.45)  near  the  airfoil  leading  edge.  For  the  oscillatory  motion  cases,  the 
computed  solution  continued  for  three  cycles.  The  second  and  third  cycle  values  for 
lift,  drag,  and  moment  coefficients  were  identical,  and  the  results  are  shown  for  the 
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third  cycle.    Therefore,  the  275x81  point  grid  was  used  for  both  the  steady  and 
unsteady  solutions  that  are  presented  in  this  thesis. 
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V.  EXPERIMENTAL  RESULTS 

This  chapter  presents  the  results  from  the  LDV  and  PDI  experimental 
methods  described  in  Chapter  III.  Measurements  were  made  to  characterize 
steady,  high-  and  low-amplitude  oscillatory,  deep  and  lightly  stalled  compressible 
flows  with  separation  bubbles  near  the  airfoil  leading  edge  that  affect  the 
inception  and  development  of  the  dynamic-stall  process.  Both  the  LDV  and  PDI 
measurements  gave  valuable  data  for  evaluation  of  computational  capabilities  to 
predict  the  leading-edge  airfoil  flow  field. 

Initially,  the  leading-edge  flow  field  over  the  NACA  0012  airfoil  is  best 
described  by  the  interferograms  from  the  PDI  measurements.  Each  interference 
fringe  gives  a  line  of  constant  density  with  an  associated  Mach  number.  The 
fringe  map  can  also  be  used  to  quantify  the  local  pressure  distributions  on  the 
airfoil,  as  shown  in  the  next  section.  The  registration  markers  for  identifying 
proper  airfoil  orientation  were  described  in  Chapter  III  and  are  seen  in  all  the  PDI 
Figures.  The  interferograms  focus  on  a  region  from  slightly  forward  of  the  airfoil 
leading  edge  to  30-45%  chord  downstream. 

A.        FLOW  OVER  STEADY  NACA  0012  AIRFOIL 

Interferograms  are  shown  for  a  =  0°  in  Figure  5.1.  The  stagnation  point  is 
at  the  leading  edge,  which  is  at  the  center  of  the  smallest  closed  density  contour. 
The  velocity  increases  from  the  stagnation  point  in  both  directions  over  the  upper 
and  lower  surfaces.  The  stagnation  point  in  the  flow,  which  moves  under  the 
airfoil  leading  edge  with  increasing  angle  of  attack,  quantifies  the  density  field 
and  related  Mach  contours  for  the  complete  interferogram.  Good  flow  symmetry 
is  seen  for  the  upper  and  lower  surfaces  of  the  airfoil,  as  expected  at  a  =  0° . 

At  a  =  6° ,  the  stagnation  point  has  moved  under  the  leading  edge  of  the 
airfoil,  and  the  number  of  fringes  has  increased  from  eight  to  16.  The  additional 
fringes  are  concentrated  over  the  upper  airfoil  leading-edge  surface.  All  density 
contours  present  a  continual  and  gradual  curvature  with  the  upper  surface,  that  is, 
no  presence  of  a  separation  bubble  is  indicated  in  Figure  5.2. 
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Figure  5.1.   Interferogram  of  the  Flow  Near  the  Airfoil  Leading  Edge.    M^  =  0.3, 
a  =  0°,  k  =  0.0,  Re  =  540,000. 


Figure  5.2.   Interferogram  of  the  Flow  Near  the  Airfoil  Leading  Edge.    M^  =  0.3 


a  =  6°,  k  =  0.0,  Re  =  540, 000. 
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In  Figure  5.3  at  a  =  8°,  three  additional  fringes  are  seen  over  the  airfoil 
leading  edge,  where  the  constant  density  contours  show  a  change  in  shape  from 
those  previously  seen  at  a  =  6°.  Just  downstream  of  and  not  belonging  to  the 
suction  peak  itself,  the  contours  close  to  the  airfoil  surface  turn  almost  parallel 
with  the  airfoil  surface  for  some  distance.  The  region  encompassing  five  fringes 
shows  very  little  gradient  in  the  streamwise  direction,  but  at  the  downstream  ends 
turn  abruptly  towards  the  airfoil  surface.  The  fringes  meet  the  airfoil  surface 
almost  perpendicularly  in  contrast  to  the  gradual  low  angles  that  are  shown  at 
a  =  6°.  The  shape  of  the  fringe  contours  suggests  the  presence  of  a  laminar 
separation  bubble  because,  subsequently,  the  flat  contours  are  shown  to 
correspond  to  a  plateau  in  the  pressure  distribution.  The  test  conditions  for  the 
NACA  0012  airfoil  at  Re  =  540,000  are  appropriate  for  laminar  separation 
bubbles.  The  formation  of  laminar  separation  bubbles  was  already  confirmed  by 
LDV  measurements  over  the  leading  edge  of  the  airfoil  for  low-frequency  forced 
oscillations  [Ref.  41  and  34].  The  contour  lines  outside  the  separation  bubble 
turn  gradually  downstream  and  are  seen  to  define  the  boundary  layer  over  the 
airfoil  surface  downstream  of  the  separation  bubble. 


Figure  5.3.   Interferogram  of  the  Flow  Near  the  Airfoil  Leading  Edge.    M^  =  0.3, 
a  =  8°,  k  =  0.0,  Re  =  540,000. 
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As  the  flow  accelerates  over  the  leading  edge  of  the  airfoil  at  a  =  10°,  five 
additional  fringes  and  higher  density  gradients  develop  as  shown  in  Figure  5.4. 
A  dark  region  appears  at  the  leading  edge  of  the  airfoil  because  much  of  the  light 
has  been  severely  deflected  by  the  steep  gradients  converging  in  the  area,  but  it 
has  been  shown  not  to  affect  the  fine  details  of  the  pressure  distribution.  Just 
downstream  of  the  suction  peak,  where  the  airfoil  surface  curvature  is  also 
reduced,  the  contours  again  are  seen  to  "plateau"  for  some  distance  along  the 
airfoil  surface  just  as  seen  in  Figure  5.3.  Also,  the  fringes  turn  nearly 
perpendicular  to  the  airfoil  at  the  aft  end  of  the  separation  bubble  and  a 
thickening  boundary  layer  is  seen  downstream  of  the  separation  bubble. 


Figure  5.4.   Interferogram  of  the  Flow  Near  the  Airfoil  Leading  Edge.    M^  =  0.3. 
a  =  10°,  k  =  0.0,  Re -540,000. 


Figures  5.5  and  5.6  show  the  unsteady  nature  of  the  flow  field  in  the 
vicinity  of  the  airfoil  stalling  angle  of  attack  at  a  =  12°.  A  separated  flow 
condition  is  shown  in  the  first  case  at  (a)  a  =  1 1.95°.  and  in  the  second  case,  a  well 
defined   separation  bubble  resides  near  the   leading  edge   at    (b)  a  =  12.05°. 
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However,  the  fringes  turning  vertical  towards  the  airfoil  surface  in  the  boundary 
layer  at  the  downstream  end  of  the  separation  bubble  are  the  first  evidence  of 
stall  [Ref.  63]. 


Figure  5.5.   Interferogram  of  the  Flow  Near  the  Airfoil  Leading  Edge.    MM  =0.3, 
a  =  1 1.95°,  k  =  0.0,  Re  =  540,000. 


Figure  5.6.   Interferogram  of  the  Flow  Near  the  Airfoil  Leading  Edge.    M^  =  0.3, 
a  =  12.05°,  k  =  0.0,  Re  =  540, 000. 


59 


1.         Pressure  Distributions  From  the  PDI  Interferograms 

The  pressure  plots  were  obtained  by  mapping  the  x-location  of  the 
intersection  of  each  fringe  with  the  airfoil  surface.  The  fringe  numbers  were  then 
converted  to  Cps  for  the  pressure  distributions.  The  resultant  pressure 
distributions  for  0  <  x/c  <  0.2  derived  from  the  interferograms  are  given  in  the 
next  series  of  plots. 

In  Figure  5.7,  the  suction  peak  (Cp)  is  shown  to  increase  by  at  least  (-1)  for 
every  2°  increase  in  angle  of  attack  from  a  =  6°  to  10°.  No  pressure  plateau  is 
seen  downstream  of  the  suction  peak  for  a  =  6° .  The  innermost  fringe  that  could 
be  defined  at  a  =  6°  extended  for  a  short  distance  along  the  airfoil  upper  surface, 
resulting  in  a  flat  suction  pressure  peak.  The  flat  pressure  plateau  aft  of  the 
suction  peak  for  a  =  8°  and  a  =  10°  is  indicative  of  a  separation  bubble, 
corresponding  to  the  density  contours  that  turned  parallel  to  the  airfoil  surface 
before  turning  abruptly  towards  the  airfoil  surface  at  the  downstream  edges.  Aft 
of  the  pressure  plateau,  a  second  region  of  adverse  pressure  gradient  is  shown  for 
0.05  <  x/c  <  0.075.  For  example,  here  the  pressure  gradient,  dCp/d(x/c),  is  125 
while  the  Cp  increases  from  -3.6  to  -2.15  for  a  =  10°.  The  adverse  pressure 
gradient  aft  of  the  suction  pressure  peak  is  generally  less  (  =  60-70).  The  pressures 
then  gradually  increase  over  the  rest  of  the  airfoil.  The  chordwise  extent  of  the 
pressure  plateau  region  may  be  used  to  estimate  the  size  of  the  separation  bubble. 
At  a  =  8°  the  bubble  begins  at  3.0%  chord  and  at  a  =  10°,  the  separation  bubble 
location  begins  at  2.0%  chord,  indicating  onset  of  the  bubble  closer  to  the 
leading  edge  at  higher  angles  of  attack  as  a  result  of  an  increased  adverse 
pressure  gradient.  The  downstream  edge  of  the  separation  bubble  ranges  from 
5.0%  chord  at  a  =  10°  to  6.0%  chord  at  a  =  8°.  Presumably,  the  separation 
bubble  extends  to  downstream  locations  that  are  part  way  down  the  adverse 
pressure  gradient  or  where  the  turbulent  flow  pressure  profile  intersects  the 
pressure  profile  with  separation  bubbles  for  this  flow,  much  the  same  as  that 
shown  by  Tani  [Ref.  5]  for  a  higher  Reynolds  number  flow. 

A  pressure  plateau  is  present  downstream  of  the  suction  peak  for  the  flow 
at  each  angle  of  attack,  a  =  7°,  9°,  and  1 1°,  presented  in  Figure  5.8.  The  suction 
pressure  peaks  are  shown  to  move  forward  slightly  on  the  airfoil  surface  with 
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increasing  angle  of  attack.  The  length  of  the  pressure  plateau  decreases 
noticeably  with  increasing  angle  of  attack.  The  adverse  pressure  gradient 
downstream  of  the  plateau  and  separation  bubble  is  125  for  a  =  9°  and  decreases 
to  75  for  a  =  11°  in  the  range  0.04  <  x/c  <  0.065.  The  suction  pressure  peak  at 
a  =  11°  is  -4.45  compared  to  -4.25  at  a  =  10°,  indicating  a  flow  that  is 
approaching  its  full  lifting  capacity. 
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Figure  5.7.    Pressure  Distributions  Derived  From  PDI  Measurements.  MM  =0.3, 
a  =  6°,  8°,  and  10°,  k  =  0.0,  Re  =  540,000. 
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Figure  5.8.    Pressure  Distributions  Derived  From  PDI  Measurements.  M^  =0.3, 
a  =  7°,  9°,  and  11°,  k  =  0.0,  Re  =  540,000. 
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Pressure  distributions  for  a  stalled  flow  condition  at  11.95°  angle  of  attack 
and  a  flow  on  the  verge  of  stall  at  12.05°  angle  of  attack  are  shown  in  Figure  5.9. 
The  pressures  correspond  to  the  interferograms  in  Figure  5.5  and  5.6.  Indicative 
of  a  flow  on  the  verge  of  stall,  the  pressure  peak  at  12.05°  angle  of  attack  has 
decreased  slightly  from  that  previously  shown  at  1 1  °  angle  of  attack.  The  plateau 
region  shows  the  beginnings  of  breakdown  as  its  downstream  end  is  lifting  away 
from  the  surface.  The  flat  (stalled)  pressure  profile  at  1 1.95  °  angle  of  attack  gives 
pressure  levels  below  those  previously  shown  for  6  °  angle  of  attack  forward  of 
5%  chord  over  the  airfoil  upper  surface. 
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Figure  5.9.    Pressure  Distributions  Derived  From  PDI  Measurements.  M^  =0.3, 
cc  =  11.95°  and  12.05°,  k  =  0.0,  Re  =  540,000. 

2.         Mach  Number  From  PDI  Density  Ratio 

Figure  5.10  gives  a  representative  case  of  fringes  that  have  been  converted 
to  digital  form  for  a  =  10°  and  a  more  detailed  investigation,  indicating  the  density 
ratio  value  on  selected  contours.  The  0.944  density  contour  (one  below  0.935) 
that  was  on  the  upper  surface  of  airfoil  at  a  =  0°  now  meets  the  airfoil  leading 
edge  from  nearly  the  horizontal.  The  flow  is  highly  accelerated  over  the  leading 
edge  to  2.34  times  the  free-stream  Mach  number  (M  =  0.701).    By  including  all 
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fringes  downstream  of  the  suction  peak  that  run  parallel  with  the  airfoil  surface 
and  then  turn  perpendicular  to  the  airfoil  surface  downstream  of  5%  chord,  the 
separation  bubble  is  estimated  to  extend  from  0.025  <  x/c  <  0.065  along  the 
surface  of  the  airfoil.  Using  the  same  fringes,  the  separation  bubble  is  estimated  to 


be  2%  chord  in  height. 
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Figure  5.10.  Density  Contours  From  Interferograms  in  Flow  Near  Airfoil  Leading 
M^=0.3,  a  =  10°,  k  =  0.0,  Re  =  540,000. 


Edge 


3.         U-Component  (LDV)  Velocity  Contours 

The  U-component  velocity  contours  in  the  proximity  of  the  airfoil  surface 
are  given  for  8°  angle  of  attack  in  Figure  5.11.  A  high- velocity  region  extends 
from  5  to  9%  chord  at  heights  of  0.021  >  y/c  >  0.017  above  the  airfoil  surface 
and  is  slightly  aft  of  the  high-velocity  region  indicated  in  the  PDI  measurements. 
The  high-velocity  region  is  forward  of  the  surface  locations  0.10  <  x/c  <  0.117, 
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where  a  minimum  distance  of  1.6%  chord  exists  between  the  measurements  and 
the  airfoil  surface.  A  minimum  measured  velocity  in  the  flowfield  of  1.07  U^ 
occurs  above  the  airfoil  leading  edge. 
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Figure  5.11.  Normalized  Velocity  Contours  (U-component,  LDV)  Near  the  Airfoil 
Leading  Edge.    M^O.3,  a  =  8°,  k  =  0.0,   Re  =  540,000. 
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Velocity  contours  for  a  =  12.05°  are  shown  in  Figure  5.12.  The 
measurement  region  relative  to  the  airfoil  position  is  similar  to  that  at  8  °  angle  of 
attack.  Slightly  higher  velocities  (1.6811^)  are  noted  because  of  the  increased 
angle  of  attack,  which  occur  in  a  slightly  larger  region  from  4.5  to  9%  chord  at 
heights  from  0.025  >  y/c  >  0.021  above  the  airfoil  surface.  Measurement  points 
closer  to  the  airfoil  surface  could  not  be  optically  accessed  because  of  laser  beam 
blockage  by  the  airfoil  surface.  The  highest  velocity  region  encompasses  the 
separation  bubble  location  near  the  airfoil  surface. 
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Figure  5.12.  Normalized  Velocity  Contours  (U-component,  LDV)  Near  the  Airfoil 
Leading  Edge.    M^O.3,  a  =  12.05°,  k  =  0.0,  Re  =  540,000. 
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B.        UNSTEADY/OSCILLATING  NACA  0012  AIRFOIL 

LDV  results  are  shown  in  this  section  for  two  separate  test  cases  in  which 
the  airfoil  oscillation  amplitude  is  varied  about  the  mean  angle  of  attack,  a  =  10°, 
such  that  a  =  10° -10°  sin  cot  and  a  =  10°  -2°  sin  cot  in  a  Moo=0.3  flow.  The 
reduced  frequency  ( k  =  2  jcf c/UM )  based  on  airfoil  chord  for  each  test  case  was 
0.1.  As  for  the  steady  flow  results,  the  test  Reynolds  number  is  540,000  and  the 
airfoil  is  the  NACA  0012. 

Several  PDI  images,  focusing  on  the  first  30  to  40%  airfoil  chord,  are  used 
initially  to  describe  the  flow  physics  for  the  high-amplitude  forced-oscillation  test. 

1.         Deep  Dynamic  Stall  ( a  =  10°  - 10°  sin  co  t ) 

This  test  case  (Moo=0.3,  cc  =  10°  -10°  sin  cot)  resulted  in  massive  flow 
separation  at  the  higher  angles  of  attack.  LDV  measurements  were  more  difficult 
in  the  highly  accelerated  and  separated  flow  that  showed  extremely  complex 
flow  physics. 

The  PDI  images  in  Figure  5.13  show  the  flow  physics  for  the  high- 
amplitude  oscillation  case  during  airfoil  upward  motion.  At  0  degree  angle  of 
attack,  the  flow  has  fully  established  to  the  point  where  the  picture  is  almost 
identical  to  that  previously  shown  for  the  steady  case  at  0  degree  angle  of  attack. 
For  the  unsteady  case,  stall  does  not  occur  at  12°  angle  of  attack  on  airfoil 
upward  motion.  As  is  well  known  from  Carr  and  Chandrasekhara  [Ref.  64],  the 
process  of  dynamic  stall  in  this  case  starts  with  the  formation  of  the  dynamic  stall 
vortex  at  angles  of  attack  slightly  higher  than  static  stall.  Until  the  vortex  is  shed 
into  the  wake,  the  airfoil  continues  to  generate  lift.  The  vortex  release  (shedding) 
angle  as  measured  earlier  by  schlieren  flow  visualization  experiments  was  shown 
to  be  a  =  15.9°  [Ref.  26].  The  PDI  image  in  Figure  5.13d  indicates  that  the  flow 
has  almost  reached  the  fully  stalled  condition. 
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(a)ct  =  0°. 


(b)  a  =  6C 


(c)a  =  12°. 


(d)a  =  15.5c 


Figure  5.13.  Interferograms  of  the  Flow  Near  the  Leading  Edge  of  the  Oscillating 
Airfoil.   ML  =0.3,   a  =  10° -10°  sin  cot,  k  =  0.1. 
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A  progression  of  flow  reattachment  is  shown  in  Figure  5.14  during  the 
airfoil  downward  motion.  The  flow  remains  highly  separated  at  12°  angle  of 
attack  on  airfoil  downward  motion,  even  though  it  is  well  below  the  dynamic  stall 
angle  of  attack  of  15.9°.  Flow  reattachment  has  begun  at  a  =  11°  with  several 
fringe  now  contacting  the  airfoil  surfa  :  over  the  leading  edge.  As  documented 
in  [Ret.  42]  and  as  shown  in  Figure  5.14d,  complete  flow  reattachment  does  not 
happen  until  the  airfoil  passes  9°  angle  of  attack.  Thus,  a  larger  hysteresis  loop  is 
present  for  the  flow. 


(a)  a  =  12°. 


(b)a  =  ll°. 


(c)  a  =  10.5°. 


(d)a  =  9°. 


Figure  5.14.  Interferograms  of  the  Flow  Near  the  Leading  Edge  of  the  Oscillating 
Airfoil.  MM=0.3,  a  =  10°  -10°  sin  cot,  k=  0.1. 
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U-component  velocities  throughout  the  oscillation  cycle  (vs.  phase  angle) 
are  shown  in  Figure  5.15  at  x/c  =  -0.100  offset  by  0.211^  at  each  increasing 
height  for  presentation  purposes.  Even  at  this  upstream  location,  the  effect  of  the 
oscillating  airfoil  on  the  flow  is  felt  at  all  heights  above  the  airfoil  surface.  At  the 
lowest  measurement  height  (y/c  =  0.083),  U  velocities  are  slightly  less  than  the 
free-stream  speed  and  data  are  not  continuous  because  the  LDV  beam  is  blocked 
by  the  airfoil.  A  sinusoidal  relationship  exists  between  the  phase  angles  (angles 
of  attack)  along  the  x-axis,  and  the  U  velocity  tends  to  follow  the  airfoil 
oscillatory  motion  up  to  the  dynamic  stall  angle  of  attack,  a  =  15.9°  or  <j>  =  216°. 
U  velocities  peak  at  10%  above  the  free-steam  value  just  before  reaching  the 
dynamic  stall  angle  of  attack  and  then  fall  up  to  20%  because  of  the  separated 
flow  as  the  angle  of  attack  increases  to  20°.  The  reduced  velocities  continue 
until  the  flow  begins  to  reattach  at  a  =  1 1°  on  the  airfoil  downstroke  to  minimum 
values  in  the  cycle  just  after  passing  the  dynamic  stall  angle  of  attack. 
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Figure  5.15.    Phase  Distribution  of  U-Component  Velocities  at  x/c  =  -0.100. 
(Velocity  offset  by  0.2  UTO  at  each  height.)  Moo=0.3,  a  =  10°  -10°  sin  cot,  k=  0.1. 
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V-component  (vertical)  velocities  forward  of  the  airfoil  surface  vary  mo- ■ 
than  the  U-component  during  the  airfoil  oscillations  (5%  to  40%  of  free-streain 
values)  before  dynamic  stall  occurs,  as  shown  in  Figure  5.16.  Again,  the  velocities 
follow  the  motion  of  the  airfoil  to  the  dynamic  stall  angle  of  attack.  Peak 
velocities  during  the  oscillation  cycle  occur  slightly  after  those  seen  for  the  U- 
component.  Velocity  reductions  after  dynamic  stall  are  not  as  extensive  as  for  the 
U-component  because  the  streamlines  are  deflected  upward  as  the  flow 
undergoes  deep  stall.  The  V  velocity  slowly  decreases  until  flow  reattachment  is 
complete  [Ref.  42].  Disregarding  the  velocity  offset,  V  velocities  overlay  each 
other  at  the  bottom  two  heights  with  average  values  between  0.05  and  0.4  of  the 
free-stream  speed.  At  increasing  heights,  V  velocities  decrease  slightly. 


8 


1.6 


10 


O  Y/C=  0.167 
A  Y/C=  0.150 
X      Y/C=  0.133 


O  Y/C=  0.117 
a  Y/C=  0.100 
♦      Y/C=  0.083 


90.0  180.0  270.0 

Phase  Angle,  <p 

0  10  20 

Angle  of  Attack,  a 


360.0 


10 


Figure  5.16.    Phase  Distribution  of  V-Component  Velocities  at  x/c  =  -0.100. 
(Velocity  offset  by  0.2  LL  at  each  height.)  Mx  =  0.3,  a  =  10°  -10°  sin  cot,  k  =  0.1. 
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More  dramatic  changes  in  the  flow  are  seen  above  the  airfoil  surface  at  x/c 
=  0.067  shown  in  Figure  5.17.  Validated  average  U-component  velocities  vary 
widely  from  0.45  to  1.5  U^  throughout  the  oscillation  cycle.  The  missing  data  at 
the  lower  heights  occur  because  of  laser-beam  blockage  by  the  airfoil,  and  at  y/c 
=  0.017  a  lack  of  validated  particles  is  the  cause  of  missing  data  in  some  cases. 
The  velocities  at  this  location  over  the  airfoil  surface  vary  more  widely  than  those 
forward  of  the  airfoil  before  reaching  the  dynamic  stall  angle  of  attack.  At  y/c  = 
0.167  and  a  =  8°,  velocities  are  at  the  free-stream  value,  but  those  at  a  =  15.9°  are 
50%  greater  than  the  free-stream  velocity.  At  the  three  uppermost  (y/c)  locations 
(above  or  out  of  the  separation  bubble)  after  passing  the  dynamic  stall  angle  of 
attack,  the  velocity  drops  up  to  50%  with  major  velocity  fluctuations  (0.5  of  the 
free-stream  velocity)  noted  in  the  measurements.  At  the  lower  heights  (y/c  = 
0.083  and  0.100),  U  velocities  fall  off  rapidly  already  at  7  and  8°  angle  of  attack 
from  1.35  and  1.4-U^,  getting  as  low  as  0.411^  before  the  angle  of  attack  reaches 
13°.  This  dramatic  drop  in  velocity  is  attributable  to  velocity  measurements 
within  the  outer  region  of  the  separation  bubble. 
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Figure  5.17.    Phase  Distribution  of  U-Component  Velocities  at  x/c  =  0.067. 
(Velocity  offset  by  0.2  UTC  at  each  height.)  M00  =  0.3,  a  =  10° -10°  sin  cot,  k=  0.1. 
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In  Figure  5.18,  V-component  velocities  at  x/c  =  0.067  near  the  dynamic 
stall  angle  of  attack  show  rapid  increases,  starting  just  before  and  then  peak  at 
dynamic  stall.  For  example,  at  y/c  =  0.133,  the  V  velocities  rise  abruptly  from  0.1 
to  0.55  U^  as  angle  of  attack  increases  from  13  to  16°.  During  this  phase  of 
oscillation,  the  dynamic  stall  vortex  grows,  pushing  the  outer  edge  of  the 
enveloping  shear  layer  outwards.  At  a  =  15.9°  ,  the  dynamic  stall  vortex  is  shed 
[Ref.  26]  and  the  shear  layer  detaches  from  the  surface  everywhere,  except  at  the 
airfoil  leading  edge.  The  large  V-component  velocities  and  the  simultaneous  drop 
in  the  U-component  velocities  are  further  evidence  of  separation  bubble  bursting, 
which  initiates  the  dynamic  stall  process. 
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Figure  5.18.    Phase  Distribution  of  V-Component  Velocities  at  x/c  =  0.067. 
(Velocity  offset  by  0.2 UM  at  each  height.)  M^  =  0.3,  a  =  10°  -10° since  t,  k=  0.1. 
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Velocity  magnitude  profiles  are  shown  at  x/c  =  0.067  and  angles  of  attack 
from  7  to  16°  in  Figure  5.19.  Velocities  at  y/c  =  0.083  continue  to  decrease  from 
0.83  to  0.65  U*.  as  the  angle  of  attack  increases  from  7  to  13°.  The  velocity 
profiles  in  the  laminar  shear  layer  initially  have  full  laminar  shapes,  and  then 
exhibit  slight  inflections  at  the  higher  angles  of  attack.  From  14  to  16°  angle  of 
attack,  inverted  "S"  profiles  appear,  indicating  unusual  shear  patterns  in  the  flow. 
At  y/c  =  0.10  and  for  a  >  14°,  velocities  are  as  low  as  0.5  U^. 
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Figure  5.19.   Velocity  Magnitude  IVI  Profiles  at  x/c  =  0.067.  (Profiles  offset  by 
0.21^  ateachAOA.)   1^=0.3,  a  =  10° -10°  sin  cot,  k=  0.1. 

The  previous  figures  describe  the  dramatic  flow  physics  that  are  part  of  the 
high-amplitude  deep  dynamic  stall  flow.  Computationally,  the  problem  becomes 
intractable  for  such  flows,  and  thus,  the  low  amplitude  case  was  measured,  which 
is  shown  to  present  the  case  of  light  dynamic  stall. 
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2.         Light  Dynamic  Stall  Flow  ( a  =  10°  -  2°  sin  co  t ) 

Light  dynamic  stall  occurs  when  the  static  stall  angle  is  exceeded  only 
slightly.  The  extent,  severity,  and  duration  of  the  separation  is  significant,  yet 
much  less  than  for  deep  stall.  A  distinguishing  feature  of  light  dynamic  stall  is 
that  the  vertical  extent  of  the  viscous  zone  remains  thin  (on  the  order  of  airfoil 
thickness)  and  the  thin-layer  N-S  equations  can  be  used  for  analysis  [Ref.  24]. 
Detailed  experimental  flow  results  are  presented  for  the  low-amplitude  forced- 
oscillation  test  case  at  M^  =0.3,  where  the  top  of  the  oscillation  cycle  is  near  the 
airfoil  static  stall  incidence  angle.  The  new  finding  in  the  measurements  show 
that  light  dynamic  stall  occurs  on  the  airfoil  downstroke  below  the  static  stall 
angle  of  attack,  providing  an  opportunity  to  evaluate  the  capability  of  the  thin- 
layer  N-S  equations  to  describe  the  flow  characteristics,  as  done  in  Chapter  VII. 

a.         PDI  Images 

The  images  in  Figure  5.20  show  the  flow  physics  throughout  the 
oscillation  cycle  for  this  test  case.  At  the  bottom  (a  =  8°)  of  the  oscillation  cycle, 
the  flow  has  established  to  the  point  that  the  picture  closely  resembles  that 
previously  shown  for  the  steady  case  at  8  °  angle  of  attack.  This  is  not  surprising 
because,  in  fact,  the  pitch  rate  is  zero  deg/s  for  an  instant  during  reversal  of  airfoil 
motion  both  at  the  bottom  and  top  (a  =  12°)  of  the  oscillation  cycle,  and  the 
reduced  frequency  is  low.  As  described  previously,  the  suction  peak  occurs  on 
the  upper  airfoil  surface  at  the  innermost  fringe.  Immediately  downstream  of  the 
suction  peak,  the  fringes  run  parallel  with  the  airfoil  surface  for  some  length  and 
then  turn  abruptly  towards  the  airfoil  surface.  As  before,  the  distance  between 
the  suction  peak  and  the  abrupt  turning  end  of  the  fringes  indicates  the  extent  of 
the  separation  bubble  along  the  airfoil  surface,  which  is  present  at  all  angles  of 
attack  on  the  airfoil  upstroke.  Later,  pressure  distributions  derived  from  the  PDI 
images  clearly  show  the  separation  bubble  at  a  =  8°,  11°,  and  11.9°.  The  density 
contours  in  the  outer  flow  where  no  separation  bubble  is  present  merge  gradually 
with  the  airfoil  upper  surface  and  turn  downstream,  defining  the  edge  of  the 
boundary  layer  downstream.  As  the  airfoil  angle  of  attack  increases,  additional 
leading-edge  fringes  develop.  However,  the  flow  everywhere  indicates  the 
bubble  is  fully  attached  to  the  top  of  the  cycle. 
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Surprisingly,  shortly  after  the  airfoil  starts  the  downward  motion  from  the 
top  of  the  oscillation  cycle,  a  stalled  condition  is  indicated  by  the  fringes  in  Figure 
5.20d.  The  flow  reattachment  process  occurs  quickly  and  is  well  underway  at 
a  =  10.5°  during  airfoil  downward  motion.  The  PDI  images  shown  in  Figure  5.20f 
at  a  =  9°  during  continued  downward  airfoil  motion  show  that  another 
separation  bubble  is  fully  formed  at  the  airfoil  leading  edge.  The  light  dynamic 
stall  condition  occurs  even  though  the  peak  angle  of  attack  in  the  oscillation 
cycle  is  at  or  slightly  below  the  static  stall  angle  of  attack.  In  addition,  a  mild  stall 
was  verified  in  the  flow  field  vorticity  levels  by  the  LDV  measurements  [Ref.  41]. 


(a)  a  =  8°  (bottom). 


(b)  a  =  11°  (up). 


Figure  5.20.     Interferograms  of  the  Flow  Near  the  Airfoil  Leading  Edge. 
MTC=0.3,  a  =  10° -2°  sin  cot,  k  =  0.1. 
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(c)  a  =  11. 9° (up). 


(d)  a  =  11.5° (down). 


(e)a  =  10.5°(down). 


(f)  a  =  9° (down). 


Figure  5.20.  (Contd.) 

b.         Pressure  Distribution  Derived  From  PDI  Images 

Pressure  distributions  derived  from  the  PDI  images  near  the  airfoil  leading 
edge  are  shown  in  the  next  two  figures  for  this  low-amplitude  oscillatory  case. 
Surface  pressure  coefficients  for  airfoil  upward  motion  are  shown  in  Figure  5.21. 
The  suction  pressure  peak  maxima  occur  in  the  first  2%  chord,  moving  forward 
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with  increasing  angle  of  attack,  as  expected.  The  suction  peak  is  forward  of  1% 
chord  at  a  =  1 1°.  The  pressure  peaks  range  from  Cp  =  -3.0  at  8°  angle  of  attack 
to  Cp  =  -4.3  at  1 1  °  angle  of  attack.  However,  as  the  top  of  the  oscillation  cycle  is 
approached  (a  =  12°),  the  pitch  rate  goes  to  zero  and  the  minimum  Cp  decreases 
to  -3.8,  as  shown  in  Figure  5.22.  The  pressure  plateaus  aft  of  the  pressure  peaks, 
that  have  been  alluded  to  numerous  times,  indicate  that  the  separation  bubble 
moves  forward  with  increasing  angle  of  attack  from  8  to  11°,  but  the  extent  along 
the  airfoil  surface  remains  nearly  constant.  Aft  of  the  pressure  plateau,  the 
adverse  pressure  gradients  are  60  to  70  (compared  to  110  to  125  for  the  steady 
results)  while  the  Cp  decreases  from  -3.5  to  -1.5  over  a  small  chordwise  distance, 
0.045  <  x/c  <  0.07. 
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Figure  5.21.    Pressure  Distributions  Near  the  Airfoil  Leading  Edge  (upstroke). 
M    =0.3,  a  =  10° -2°  sin  cot,  k  =  0.1. 


Surface  pressure  distributions  during  the  airfoil  downward  motion  are 
given  in  Figure  5.22.  The  vortex  shedding  and  light  dynamic  stall  process  shown 
in  Figure  5.20d  is  also  indicated  by  the  flat  pressure  distribution  and  a  minimum 
Cp  of  -0.9  at  a=  11.5°.  The  PDI  images  showed  that  flow  reattachment  had 
already  begun  at  a-  10.5°  angle  of  attack.    This  is  also  shown  by  the  suction 
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pressure  peak  for  a  -  10°  that  has  increased  appreciably  (Cp  =  -2.0),  but  is  well 
below  the  Cp  =  -3.7  during  the  airfoil  upstroke.  At  a  =  9°,  the  separation  bubble 
plateau  is  again  evident  in  the  pressure  distribution.  However,  the  peak  suction 
pressure  is  lower  than  at  the  same  angle  of  attack  during  airfoil  upward  motion.  If 
the  stall  had  not  occurred  during  the  oscillation  cycle,  pressure  lag  effects  would 
give  higher  suction  peaks  on  airfoil  downward  motion,  much  like  those  calculated 
with  both  the  potential  (UPOT)  and  N-S  codes,  as  shown  in  Chapter  VI. 
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Figure  5.22.  Pressure  Distributions  Near  the  Airfoil  Leading  Edge  (downstroke). 
M^O.3,  a  =  10° -2°  sin  cot,  k  =  0.1. 

c.         Phase  Distribution  of  LDV  Velocities  at  (X/C)  Locations 

The  next  series  of  figures  present  the  results  from  an  extensive  array  of 
LDV  measurements.  To  review,  the  LDV  probe  volume  diameter  for  the  tests  is 
0.191  millimeter  (0.25%  chord  in  height),  whereas  the  height  of  either  a  laminar  or 
turbulent  boundary  layer  without  a  separation  bubble  is  0.121  millimeter  and 
0.184  millimeter  at  5%  chord  downstream,  respectively.  Thus,  the  LDV 
measurement  of  reverse  flows  at  the  airfoil  leading  edge  is  nearly  impossible 
because  of  the  relative  size  of  the  probe  volume  diameter  and  the  thin  reverse- 
flow  region  in  the  separation  bubble.    In  addition,  the  particles  tend  to  deflect 
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away  from  the  airfoil  surface  along  with  the  separating  shear  layer  rather  than 
follow  the  flow  into  the  separation  bubble.  Particle  acceleration  rates  in  the  outer 
region  of  the  separation  bubble  are  on  the  order  of  103  m/s2  [Ref.  34]. 

As  the  airfoil  rotates  upward  from  8  to  12°  angle  of  attack,  the  minimum 
vertical  distance  between  the  measurements  and  the  surface  of  the  airfoil  moves 
forward  from  12.3  to  9.3%  chord  on  the  airfoil  surface,  ranging  in  height  from 
1.58  to  0.58%  chord.  The  highest  velocities  in  the  flow  do  not  coincide  with 
these  locations  as  might  be  expected,  but  rather  are  forward  of  the  locations  and 
reside  in  the  region  over  the  separation  bubble. 

The  U-component  velocity  at  x/c  =  -0.100  shows  only  minor  effects  of  the 
downstream  oscillating  airfoil  in  Figure  5.23.  Again,  velocities  at  each  y/c 
location  are  offset  by  0.2  U^  for  presentation  purposes.  The  velocities  vary  from 
0.9  UM  to  1.1  UM  at  the  lowest  and  uppermost  heights  above  the  airfoil  surface. 
The  V-component  velocity  at  this  upstream  location  rises,  as  expected,  as  the 
airfoil  angle  of  attack  increases  from  8  to  12°.  V  velocities  throughout  the 
oscillation  cycle  are  about  one-third  of  the  U  component,  varying  from  0.2  U^  at 
the  upper  heights  to  as  much  as  0.351^  at  y/c  =  0.083.  The  measurements  show 
that  the  flow  remains  parallel  upstream  of  the  airfoil  and  that  the  streamlines  are 
curved  through  the  oscillation  cycle  as  noticed  in  the  V-component. 

As  the  flow  passes  over  the  airfoil  leading  edge  in  Figure  5.24  at  x/c  = 
0.017,  U-component  velocities  are  distinctly  different  at  the  three  lower  and  three 
upper  heights  above  the  airfoil  surface.  The  velocities  do  not  increase  as 
expected  as  the  airfoil  angle  of  attack  increases  to  the  top  of  the  oscillation  cycle 
like  at  the  upper  heights,  instead  the  velocities  show  prominent  decreases  at  the 
lowest  height.  At  y/c  =  0.083,  average  U  velocities  decrease  by  more  than  20% 
at  the  higher  angles  of  attack.  Velocities  vary  from  0.96 U^  at  a-  11.7°  during 
airfoil  upward  motion  to  1.35  LL  at  the  bottom  of  the  oscillation  cycle,  the 
opposite  of  what  one  might  expect.  The  decelerated  flow  is  a  result  of 
measurements  in  the  separation  bubble.  At  y/c  =  0.083,  the  reduced  velocity 
region  continues  from  a  -  9.3°  on  airfoil  upward  motion  to  a  -  11.0°  on  airfoil 
downward  motion  with  the  minimums  occurring  at  1 1°  <  a  <  12°.  The  maximum  U 
velocities  of  1.42 U^  occur  at  y/c  =  0.117  and  a=  10°,  and  also  at  y/c  =  0.133 
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and    a=   11.7°,  both  on  the  airfoil  upstroke  in  the  flow  region  above  the 
separation  bubble. 
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Figure  5.23.  Phase  Distribution  of  U-  and  V-Component  Velocities  at  x/c  =  -0.100. 
(Velocity  offset  by  0.2  U^  at  each  height.)   1^  =  0.3,  a  =  10° -2°  sin  cot,  k  =  0.1. 


The  V-component  velocities  at  upper  heights  are  much  like  that  forward  of 
the  airfoil  surface,  varying  from  0.2  to  0.25  U^  throughout  the  oscillation  cycle. 
However,  the  velocities  at  the  lowest  height  are  nearly  constant  at  all  angles  of 
attack  and  have  increased  to  0.4  U^.   The  higher  velocity  is  reasonable  because 
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the  airfoil  surface  at  the  x/c  =  0.017  measurement  location  presents  the  steepest 
angle  to  the  oncoming  flow,  except  at  the  airfoil  leading  edge.  V  velocities  in  the 
flow  field  reach  a  maximum  of  0.5  U^  at  a  =  12°  and  y/c  =  0.10. 
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Figure  5.24.  Phase  Distribution  of  U-  and  V-Component  Velocities  at  x/c  =  0.017. 
(Velocity  offset  by  0.21L  at  each  height.)   M^O.3,  a  =  10° -2°  sin  cot,  k  =  0.1. 

At  x/c  =  0.033  in  Figure  5.25,  the  U  velocity  decreases  in  an  expanding 
angle  of  attack  range  at  heights  closer  to  the  airfoil  during  airfoil  upward  motion 
even  though  the  measurements  are  0.533  mm,  or  2  to  3  probe  volumes  closer  to 
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the  airfoil  surface  than  at  the  previous  upstream  location.  U-component  velocities 
for  this  (x/c)  location  at  the  three  lowest  heights  were  on  the  order  of  20%  less 
than  those  shown  at  x/c  =  0.017.  The  velocities  at  y/c  =  0.083  drop  30%  from 
the  peak  velocity  of  1.25  UM  at  a  =  8°.  The  velocity  decrease  starts  with  the 
airfoil  upstroke  and  persists  until  a  =  10°  on  the  downstroke.  The  V-component 
velocities  at  the  upper  heights  followed  the  motion  of  the  airfoil  surface  as  before, 
but  at  the  two  lowest  heights  were  greatly  reduced  down  to  0.2  U^,  giving 
further  evidence  of  measurements  into  the  separation  bubble  region. 

d.         LD V  Velocity  Profiles  ( U  /  U  J 

Velocity  profiles  are  shown  for  angles  of  attack  in  1  degree  intervals  from 
8°  to  the  top  of  the  oscillation  cycle  and  then  back  down  to  9°.  Velocities  are 
offset  by  0.211^  at  each  angle  of  attack  for  better  presentation  and  are 
normalized  by  the  free  stream  velocity,  UM . 

Velocity  profiles  in  Figure  5.26  are  for  locations,  x/c  =  -0.10  and  x/c  =  0.00. 
At  x/c  =  -0.1,  U  velocities  ranged  from  0.96  to  1.08  U^.  At  the  lowest  measured 
height,  velocities  were  10%  lower  than  at  the  uppermost  height.  All  the  profiles 
overlie  each  other  if  the  offset  were  removed,  except  for  slight  velocity  reductions 
at  the  top  of  the  airfoil  downstroke  at  1 1  °  angle  of  attack. 

The  number  of  data  samples  at  x/c  =  -0.100  varied  from  137  on  airfoil 
downward  motion  at  a-  11.9°  to  383  at  a-  8.7°  during  airfoil  upward  motion. 
The  average  number  of  samples  is  274  if  evenly  distributed.  Generally,  less 
validated  samples  were  measured  near  the  top  of  the  oscillation  cycle.  However, 
there  were  always  more  than  50  to  give  an  acceptable  mean  value. 

In  general,  over  the  airfoil  leading  edge,  U  velocities  have  increased 
appreciably  (1.3 U^),  especially  at  the  lower  heights  in  closer  proximity  to  the 
airfoil  surface  such  that  the  profiles  are  now  more  uniform.  Velocities  at  mid- 
heights  and  12°  angle  of  attack  reach  as  high  as  1.411^. 
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Figure  5.25.  Phase  Distribution  of  U-  and  V-Component  Velocities  at  x/c  =  0.033. 
(Velocity  offset  by  0.2  UM  at  each  height.)   Meo=0.3,  a  =  10°  -2°  sin  cot,  k  =  0.1. 
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Figure  5.26.    U-Component  Velocity  Profiles  at  x/c  =  -0.10  and  x/c  =  0.00. 
(Velocity  offset  by  0.2  U^  at  each  AO A.)   MM=0.3,  a  =  10°  -2°  sin  cot,  k  =  0.1. 

U  velocity  profiles  for  x/c  =  0.017  and  0.033  are  shown  in  Figure  5.27. 
Obvious  in  both  plots  are  the  reduced  velocities  at  lower  heights,  especially  from 
10  to  12°  angle  of  attack  (up)  and  continuing  to  occur  during  the  airfoil 
downstroke  at  x/c  =  0.033.  Measured  U  velocities  span  a  larger  range  from  0.92 
to  1.45  UM  at  these  (x/c)  locations.  Velocity  profiles  with  slight  inflections  begin 
at  x/c  =  0.017  and  a  =  1 1°  (up),  and  downstream  at  x/c  =  0.033  inflections  are  first 
seen  at  a  =  10°.  With  the  sole  use  of  substantial  velocity  reductions  as  a  criteria, 
first  indications  of  the  separation  bubble  appear  at  x/c  =  0.017  and  10°  angle  of 
attack  (up). 
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The  number  of  validated  data  samples  was  widely  spread  at  x/c  =  0.017, 
varying  from  27  at  a  =  12°  and  y/c  =  0.100  to  474  at  a=  8.3°  and  y/c  =  0.083. 
Again,  the  larger  number  of  validated  samples  occurs  near  the  bottom  of  the 
oscillation  cycle  on  the  airfoil  upstroke. 

At  x/c  =  0.033,  U  velocity  gradients  are  large  from  9  to  11°  angle  of  attack 
on  the  upstroke,  and  the  profiles  cross  each  other  at  the  lower  heights  even  with 
the  added  offset.  At  this  downstream  location,  the  expanded  region  of  reduced 
velocities  indicates  that  a  larger  portion  of  the  flow  in  the  separation  bubble  is 
being  measured. 
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Figure  5.27.    U-Component  Velocity  Profiles  at  x/c  =  0.017  and  x/c  =  0.033. 
(Velocity  offset  by  0.2  UM  at  each  AOA.)   M^  =0.3,  a  =  10°  -2°  sin  cot,  k  =  0.1. 
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In  Figure  5.28,  velocities  are  given  for  x/c  =  0.067  and  0.083.  Again,  steep 
velocity  gradients  occurred  at  the  lower  heights  from  10  to  12°  angle  of  attack 
(up)  at  both  (x/c)  locations.  The  maximum  U  velocity  (1.5  U^)  over  the  airfoil 
surface  occurred  at  y/c  =  0.133  and  12°  angle  of  attack.  Indications  of  separation 
bubble  presence  are  slightly  greater  at  x/c  =  0.067  than  at  x/c  =  0.033.  However, 
velocities  at  the  lowest  height  (y/c  =  0.083)  and  angles  of  attack  of  1 1  (up)  and 
12°  begin  to  increase  slightly  from  those  at  the  upstream  location. 

At  x/c  =  0.083,  velocities  at  upper  heights  and  all  angles  of  attack  are 
slightly  lower  than  at  x/c  =  0.067  and  become  more  uniform  in  nature.  On  the 
other  hand,  velocities  at  the  lowest  height  and  all  angles  of  attack  are  greater  than 
at  x/c  =  0.067.  The  velocity  profile  at  10°  angle  of  attack  (up)  is  the  only  one 
that  resembles  the  profile  at  the  immediate  upstream  location  and  same  angle  of 
attack.  For  angles  of  attack  at  1 1  (up)  and  12°,  velocities  at  the  lowest  height  are 
greater  than  at  the  measurement  location  immediately  above,  showing  an  unusual 
shear  pattern  in  the  flow.  However,  notice  that  the  velocity  profile  at  11  °  angle 
of  attack  (down)  has  no  inverted  characteristic,  indicating  a  rapid  flow  recovery. 

The  velocity  profiles  shown  in  Figure  5.29  at  x/c  =  0.100  and  0.133 
indicate  that  the  measurements  are  approaching  the  downstream  edge  of  the 
separation  bubble.  At  x/c  =  0.100,  the  velocity  profiles  are  similar  to  those  shown 
previously  at  x/c  =  0.083  with  the  absence  of  the  inverted  "S"  profile 
characteristic. 

Velocity  profiles  at  x/c  =  0.133  are  similar  to  those  previously  shown  at  the 
airfoil  leading  edge.  No  dramatic  velocity  reductions  are  shown  at  any  angle  of 
attack.  The  right  edge  of  the  separation  bubble  appears  to  be  forward  of  13.3% 
chord  in  the  downstream  direction  if  based  on  the  velocity  profiles. 
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Figure  5.28.    U-Component  Velocity  Profiles  at  x/c  =  0.067  and  x/c  =  0.083. 
(Velocity  offset  by  0.21^  at  each  AOA.)   Moo=0.3,  cx  =  10°  -2°  sin  cot,  k  =  0.1. 
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Figure  5.29.    U-Component  Velocity  Profiles  at  x/c  =  0.100  and  x/c  =  0.133. 
(Velocity  offset  by  0.21^  at  each  AOA.)   M^O.3,  a  =  10°  -2°  sin  cot,  k  =  0.1. 

In  summary,  over  the  airfoil  leading  edge  surface,  the  velocities  greatly 
exceed  the  free  stream,  but  from  2  to  10%  chord  downstream,  the  profiles  show 
that  measurements  are  made  in  a  portion  of  the  separation  bubble. 


88 


e.  LDV  Velocity  Magnitude  Contours 

An  example  of  flowfield  velocity  magnitude  contours  is  given  for  a  =  9° 
(up/down)  in  Figure  5.30.  The  velocity  magnitude  range  in  the  flowfield  is 
indicated  in  the  legend  with  velocities  normalized  by  the  free  stream  and  shown 
in  0.06  U^  intervals. 

The  oncoming  flow  velocity  is  slightly  higher  during  the  airfoil  upstroke. 
Also,  the  maximum  flow  speed  (1.4411^)  region  is  more  extensive,  starting  slightly 
forward  of  the  airfoil  leading  edge  and  continuing  around  the  steep  velocity 
gradient  area  pointing  towards  the  location  of  the  separation  bubble.  It  is 
interesting  to  note  that  all  the  measured  velocities  shown  in  the  flow  field  are 
greater  than  the  free  stream.  For  airfoil  oscillations  in  both  directions,  a 
decelerated  flow  region  exists  over  the  airfoil  leading  edge  with  the  minimum 
velocity  contour  in  each  case  being  1.02  U^.  The  lower  velocity  region  is  wider 
and  extends  to  greater  heights  above  the  airfoil  surface  during  the  airfoil 
upstroke.  For  example,  velocity  contour  number  three  (1.0211^)  during  airfoil 
upward  motion  exceeds  the  lateral  boundaries  of  contour  number  four  (1.0811^) 
on  airfoil  downward  motion. 

f.  LDV  Vorticity  Magnitude  Contours 

The  vorticity  was  obtained  by  first  interpolating  the  velocity  component 
data  and  mapping  the  data  with  a  spline  curve  fit.  Then  the  z-component  of 
vorticity  was  calculated  from  the  velocity  gradients  of  the  fitted  curves.  The 
vorticity  was  normalized  by  c/LL  and  is  presented  in  the  next  three  figures. 

Flow  field  vorticity  is  given  for  angles  of  attack  of  8  and  12°  in  Figure 
5.31.  The  vorticity  levels  in  the  flow  over  the  airfoil  at  8°  angle  of  attack  are  low 
compared  with  those  at  the  top  of  the  oscillation  cycle.  Only  one  region  of 
counterclockwise  vorticity  occurs  in  the  flowfields  at  the  lowest  measurement 
location  and  aft  of  8%  chord  near  the  downstream  edge  of  the  separation  bubble. 
Clockwise  vorticity  (-24  units  maximum)  exists  only  near  the  airfoil  leading  edge 
from  0.00  <  x/c  <  0.08  and  at  the  lowest  heights  above  the  separation  bubble  for 
a  =  8°.     At  a  =  12°,  the  higher  levels  of  clockwise  vorticity  (-30  units)    occur 

89 


-0.06 


NO  CONTOUR  LINES  FOUND 


J L 


Figure  5.30.     Velocity  Magnitude  IVI  Contours  Near  Airfoil  Leading  Edge. 
Moo=0.3,  a=  10o-2°sina)t  =  9°(up/down),  k  =  0.1. 
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between  0.01  <  x/c  <  0.03,  extending  to  1.7%  chord  above  the  airfoil  surface.  In 
addition,  two  pockets  of  clockwise  vorticity  (-30  units)  occur  at  4.5  and  8% 
chord  downstream,  both  at  5.5%  chord  above  the  airfoil  surface. 
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Figure  5.31.    Vorticity  Contours  (z-Component)  Near  the  Airfoil  Leading  Edge. 
1^  =  0.3,  cc=  100-2°sincot  =  80(bottom)andl2°(top),  k  =  0.1. 
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In  Figure  5.32,  vorticity  in  the  flow  field  is  shown  for  11°  angle  of  attack 
during  airfoil  upward  motion  and  11.53°  angle  "of  attack  on  the  airfoil 
downstroke.  The  vortical  flow  at  1 1  °  angle  of  attack  (up)  is  much  like  that  given 
for  12°  angle  of  attack  in  the  previous  figure  with  a  maximum  clockwise  vorticity 
level  of  -30  units.  Although,  no  counterclockwise  vorticity  is  shown  aft  of  the 
separation  bubble  location  as  at  the  top  of  the  oscillation  cycle.  In  each  case  a 
very  small  region  of  clockwise  vorticity  (+6.0  units)  is  indicated  just  forward  of 
the  airfoil  leading  edge  between  two  clockwise  vorticity  regions.  On  the  airfoil 
downstroke  at  11.53°  angle  of  attack,  the  clockwise  vorticity  over  the  airfoil 
leading  at  2%  chord  downstream  has  increased  by  at  least  6  units  to  -36  units 
(shown  earlier  at  -40  units,  [Ref.  41]),  greater  than  previously  measured  at  1 1  (up) 
and  12°  angle  of  attack.  The  two  pockets  of  clockwise  vorticity  discussed  for 
a  =  12°  have  now  coalesced  into  one  region  of  increased  vorticity  and  extends  to 
slightly  higher  than  y/c  =  0.100  or  4.0%  chord  above  the  airfoil  surface. 

At  10.68°  angle  of  attack  (down)  in  Figure  5.33,  the  maximum  clockwise 
vorticity  has  decreased  to  -24  units.  The  concentrated  areas  of  high  vorticity  that 
existed  at  a  =  11.53°  (down)  are  now  seen  to  move  downstream.  At  a  =  10°,  the 
vorticity  levels  increase  again  in  the  region  over  which  the  separation  bubble 
develops  again  over  the  airfoil  leading  edge.  Flow  reattachment  in  the  PDI 
measurements  was  seen  to  begin  already  at  a  =  10.5°  (down),  and  the  increased 
vorticity  indicated  by  the  LDV  measurements  above  what  was  measured  at  the 
higher  angle  of  attack  indicates  that  the  vorticity  was  partially  shed  between 
11.53  and  10°  angle  of  attack.  The  vorticity  distributions  corroborate  the 
occurrence  of  light  dynamic  stall  that  was  shown  earlier  in  the  PDI  images.  The 
vorticity  apparently  can  build  up  only  to  certain  levels  before  it  must  be  shed  by  a 
means  other  than  diffusion  through  the  boundary  layer.  By  having  the  airfoil 
oscillate  near  the  stall  angle  of  attack,  the  vorticity  in  the  flow  reached  levels  that 
require  the  vorticity  balance  to  be  preserved  as  the  flow  attempts  to  adapt  to  the 
decreasing  angle  of  attack  and  was  shed  on  the  airfoil  downstroke.  The  vorticity 
shedding  process  occurs  over  a  small  range  of  angle  of  attack.  This  surprising 
result  of  light  dynamic  stall  is  documented  by  both  the  PDI  and  LDV 
measurements  at  an  angle  of  attack  below  the  static  stall  on  the  airfoil 
downstroke,  enhancing  the  confidence  in  the  measured  data. 
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Figure  5.32.   Vorticity  Contours  (z-Component)  Near  the  Airfoil  Leading  Edge. 
Moo=0.3,  a=  10°-20sincot  =  ll0(up)andll.53°(down),  k  =  0.1. 
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Figure  5.33.    Vorticity  Contours  (z-Component)  Near  the  Airfoil  Leading  Edge. 
MM=0.3,  a=  10°-2osina)t  =  10.68o(down)andl0°(down),  k  =  0.1. 


94 


Previously,  it  was  thought  that  the  separation  bubble  remained  on  the 
airfoil  surface,  expanding  and  contracting  in  size  during  the  oscillation  cycle  for 
the  low-amplitude  test  case  [Ref.  65].  Both  the  PDI  and  LDV  test  results 
confirmed  the  incipient  stall,  as  previously  reported  in  Reference  41.  The 
measurements  show  that  the  vertical  extent  of  the  viscous  zone  remains  thin  and 
on  the  order  of  the  airfoil  thickness,  giving  valuable  data  for  making  comparisons 
with  computational  methods.  Thus,  the  capabilities  of  the  computational  methods 
to  reproduce  the  experimental  results  are  evaluated,  using  the  N-S  approach  as 
the  main  method  of  analysis  supplemented  by  potential  flow  and  boundary  layer 
methods  as  discussed  in  Chapter  VI.  The  computed  and  measured  comparisons 
are  reserved  for  Chapter  VII. 
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VL       COMPUTED  RESULTS 

In  this  chapter,  computational  results  are  presented  using  the  potential 
flow,  boundary  layer,  and  N-S  analysis  methods  described  in  Chapter  IV.  The 
potential  flow  method  provides  a  very  rapid  means  of  visualizing  the  effect  of 
angle  of  attack,  amplitude,  and  frequency  of  airfoil  oscillation  on  the  pressure 
distribution  for  both  the  lower  and  upper  airfoil  surfaces.  In  particular,  it  permits 
quantification  of  the  phase  lag  between  the  airfoil  motion  and  the  pressures 
induced  at  various  points  on  the  airfoil.  The  combination  of  the  unsteady 
potential  flow  code  (panel)  with  the  boundary-layer  code  provides  another 
efficient  and  rapid  tool  to  visualize  the  response  of  the  boundary  layer  to  changes 
in  angle  of  attack,  amplitude,  and  frequency  of  oscillation.  This  combined 
panel/boundary  layer  approach  is,  of  course,  limited  to  the  analysis  of  attached 
flows  only  because  the  boundary  layer  is  computed  in  the  direct  mode. 
Furthermore,  the  calculations  are  limited  to  relatively  low  reduced  frequencies  of 
oscillation  because  the  boundary  layer  code  is  based  on  the  steady  boundary 
layer  equations.  This  latter  limitation  is  likely  to  have  only  a  minor  effect  on  the 
accuracy  of  the  results  for  reduced  frequencies  less  than  one  because  the 
boundary  layer  lag  effects  are  known  to  be  much  less  than  the  pressure  lag 
effects.  The  thin-layer  N-S  code  is  used  to  analyze  the  development  of 
separations  bubbles  and  the  formation  and  shedding  of  the  dynamic  stall  vortex. 

A.        FLOW  OVER  STEADY  NACA  0012  AIRFOIL 

1.         Pressure  Distributions 

The  experiments  have  shown  that  steady-state  stall  occurs  at  an  angle  of 
attack  slightly  above  12°  at  a  free-stream  Mach  number  of  0.3  and  a  Reynolds 
number  of  540,000.  Therefore,  it  is  instructive  to  study  the  change  in  pressure 
distribution  as  the  angle  of  attack  is  increased  from  0°  to  12°.  This  is  shown  in 
Figures  6.1  and  6.2  for  the  first  quarter  chord  of  the  airfoil,  using  UPOT.  It  is  seen 
that  the  leading-edge  suction  pressures  increase  dramatically  with  increasing 
angle  of  attack  while  the  suction  peaks  move  slightly  forward.  The  stagnation 
point  moves  downstream  on  the  lower  surface  from  1  to  4%  chord.  The  maximum 
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adverse  pressure  coefficient  gradient  reaches  a  value  of  40  at  a  =  8°  and  115  at 
a  =  12°. 
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Figure  6.1.   Surface  Pressure  Distribution  Near  the  Airfoil  Leading  Edge,  UPOT. 
a  =  0°,  4°,  8°  and  12°. 
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Figure  6.2.   Surface  Pressure  Distribution  Near  the  Airfoil  Leading  Edge,  UPOT. 
a  =  2°,  6°  and  10°. 
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In  Figure  6.3,  the  UPOT-predicted  pressure  distributions  are  compared  with 
the  N-S  code  predictions.  In  the  N-S  calculations,  the  effect  of  boundary  layer 
transition  from  laminar  to  turbulent  flow  is  taken  into  account  by  specifying  the 
transition  onset  location  and  by  computing  the  transition  length  using  the  Chen- 
Thyson  transition  model.   The  transition  value  Gy    in  the  Chen-Thyson  model 

was  chosen  to  be  200,  the  basis  of  which  is  discussed  in  detail  in  Chapter  VII.  As 
expected,  complete  agreement  between  the  two  codes  is  achieved  at  oc  =  0°, 
whereas  at  a  =  6°,  a  slight  deviation  occurs  near  the  suction  peak.  The  N-S 
predicted  suction  is  slightly  lower  than  the  UPOT  prediction.  However,  at 
a  =  12°,  the  suction  predictions  are  distinctly  different.  The  N-S  computation 
shows  a  much  lower  suction  pressure  and  a  pressure  plateau  typical  for  the 
formation  of  a  small  separation  bubble;  whereas,  of  course,  the  UPOT  code  cannot 
account  for  this  viscous-dominated  phenomenon.  It  is  also  important  to  note  that 
a  purely  laminar  N-S  calculation  fails  to  converge,  while  purely  turbulent  N-S 
calculations  produce  excessive  suction  peaks  much  like  the  UPOT-predicted 
peaks  for  flow  conditions  in  the  transitional  Reynolds  number  range. 
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Figure  6.3.  Surface  Pressure  Distribution  Near  Airfoil  Leading  Edge,  UPOT  and 
N-S  Codes.  N-S  Analysis,  M^  =  0.3,  a  =  0°,  6°  and  12°,  Re  =  540,000,  Glu  =  200. 
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Figure  6.4  shows  that  the  N-S  code  first  predicts  a  small  separation  bubble 
at  a  =  8°,  which  is  in  accordance  with  the  measurements.  At  a  =  10°,  the 
predicted  separation  bubble  is  significantly  bigger  and  it  is  located  further 
forward.  At  a  =  12°,  the  N-S  calculations  showed  remarkable  sensitivity  to  the 
choice  of  the  transition  onset  location  and  the  transition  length.  Figure  6.4 
shows  two  calculations  for  a  =  12°,  one  with  a  separation  bubble  similar  to  the 
bubble  obtained  at  a  =  10°,  the  other  one  showing  a  fully  stalled  flow.  These  two 
dramatically  different  solutions  were  obtained  by  specifying  a  slightly  earlier- 
transition  onset,  moving  it  from  1.8%  chord  to  1.2%  chord,  and  by  specifying  a 
slightly  longer  transition  length  by  choosing  a  GYtr  value  of  200  rather  than  100. 

The  details  of  transition  onset  and  length  effects  are  covered  in  Chapter  VII.  Aft 
of  the  transition  zone,  the  Baldwin-Barth  turbulence  model  is  used  in  the  N-S 
code  for  all  the  results  given  in  this  thesis. 
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Figure  6.4.  Surface  Pressure  Distribution  Near  Airfoil  Leading  Edge,  N-S  Code, 
Transition  Onset  Locations  Specified.  Mx  =  0.3,  a  =  8°,  10°  and  12°,  Re  =  540,000, 
Gv    =  200,  or  100  for  (x  /  c)  transition  =  0.0177 . 
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2.  Skin-Friction  Coefficient  Distribution  on  Suction  Surface 

The  panel/boundary-layer  code  analysis  used   100  and  214  panels  to 
predict  the  skin-friction  coefficient  distribution  on  the  suction  surface.    The 
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program  uses  the  Chen-Thyson  transition  model  and  the  Cebeci-Smith  eddy- 
viscosity  model.  The  100-panel  results  are  shown  in  Figure  6.5,  and  the  214  panel 
results  are  shown  in  Figure  6.6.  At  the  lower  angles  of  attack  of  0°,  2°,  and  4°, 
differences  are  noticed  in  the  slightly  smoother  curves  for  the  denser  panel  (214) 
arrangement  and  in  the  transition  onset  locations.  Initial  transition  onset 
locations  were  estimated  by  using  Michel's  criterion,  and  in  all  cases  Michel's 
onset  location  was  aft  of  the  specified  transition  onset  locations  used  for  the  final 
converged  solutions  in  the  boundary  layer  code.  The  214-panel  solution  shows  a 
forward  progression  of  the  transition  onset  location  for  increasing  angles  of 
attack  up  to  8°,  which  is  in  agreement  with  experimental  observations.  Transition 
onset  locations  for  the  100-panel  solution  are  similar,  except  from  a  =  6°  to  a  =  8° 
where  the  converged  solution  incorrectly  predicted  onset  to  move  aft  on  the 
airfoil  surface. 

Flow  separation  is  indicated  when  Cf  passes  through  0  to  negative  values. 
Negative  skin-friction  values  occur  only  for  8  °  angle  of  attack  in  a  small  region 
near  the  leading  edge  of  the  airfoil,  which  is  in  the  same  general  area  as  the 
measurements.  Specifically,  the  boundary-layer  equations  predict  negative  CfS  at 
only  three  airfoil  surface  locations  from  x/c  =  0.017  to  0.021,  which  is  forward  of 
the  measured  results.  Figures  6.5  and  6.6  show  that  transition  onset  was  set  near 
locations  of  minimum  Cf  and  then  the  flow  immediately  responds  to  the  transition 
model.  The  rapid  increase  in  positive  Cf  values  after  the  small  negative  region 
indicates  that  the  laminar  flow  has  transitioned  to  turbulent  flow  and  the  flow  has 
reattached  downstream  of  the  small  leading-edge  separation  bubble. 

It  can  also  be  seen  that  the  well  known  breakdown  of  the  boundary-layer 
equations  at  the  flow  separation  point  does  not  cause  an  immediate  divergence  of 
the  finite-difference  solution  of  these  equations.  Instead  the  code  is  able  "to 
predict"  a  small  region  of  negative  skin  friction  at  a  =  8°,  using  either  100  or  214 
panels.  Hence,  the  code  can  be  used  to  identify  the  angle  of  attack  at  which  a 
tiny  separation  bubble  will  occur.  Obviously,  the  direct  boundary-layer 
computations  are  expected  to  diverge  as  the  region  of  negative  skin  friction 
increases  and,  indeed,  no  flow  solutions  were  found  for  either  panel  arrangement 
at  angles  of  attack  above  8°.     Table  6.1  gives  an  indication  of  suggested 
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transition  locations  and  convergence  ranges  for  the  panel/boundary-layer  code, 
using  either  100  or  214  panels. 
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Figure  6.5.    Suction  Surface  Skin-Friction  Distributions,  Panel/Boundary- Layer 
Code  (100  Panels),  Transition  Onset  Specified.  Re  =  540,000. 


u 


c 

— 
u 
o 
U 

c 
_o 

o 

£ 


on 


0.006 


0.004 


0.002 


0.000 


-0.002 


pil 

AOA.  (x/c)transition 

0  degrees.  0.515 

2  degrees,  0.350 

-.-  4  degrees.  0.177 

6  degrees.  0.023 

8  degrees.  0.022 

t!l  \                                     i 

f! '.  \ 

\\\\                \ 

h\k        *"" -':-^'       " "" "  '"'- 

i:'    '.,'  \                                   /~-- -U;;-- 

■t:          V                  :     /*»>                   "^-i^                              > 

ii-             ■■*.^   :/            --..     y ~^__/ 

!' 
ii 

— 

**  "-T           *  *  "  "t  *  » 

.....yf. r...^.. ....... ...  ( 

/                     W  V 

;          V. 

\  \ 
\: 

1 
j 

:                        j 

j                i 

t 

1 

■1 

I  ' 

ii 

li 
Ii 

i! 
2 

0.0 


0.2 


0.4 


0.6 


0.8 


1.0 


Chordwise  Location,  x/c 


Figure  6.6.    Suction  Surface  Skin-Friction  Distributions,  Panel/Boundary-Layer 
Code  (214  Panels),  Transition  Onset  Specified,  Re  =  540,000. 
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Previously,  Nowak  [Ref.  66]  used  100  panels  to  investigate  the  effect  of 
Reynolds  number.  Consistent  with  the  experimental  evidence,  the  predictions 
showed  that  as  the  Reynolds  number  increased,  flow  separation  was  less  likely  to 
happen. 

Table  6.1.  Panel/Boundary-Layer  Code  Transition  Onset  Criteria. 


Angle  ol  attack. 
Degrees 

0 

2 

4 

6 

8 

Number  of  panels 

100 

214 

100 

214 

100 

214 

100 

214 

100 

214 

Stagnation  panel 
number 

51 

101 

50 

96 

49 

91 

48 

87 

47 

82 

Michel's  transition 
onset  prediction. 

Upper  (x/c) 

0.598 

0.560 

0.380 

0.393 

0.244 

0.258 

0.099 

0.112 

0.164 

0.057 

Lower  (x/c) 

0.598 

0.560 

0.769 

0.744 

0.839 

0,836 

0.898 

0.877 

0.909 

0.897 

Convergence 
range 

Upper  (x/c) 

0.57- 
0.61 

0.49- 
0.53 

0.33- 
0.37 

0.33- 
0.35 

0.17- 
0-19 

0.17- 
0.20 

0.025- 
0.029 

0.019  - 
0.024 

0.029  - 
0.041 

0.022  - 
0.023 

Transition  onset 
specified,  (best 
convergence) 

Upper  (x/c) 

0.59 

0.515 

0.36 

0.35 

0.18 

0.177 

0.028 

0.023 

0.035 

0.022 

N-S  solutions  for  the  same  of  attack  range  are  shown  in  Figure  6.7.  The 
skin-friction  curves  show  the  same  trends  that  were  obtained  with  the  boundary- 
layer  solutions,  including  negative  values  only  for  the  8  °  case.  Transition  onset 
locations  for  0  and  4°  were  set  to  values  used  for  the  boundary-layer  solutions 
with  214  panels.  For  other  angles  of  attack,  the  N-S  transition  onset  values  were 
specified  such  that  the  solutions  correlated  best  with  the  measured  results.  The 
initial  time  steps  of  the  N-S  solution  were  stabilized  by  using  laminar  flow 
calculations.  Thereafter,  the  Chen-Thyson  transition  and  Baldwin-Barth 
turbulence  models  are  initiated  at  100  to  300  time  steps,  and  the  predicted  onset 
of  flow  transition  moves  dramatically  forward  on  the  airfoil  surface.  However, 
Michel's  criterion  for  transition  onset  gave  a  separation  bubble  that  was  much 
smaller  than  the  measured  ones.   Therefore,  transition  onset  was  specified  aft  of 
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Michel's  criterion  to  give  predicted  results  that  agreed  more  closely  with  the 
measurements. 

As  noted  before,  a  Gy    value  of  200  was  used  to  set  the  transition  region 

length.  Minimum  CfS  move  forward,  remain  positive  and  increase  with  increasing 
angles  of  incidence  for  both  the  boundary  layer  and  N-S  solutions  until  both 
predict  a  negative  Cf  region  near  the  airfoil  leading  edge  for  a  =  8°.  Obvious 
differences  between  the  two  solutions  include  a  larger  separation  bubble 
predicted  by  the  N-S  solution  at  a  =  8°.  Negative  CfS  are  predicted  at  nine 
surface  locations  from  x/c  =  0.028  to  0.058  for  the  N-S  solutions,  which  indicates 
a  laminar  separation  bubble  more  in  line  with  the  measured  results  than  that 
predicted  by  the  boundary-layer  code.  Higher  CfS  are  predicted  along  the  airfoil 
upper  surface  by  the  N-S  code  downstream  of  transition  onset  and  flow 
reattachment  that  are  in  part  attributable  to  differences  in  turbulent  modeling 
(Baldwin-Barth  vs.  Cebeci-Smith)  and  compressible  flow  considerations  that  are 
not  considered  in  the  direct  boundary-layer  analysis. 
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Figure  6.7.    Suction  Surface  Skin-Friction  Distributions,  N-S  Code,  275x81, 
Transition  Onset  Locations  Specified,  MM  =  0.3,  Re  =  540,000,  GYtr  =  200. 
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3.  Displacement  Thickness  Distribution  on  Suction  Surface 

The  boundary-layer  displacement  thickness  distributions  on  the  suction 
surface  as  predicted  by  the  boundary  layer  (214  panels)  and  N-S  codes  are  shown 
in  Figures  6.8  and  6.9.  Both  computations  yield  the  well  known  dip  in 
displacement  thickness  due  to  boundary  layer  transition.  The  movement  of  the 
transition  region  toward  the  leading  edge  is  clearly  visible  as  the  angle  of  attack  is 
increased  from  0  to  8°.  Also,  it  is  seen  that  the  computations  agree  reasonably 
well. 
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Figure  6.8.  Suction  Surface  Displacement  Thickness  Distributions,  Panel/ 
Boundary-Layer  Code  (214  Panels),  Transition  Onset  Locations  Specified. 
Re  =  540,000. 
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Figure  6.9.    Suction  Surface  Displacement  Thickness  Distributions,  N-S  Code, 
275x81,    Transition    Onset   Locations    Specified.    Moo=0.3,     Re  =  540,000, 


Gv    =  200. 


4.  Boundary-Layer  Velocity  Profiles  on  Suction  Surface 

The  flow  behavior  can  be  understood  in  more  detail  by  examining  the 
change  in  boundary-layer  profiles  as  a  function  of  chordwise  locations  on  the 
suction  surface.  Therefore,  the  profiles  computed  with  the  214-panel  boundary- 
layer  code  are  shown  for  0,  2,  4,  6,  and  8°  angle  of  attack  in  Figures  6.10  through 
6.14.  It  is  readily  seen  that,  in  each  case,  the  laminar  flow  profiles  near  the  leading 
edge  change  to  profiles  with  inflection  points  indicative  of  flow  into  an  adverse 
pressure  gradient  region.  This  is  consistent  with  the  previously  shown  skin 
friction  and  displacement  thickness  distributions.  Laminar  flow  separation  is 
averted  by  virtue  of  the  fact  that  the  flow  starts  to  transition  to  turbulent  flow 
and  thereby  retains  the  ability  to  remain  attached.  This  transition  occurs  at 
approximately  50%  chord  for  the  zero  angle  of  attack  flow  and  rapidly  moves 
toward  the  leading  edge  as  the  angle  of  attack  is  increased.     Note  that  the 


106 


velocity  profiles  are  non-dimensionalized  by  the  local  external  inviscid  flow 
velocity,  Ue,  giving  an  indication  of  the  growth  of"  the  boundary  layer  in  the 
streamwise  direction.  In  the  Figures,  UG  =  0.995  xUe. 
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Figure  6.10.    Boundary-Layer  Velocity  Profiles  (Panel/Boundary-Layer  Code 
214  Panels).  a  =  0°,  (x/c)tr=  0.515,  Re  =  540,000. 
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Figure  6.11.    Boundary-Layer  Velocity  Profiles  (Panel/Boundary-Layer  Code 
214  Panels),  a  =  2°,  (x/c)tr  =  0.35,  Re  =  540,000. 
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Figure  6.12.    Boundary-Layer  Velocity  Profiles  (Panel/Boundary-Layer  Code  - 
214  Panels).  a  =  4°,  (x/c)tr=  0.177,  Re  =  540,000. 

In  Figures  6.13  and  6.14,  velocity  profiles  are  shown  for  angles  of  attack  of 
6  and  8°,  which  are  much  farther  forward  than  those  shown  in  the  previous 
figures.  Solution  convergence  at  the  higher  angles  of  attack  was  very  sensitive 
to  flow  transition  onset  locations,  which  was  specified  at  x/c  =  0.023  for  a  =  6° 
and  at  x/c  =  0.022  for  a  =  8°.  Again,  velocity  profiles  with  inflection  points  are 
shown  forward  of  the  specified  transition  onset  locations,  but  only  slightly.  The 
slight  forward  movement  of  the  transition  onset  location,  and  no  indication  of  full 
velocity  profiles  aft  of  flow  transition  in  the  turbulent  flow  region,  points  toward 
prediction  of  an  impending  flow  separation.  Flow  separation  is  predicted  by  the 
velocity  profiles  at  x/c  =  0.0298  and  0.0413  at  8°  angle  of  attack.  Thus,  the 
panel/boundary-layer  code  predicts  flow  separation  at  a  =  8°  with  no  solutions 
found  above  8°. 
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Figure  6.13.  Boundary-Layer  Velocity  Profile  Predictions  (Panel/Boundary-Layer 
Code  -  214  Panels),  a  =  6°,  (x/c)tr  =  0.023,  Re  =  540,000. 
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Figure  6.14.    Boundary  Layer  Velocity  Profiles  (Panel/Boundary-Layer  Code 
214  Panels),  a  =  8°,  (x/c)tr  =  0.022,  Re  =  540,000. 
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Figure  6.15  gives  the  velocity  profiles  for  both  the  100-  and  214-panel 
solutions  at  6°  angle  of  attack  and  the  same  x/c  locations  as  previously  given  in 
Figure  6.13.  Only  minor  differences  occur  in  the  computed  velocity  profiles. 
Panel  densities  up  to  400  were  used  to  establish  grid  independence  and  no 
differences  in  velocity  profiles  were  found.  However,  panel  numbers  greater  than 
400  caused  a  breakdown  in  the  solution  with  unrealistic  multiple  dips  in  the  skin- 
friction  coefficient  and  multiple  vortices  predicted  for  a  =  6  and  8°. 
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Figure  6.15.    Boundary-Layer  Velocity  Profiles  (Panel/Boundary-Layer  Code  - 
100  and  214  Panels),  a  =  6°,  (x/c)tr  =  0.023,  Re  =  540,000. 

An  advantage  of  the  N-S  analysis  method  over  the  boundary-layer  code  is 
the  capability  to  compute  extensive  regions  of  flow  separation.  The  next  series  of 
figures  give  the  N-S  predictions  for  angles  of  attack  from  0  to  12°  angle  of  attack. 
Here  the  velocity  profiles  are  non-dimensionalized  with  the  free-stream  speed. 
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Figure  6.16  shows  the  predicted  velocity  profiles  for  a  =  0°  with  transition 
set  to  x/c  =  0.515,  just  as  in  the  boundary  layer  code.  Velocities  are  seen  to 
exceed  the  free  stream  by  up  to  18%  forward  of  the  adverse  pressure  gradient 
and  the  transition  onset  location. 
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Figure  6.16.    Boundary-Layer  Velocity  Profiles  (N-S  Code,  275x81).  MM  =0.3, 
a  =  0°,  Re  =  540,000. 

Next,  velocity  profiles  are  shown  at  similar  x/c  locations  for  a  =  4°  angle  of 
attack  in  Figure  6.17.  Transition  onset  moved  forward  and  again  was  specified  at 
the  same  location  (x/c  =  0.18)  as  the  boundary-layer  code  solution.  No  reverse 
flow  regions  are  predicted  and  much  higher  velocities  (U  /  UM  =  1.33)  are  reached 
forward  of  the  adverse  pressure  gradient  region  because  of  the  higher  angle  of 
attack.  Two  velocity  profiles  show  points  of  inflection.  Fully  turbulent  velocity 
profiles  are  evident  aft  of  the  flow  transition  region. 
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Figure  6.17.    Boundary-Layer  Velocity  Profiles  (N-S  Code,  275x81).  MM  =0.3, 
a  =  4°,  Re  =  540, 000. 

The  next  three  Figures  give  the  N-S  computed  velocity  profiles  for 
a=  8°,  10°,  and  12°  near  the  airfoil  leading  edge  because  transition  onset  has 
moved  to  the  first  4%  chord.  As  shown  in  Figure  6.18  for  a  =  8°,  the  local 
velocity  is  predicted  to  exceed  the  free-stream  velocity  by  up  to  77%.  Transition 
onset  was  set  at  x/c  =  0.035  to  give  a  N-S  solution  that  most  closely  matched  the 
experimental  measurements  (Cps  from  PDI  and  velocity  profiles  from  LDV). 
Points  of  inflection  are  shown  in  five  profiles  around  the  specified  transition  onset 
point.  Flow  separation  is  evident  at  x/c  =  0.0298,  0.0413,  and  0.0547  with  an 
extensive  reverse  flow  region  (U/U^  =  -0.09)  shown  at  x/c  =  0.0413  and  y/c  = 
0.0002.  Thus,  the  N-S  code  predicts  a  laminar  flow  separation,  which  is  evident  at 
x/c  =  0.0298  and  forward  of  the  specified  transition  onset  location. 
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Figure  6.18.   Boundary-Layer  Velocity  Profiles  (N-S  Code,  275x81).  MM  =0.3, 
a  =  8°,  Re  =  540,000. 

Figure  6.19  gives  the  velocity  profiles  for  a  =  10°.  Transition  onset  moves 
closer  to  the  leading  edge  of  the  airfoil  at  x/c  =  0.0225.  All  profiles  in  the  first  7% 
chord  have  points  of  inflection,  except  the  first  two  near  the  airfoil  leading  edge. 
Again,  laminar  separation  occurs  just  prior  to  the  specified  transition  onset 
location.  By  examining  all  the  velocity  profiles  in  the  vicinity  of  separation, 
laminar  flow  separation  is  first  predicted  at  x/c  =  0.0117,  and  separated  flow 
velocity  profiles  are  shown  to  7%  chord.  The  predicted  separation  bubble 
extends  from  0.01 17  <  x/c  <  0.08  along  the  airfoil  surface  and  0.0  <  y/c  <  0.00085 
for  the  reverse  flow  region  above  the  airfoil  surface.  Local  velocities  up  to 
1.95x11^  are  shown  around  the  indicated  separation  bubble. 
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Figure  6.19.   Boundary-Layer  Velocity  Profiles  (N-S  Code,  275x81).  MM  =0.3, 
a  =  10°,  Re  =  540,000. 

Figure  6.20  gives  the  velocity  profiles  for  a  =  12°.  Transition  onset  moved 
forward  to  x/c  =  0.0177  and  Gv  was  set  to  100.  At  this  angle  of  attack  near  stall 
conditions,  the  local  flow  reaches  values  of  2.22  xUM.  Flow  separation  begins  at 
x/c  =  0.0085  and  appreciable  reverse  flow  extends  to  1%  chord,  indicating  a  flow 
on  the  verge  of  stall.  A  maximum  negative  velocity  ( U  /  UM  =  -  0.48)  is  predicted 
at  x/c  =  0.0201,  which  again  is  aft  of  the  specified  transition  onset  location.  The 
extent  of  the  predicted  separation  bubble  ranges  from  0.0085  <  x/c  <  0.10  in  the 
downstream  direction  and  from  0.0  <  y/c  <  0.0009  above  the  airfoil  surface. 
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Figure  6.20.    Boundary-Layer  Velocity  Profiles  (N-S  Code,  275x81).  M^  =0.3, 
a  =  12°,  Re  =  540,000. 

It  is  instructive  to  compare  the  boundary-layer  and  N-S  code  velocity 
profiles  directly.  To  this  end,  velocities  are  normalized  by  using  the  free-stream 
velocity  as  the  normalization  factor  in  both  cases.  For  these  runs,  the  boundary 
layer  code  included  a  test  for  laminar  separation.  If  laminar  separation  is  predicted 
before  Michel's  transition  prediction,  then  flow  transition  onset  is  specified  at  the 
panel  forward  of  laminar  separation.  The  aforementioned  method  precludes  the 
computation  of  a  separation  bubble,  which  is  not  done  well  by  the  boundary- 
layer  code.  The  comparisons  are  based  on  the  choices  for  transition  onset 
specified  in  Figures  6.10  and  6.12  for  the  boundary  layer  calculations  and  Figures 
6.16  and  6.17  for  the  N-S  calculations. 

Figures  6.21  and  6.22  show  both  predicted  velocity  profiles  at  a  =  0°  and 
a  =  4°.  The  boundary  layer  velocity  profiles  are  output  until  minimal  changes  are 
seen  in  the  outer  velocity  ratio.    Excellent  agreement  is  shown  between  the 
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calculations,  keeping  in  mind  that  the  boundary-layer  code  is  based  on  the 
incompressible  flow  equations,  whereas  the  N-S  calculations  consider  a  M^  =  0.3 
compressible  flow. 
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Figure  6.21.    Boundary-Layer  Velocity  Profile  Comparisons  (N-S  and  Panel/ 
Boundary-Layer  Codes),  a  =  0°,  Re  =  540,000;  N  -  S  Analysis,  MM  =  0.3. 
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Figure  6.22.  Boundary-Layer  Velocity  Profile  Comparisons  (N-S  and  Panel/ 
Boundary-Layer  Codes),  a  =  4°,  Re  =  540,000;  N  -  S  Analysis,  IVL  =  0.3. 
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For  a  =  8°,  the  boundary-layer  solution  predicts  a  laminar-flow  separation 
forward  of  Michel's  transition  prediction,  at  x/c  =  0.015,  which  is  markedly 
forward  of  the  N-S  laminar  separation  prediction  at  x/c  =  0.0298.  Thus,  no 
comparisons  are  shown  for  a  =  8°  because  of  the  early  laminar  flow  separation 
predicted  by  the  boundary  layer  code. 

B.        OSCILLATORY  FLOW  OVER  NACA  0012  AIRFOIL 

For  the  analysis  of  the  oscillatory  pitch  tests  of  the  NACA  0012  airfoil 
described  in  Chapter  V,  it  is  useful  to  develop  an  understanding  of  the  inviscid 
flow  effects  prior  to  delving  into  the  unsteady  viscous  flow  phenomena. 
Therefore,  the  inviscid,  incompressible  UPOT  code  is  first  used  to  study  the 
pressure  lag  effects  that  are  induced  by  pitch  oscillations  about  the  quarter  chord 
point  of  the  airfoil.  UPOT  results  are  presented  for  oscillation  amplitudes  of  10 
and  2°  about  a  mean  angle  of  attack  of  10°.  The  oscillations  are 
a  =  10°  - 10°  sin  a  t  and  a  =  10°  -  2°  sinco  t  .  In  addition,  pressure  distributions  are 
given  for  unsteady  reduced  frequencies  (k  =  27tfc/U00)  of  0.1,  0.2,  0.5,  and  1.0. 

Secondly,  the  unsteady,  compressible,  and  viscous  N-S  code  is  used  in 
conjunction  with  the  previously  described  flow  transition  and  turbulence  models 
to  predict  the  surface  and  flow  field  behavior  during  the  forced  oscillations. 
Specifically,  the  code  is  used  to  predict  flow  behavior  in  the  separation  bubble 
region  near  the  airfoil  leading  edge,  and  the  subsequent  vortex-shedding 
phenomena  for  the  low-amplitude  and  low  reduced  frequency  test  case. 

1.         UPOT  Pressure  Distributions  ( a  =  10°  - 10°  sin  co  t ) 

Surface  pressure  coefficients  predicted  by  the  UPOT  code  are  shown  for 
a  =  10°  -  2° sinco  t  and  k  =  0.1  in  Figure  6.23.  Pressure  distributions  are  given  for 
airfoil  motion  in  opposing  directions  (up/down),  but  at  the  same  angle  of  attack. 
No  results  are  given  for  angles  of  attack  greater  than  12°,  since  the  unsteady  stall 
was  known  a  priori  to  be  only  slightly  higher  than  12°.  The  suction  pressure 
peaks  move  slightly  forward  with  increasing  angle  of  attack  as  was  the  case  for 
the  steady  angles  of  attack,  but  pressure  lag  effects  attributable  to  the  high- 
amplitude  forced  oscillation  are  very  much  evident  even  at  the  low  reduced 
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frequency  (k  =  0.1).  During  the  upstroke,  the  suction  pressure  peaks  are 
significantly  smaller  than  the  steady  suction  pressures  at  8,  10,  and  12°  angle  of 
attack,  as  shown  in  Figures  6.1  and  6.2.  During  the  downstroke,  on  the  other 
hand,  the  suction  pressure  peaks  are  significantly  larger.  This  lag  effect  in  the 
build  up  of  the  pressure  attributable  to  the  pitch  oscillation  is  seen  to  be  largest  at 
angles  of  attack  near  10°  where  the  pitch  rate  reaches  its  maximum  value. 
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Figure  6.23.    Unsteady  Surface  Pressure  Distributions  Near  the  Airfoil  Leading 

Edge  (UPOT),  a  =10°  -10°  sin  cot  =  8°,  10°  and  12°  (up /down),    k  =  0.1,  pitch  axis 
at  quarter  chord  point. 


The  effect  of  increased  oscillation  rates  or  reduced  frequencies  on  the 
pressure  distribution  during  airfoil  upward  motion  are  shown  in  Figure  6.24.  The 
example  is  for  12°  angle  of  attack  (up),  indicating  smaller  suction  peaks  with 
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increased  oscillation  rates.  However,  the  suction  peak  lag  effect  is  seen  to 
decrease  with  increasing  frequency.  At  k  =  1.0,  the  pressure  peak  is  3%  less  than 
at  a  reduced  frequency  of  0.5  and  only  slight  pressure  reductions  are  predicted  in 
the  first  2%  chord  downstream.  The  most  dramatic  pressure  peak  reduction 
relative  to  the  steady  state  pressure  value  occurs  for  the  lowest  reduced 
frequency  of  k  =  0.1  near  the  leading  edge  with  the  trend  continuing  to  10% 
chord  downstream. 
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Figure  6.24.  Steady  and  Unsteady  Surface  Pressure  Distributions  Near  the  Airfoil 
Leading  Edge  (UPOT),  a  =  10°-  10° sincot  =  12°(up),  k  =  0.1,  0.2,  0.5  and  1.0, 
pitch  axis  at  quarter  chord  point. 
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2.         UPOT  Pressure  Distributions  ( a  =  10°  -  2°  sin  co  t ) 

Surface  pressure  coefficients  predicted  by  the  UPOT  code  for  the  low- 
amplitude  test  case  (a  =  10°  -2°sinco  t)  are  shown  in  Figure  6.25a,  b,  and  c,  each 
with  increasing  airfoil  pitch  rates  or  reduced  frequencies  of  k  =0.1,  0.2,  and  0.5. 
Results  are  presented  for  airfoil  motion  in  opposing  directions  (up/down)  at  10° 
angle  of  attack  and  at  the  top  (a  =  12°)  and  bottom  (a  =  8°)  of  the  oscillation 
cycle.  As  pitch  rates  increase  from  k  =  0.1  to  0.5,  suction  pressure  peaks  decrease 
slightly  at  12°  angle  of  attack  and  increase  slightly  at  8°  angle  of  attack. 
However,  the  pressure  lag  effect  is  still  clearly  visible  at  the  three  reduced 
frequencies  in  the  pressure  differences  between  the  upstroke  and  downstroke  at 
10°  angle  of  attack. 
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(a)a  =  10°-2°sincot  =  8°  and  10°(up/ down),  and  12° (top),  k  =  0.1. 

Figure  6.25.    Unsteady  Surface  Pressure  Distributions  Near  the  Airfoil  Leading 
Edge  (UPOT),  pitch  axis  at  quarter  chord  point. 
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(b)  a  =  10° -2°  sin  cot  =  8°  and  10°  (up  /down),  and  12°  (top),  k  =  0.2. 
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(c)  a  =  10°  -2°  sin  cot  =  8°  and  10°  (up /down),  and  12°  (top),  k  =  0.5. 


Figure  6.25.  (Contd.) 
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The  magnitude  of  the  lag  effect  is  quantified  even  better  in  Figure  6.26a,  b, 
c,  and  d  where  the  pressures  induced  by  the  pitch  oscillations  at  k  =  0.1,  0.2,  0.5, 
and  1.0  are  compared  with  the  steady  pressures  at  8°,  12°,  and  10°  angle  of 
attack  during  the  up-  and  downstroke. 
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Figure  6.26.  Steady  and  Unsteady  Surface  Pressure  Distributions  Near  the  Airfoil 
Leading  Edge  (UPOT),  pitch  axis  at  quarter  chord  point. 
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(c)a  =  10°-2°sincot  =  10°(up),  k  =  0.1,  0.2,  0.5  and  1.0. 
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(d)a  =  10°-2osinG)t  =  10°(down),  k  =  0.1,  0.2,  0.5  and  1.0. 


Figure  6.26.    (Contd.) 
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3.         N-S  Pressure  Distributions  ( a  =  10°  -  2°  sin  co  t ) 

Pressure  distributions  are  computed  with  the  N-S  code  for  the  low- 
amplitude  oscillation  (a  =  10°  - 2° sinco  t  )  at  a  reduced  frequency  ( k  =  0.1).  As  in 
the  steady  flow  computations,  the  Chen-Thyson  transition  model  is  used  and  the 
transition  onset  locations  were  specified  during  the  up-  and  downstroke,  as 
shown  in  Figure  6.27. 
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Figure  6.27.    Transition  Onset  Locations  for  Steady  and  Oscillatory  Airfoil 
Computations.    Moo=0.3,  cc  =  10°  -2°sincot,  k  =  0.0  and  0.1,   Re  =  540,000. 

A  comparison  of  the  N-S  pressure  distributions  in  Figure  6.28a  and  b  with 
the  potential  flow  computations  in  Figure  6.25a  show  that  the  inclusion  of 
viscous  and  compressible  flow  effects  produce  significantly  smaller  suction  peaks, 
ranging  from  -3.5  to  -6.3  for  the  N-S  results  versus  -4.3  to  -8.25  for  the  UPOT 
results  as  the  angle  is  changed  from  8°  to  12°.  Similar  to  the  steady  flow  results, 
the  formation  of  a  leading-edge  separation  bubble  is  clearly  visible  throughout 
most  of  the  oscillation  cycle,  changing  from  a  very  small  bubble  at  the  bottom  of 
the  cycle  at  a  =  8°  to  a  bubble  that  extends  from  about  1  %  of  chord  to  4%  of 
chord  at  a  =  10°  and  then  becomes  somewhat  smaller  at  the  top  of  the  cycle  at 
oc  =  12°. 
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Figure  6.28.  Unsteady  Surface  Pressure  Coefficients  Near  the  Airfoil  Leading 
Edge  (N-S  Code,  275x81).  M00=  0.3,  oc  =  10°  -2°  sin  cot,  k  =  0.1,  Re  =  540,000, 
Gy    =  200. 
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If  one  overlays  the  computed  results  at  a  =  10°  for  both  the  up-  and 
downstroke  as  done  in  Figure  6.29,  one  sees  that  both  the  viscous  and  inviscid 
pressure  distributions  show  the  pressure  lag  effect. 

In  Chapter  V,  reference  is  made  to  the  experimental  flow  visualization 
results,  which  clearly  show  the  shedding  of  a  vortex  during  the  airfoil 
downstroke  shortly  after  passing  the  top  of  the  oscillation  cycle.  It  is  important 
to  draw  attention  to  the  fact  that  the  N-S  computations  similar  to  the  UPOT 
calculations  do  not  predict  the  vortex-shedding  occurrence  evident  in  the 
measurements  on  the  airfoil  downstroke,  but  rather  indicate  a  gradual  build-down 
of  the  separation  bubble. 
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Figure  6.29.  Unsteady  Surface  Pressure  Coefficient  Near  the  Airfoil  Leading 
Edge  (UPOT  and  N-S  Codes).  ot=  10°-2°sincot  =  8°(bottom),  10°(up/down), 
12°(top),  k  =  0.1.  N-S  analysis,  MM=  0.3,  Re  =  540,000,  Gy    =200. 
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4.         N-S  Skin-Friction  Coefficients 

Skin-friction  coefficient  predictions  by  the  N-S  code  during  the  airfoil 
upstroke  are  given  in  Figure  6.30.  The  beginning  of  the  negative  Cf  region 
moves  forward  on  the  airfoil  surface  with  increasing  angle  of  attack  from  x/c  = 
0.025  at  a  =  8°  to  x/c  =  0.01  at  a  =  1 1.2°.  A  minimum  Cf  of  -0.0036  occurs  at  x/c 
=  0.038  and  a  =  8°,  whereas  at  a  =  1 1.2°,  Cf  reaches  a  minimum  of  -0.023  at  x/c  = 
0.028.  Positive  CfS  are  predicted  again  at  x/c  =  0.054  and  a  =  8°  and  at  x/c=  0.50 
for  a  =  10°,  which  are  downstream  of  the  trailing  edge  of  the  predicted  pressure 
plateaus.  Thus,  the  lateral  extent  of  negative  CfS  over  the  airfoil  surface  is  very 
similar  in  location  to  the  measured  separation  bubble  area. 
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Figure  6.30.     Skin-Friction  Coefficient  Distributions  (N-S  Code  -  275x81). 

Moo  =  0.3,  a  =  10°  -  2°sincot  =  8°,  8.8°(up),  10°  and  1 1.2° (up),  k  =  0.1,  Re  =  540,000, 
G7    =200. 
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At  a  =  1 1.2°,  CfS  oscillate  between  positive  and  negative  values.  In  effect, 
the  N-S  code  predicts  multiple  reverse  flow  regions  or  multiple  leading  edge 
vortices,  which  were  not  apparent  in  the  experiments.   In  addition,  the  N-S  code 
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does  not  predict  the  vortex-shedding  process  on  the  airfoil  downstroke  and,  thus, 


no  skin-friction  distributions  are  given. 


5.        N-S  Velocity  Profiles 


Velocity  profiles  are  shown  for  only  the  airfoil  upstroke  because  of  the 
inability  of  the  N-S  analysis  to  predict  the  vortex-shedding  and  light  dynamic  stall 
process  that  was  evident  in  the  measured  data.  Unsteady  velocity  profile 
predictions  from  the  airfoil  leading  edge  to  7%  chord  in  Figure  6.31a  are  similar  to 
the  steady  results  given  previously  at  8°  angle  of  attack  in  Figure  6.18.  This 
should  not  be  surprising  in  light  of  the  temporary  zero  pitch  rate  at  the  bottom  of 
the  oscillation  cycle.  The  flow  is  predicted  to  reach  slightly  higher  velocities 
(1.79  UM  versus  1.77  U^)  during  airfoil  upward  motion  at  the  same  forward  airfoil 
location  of  x/c  =  0.0123.  The  unsteady  airfoil  motion  lessens  reverse  flow 
predictions  slightly  (-0.04  U^  vs.  -0.09  U^)  at  the  same  downstream  (4.1%)  chord 
location.  Negative  velocities  are  shown  at  two  downstream  locations  (x/c  = 
0.0298  and  0.0413)  versus  three  surface  locations,  including  one  node 
downstream  for  the  steady  angle  of  attack  at  8°. 
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Figure  6.31.   Boundary-Layer  Velocity  Profiles  (N-S  Code,  275  x  81).  Moo  =  0.3, 
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a  =  10°-2°sincot,  k  =  0.1,  Re  =  540,000,  Gv     =200. 


Figure  6.31b  gives  velocity  profiles  for  8.8°  (up)  angle  of  attack. 
Comparing  results  with  the  predicted  flow  at  the  bottom  of  the  oscillation  cycle, 
the  extent  of  negative  velocities  has  increased  slightly.  Negative  velocities  are 
now  shown  one  node  forward  on  the  airfoil  surface  at  x/c  =  0.0201,  even  though 
the  specified  transition  onset  location  remained  at  x/c  =  0.025.  Therefore,  the  N-S 
analysis  predicts  a  laminar  flow  separation  at  8.8°  angle  of  attack.  The  maximum 
flow  velocity  (1.85 U^)  was  predicted  at  x/c  =  0.0123  and  y/c  =  0.001.  The 
maximum  negative  flow  velocity  (-0.06  UM)  was  predicted  at  x/c  =  0.0298  and 
y/c  =  0.00016. 
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Figure  6.31  (Contd.) 

The  velocity  profiles  for  a  =  10°  (up)  are  given  in  Figure  6.32a.  Transition 
onset  was  specified  at  x/c  =  0.025,  which  was  slightly  aft  of  the  location  used  in 
steady  flow  analysis.  Velocity  profiles  with  inflection  points  and  the  extent  of  the 
negative  velocity  region  are  less  than  those  predicted  for  the  steady  flow  at  10° 
angle  of  attack,  given  previously  in  Figure  6.19.    The  negative  velocity  region 
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extends  from  0.0201  <  x/c  <  0.051  along  the  airfoil  surface  and  to  a  lesser  height 
of  y/c  =  0.00079  above  the  airfoil  surface.  Slightly  higher  velocities  (1.97U.. 
versus  1.95  U^)  than  for  the  steady  analysis  occur  at  x/c  =  0.0201  and  y/c  = 
0.00158. 
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(a)  a  =  10°(up). 

Figure  6.32.   Boundary-Layer  Velocity  Profiles  (N-S  Code  -  275x  81).  M^  =0.3, 
a  =  10°-2°sincot,  k  =  0.1,  Re  =  540,000,  Gv    =200. 

Mr 

The  transition  onset  location  was  moved  forward  to  x/c  =  0.0177  at  1 1.2° 
(up)  angle  of  attack.  Figure  6.32b  shows  that  the  predicted  reverse  flow  region 
has  broadened  and  moved  forward  on  the  airfoil  surface  compared  with  the 
results  for  a  =  10°  in  Figure  6.32a.  At  this  angle  of  attack,  the  flow  reached  even 
higher  velocities  (2.13UJ  in  the  boundary  layer  at  2%  chord  downstream    and 
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y/c  =  0.0021.  A  maximum  negative  velocity  (U/  U„„  =-0.41)  is  predicted  at  x/c  = 
0.0298  and  y/c  =  0.00016.  The  negative  velocity  region  extends  from  0.01  <  x/c 
<  0.08  in  the  downstream  direction  and  from  0.0  <  y/c  <  0.00077  above  the  airfoil 
surface.  Because  of  the  low-frequency  oscillation,  a  reasonable  assumption  is 
that  the  N-S  velocity  profiles  predict  a  separation  bubble  from  1  to  8%  chord 
along  the  surface  and  0.2%  chord  in  height  above  the  airfoil  surface. 
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(b)  <x  =  11.2°(up). 

Figure  6.32.  (Contd.) 

6.         N-S  Lift,  Drag,  and  Moment 

Airfoil  section  lift,  drag,  and  pitching  moment  predicted  by  the  N-S  analysis 
during  the  oscillation  cycle  are  given  in  Figure  6.33  a,  b,  and  c.  Transition  onset 
was  specified  in  the  range  from  0.012  <  x/c  <  0.025  with  forward  movement  at 
the  higher  angles  of  attack.  Lift  and  pitching  moment  results  from  the  UPOT  code 
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are  included  in  Figure  6.33a  and  c.  It  is  seen  that  the  N-S  code  predicts  a 
clockwise  lift  curve  progression  with  more  lift  predicted  on  the  airfoil  upstroke. 
In  contrast,  the  potential  flow  code  (UPOT)  predicts  a  counterclockwise  lift 
hysteresis  loop.  However,  as  expected,  the  N-S  predicted  lift  levels  are  lower  than 
the  inviscid  lift  predictions,  ranging  from  11%  to  17%  less  at  the  top  of  the 
oscillation  cycle.  The  pitching  moment  predictions  show  counterclockwise  loops, 
using  either  the  N-S  or  UPOT  codes.  However,  the  viscous  analysis  exhibits 
positive  values  throughout  the  entire  cycle.  The  N-S  drag  hysteresis  loop  is 
clockwise  like  the  lift  loop  with  more  drag  predicted  during  the  airfoil  upstroke 
and  at  the  higher  angles  of  attack.  Since  the  focus  of  this  work  was  in  the 
separation  bubble  region  near  the  airfoil  leading  edge,  no  PDI  measurements  were 
made  for  the  entire  region  surrounding  the  airfoil  surface  to  make  comparisons 
with  the  predicted  loads.  However,  future  work  could  easily  provide  such 
comparisons. 
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(a)  Computed  Lift,  a  =  1 0°  -  2°  sin  o  t,  k  =  0. 1,  Re  =  540, 000 . 

Figure  6.33.    Unsteady  Computed  Lift,  Drag,  and  Pitching  Moment  Coefficients 
(Panel/Boundary  Layer  and  N-S  Codes).  N-S  analysis,  MM  =  0.3,  Gy    =  200 . 
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(b)  Computed  Drag,  a  =  10°  -  2°  sinco  t,  k  =  0.1,  Re  =  540,000. 
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(c)  Computed  Pitching  Moment,  a  =  10°  -  2°  sinco  t,  k  =  0.1,  Re  =  540,000. 

Figure  6.33.  (Contd.) 
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VH.      COMPARISONS  OF  MEASURED  AND  COMPUTED  RESULTS 

Laser-Doppler  velocimeter  (LDV)  velocity  data  and  Point  Diffraction 
Interferometry  (PDI)  density  contours  with  resultant  pressure  distributions  are 
used  to  validate  the  computational  methods.  The  flow  considered  for  these 
comparisons  is  (1)  compressible  and  in  a  transitional  Reynolds  number  range,  (2) 
predominated  by  a  leading-edge  separation  bubble  at  higher  angles  of  attack,  and 
(3)  undergoing  light  dynamic  stall  on  the  airfoil  downstroke  below  the  static  stall 
angle  of  attack  for  the  unsteady  analysis. 

The  PDI  and  LDV  measurements  both  showed  the  presence  of  a  separation 
bubble  and  its  subsequent  bursting  on  the  downstroke  in  light  dynamic  stall, 
pointing  to  the  need  for  transitional  flow  analysis.  The  laminar  separation  bubble 
near  the  leading  edge  alters  the  flow  characteristics  considerably,  as  dynamic  stall 
initiates  from  the  bursting  bubble  as  opposed  to  trailing-edge  flow  reversal. 

In  Chapter  VI,  the  potential  and  boundary-layer  analyses  were  found 
useful  and  comparable  to  the  N-S  results  at  the  lower  angles  of  attack.  However, 
beginning  at  a  =  8°  with  a  separation  bubble  present,  (1)  the  boundary- layer  code 
predicted  only  a  limited  amount  of  reverse  flow  and  no  solution  at  higher  angles 
of  attack,  and  (2)  the  potential  code  predicted  excessive  suction-pressure  peaks 
and,  of  course,  no  separation  bubble.  Thus,  the  summary  comparisons  are  made 
between  the  measurements  and  the  N-S  analysis,  which  accounts  for  the  viscous, 
transition,  turbulent,  and  compressible  flow  effects. 

For  the  light  dynamic  stall  case  being  considered,  viscous  effects  in  the 
direction  normal  to  the  airfoil  surface  are  on  the  order  of  the  airfoil  thickness. 
Therefore,  the  prediction  of  the  flow  characteristics  was  performed  by  the  thin- 
layer,  compressible,  Reynolds-averaged,  N-S  equation  formulation. 

The  unsteady  LDV  and  PDI  measurements  were  among  the  first  attempted 
and  successfully  made  for  the  oscillatory  cases  that  are  presented  in  this  study. 
The  low-amplitude  case  with  light  dynamic  stall  gives  an  excellent  data  set  for 
evaluation  of  the  capability  of  computational  schemes  to  predict  the  unsteady 
flow  phenomena.  The  dynamic  stall  process  is  very  complex,  requiring  a  more  in- 
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depth  study  of  the  underlying  flow  physics.  To  this  end,  the  analysis  is  the  fust  to 
critically  evaluate  the  effect  of  transition  onset  locations  and  transition  length  on 
the  flow  physics  as  determined  by  making  comparisons  with  experimental 
measurements  for  identical  flow  conditions. 

A.        INTRODUCTION 

Two  predominant  factors  were  found  to  be  critical  in  predicting  separation 
bubbles  that  were  similar  to  the  measured  ones.  The  two  factors  were  (1)  the 
transition  onset  location  and  (2)  the  flow  transition  length  over  the  airfoil  surface. 
As  stated  before,  the  Chen-Thyson  transition  length  model  [Ref.  20]  was  used  in 
the  present  analysis,  where  flow  intermittency  is  defined  to  change  from  0  to  1.0 
during  transition  from  a  fully  laminar  to  fully  turbulent  flow.  In  the  region  of 
intermittency,  turbulent  spots  appear  with  progression  in  the  streamwise 
direction.  For  pressure  distributions  and  skin-friction  predictions  that  are  in  better 
agreement  with  the  measurements,  the  transition  length  is  changed  by  varying  the 
value  (Gy  )  in  the  formulation.    The  transition  model  is  repeated  here  for 

completeness. 
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(7.1) 


Ytr(x)  =  l-exp 


(--iT)Rei34(x-xtr)£^ 


G7  i/'      Aff  "JX,rue 

Itr  c  / 


The  Chen-Thyson  model  had  already  been  used  successfully  in 
conjunction  with  a  viscous/inviscid  interaction  analysis  method  by  Walker  et  al. 
[Ref.  30]  to  predict  some  features  of  a  separation  bubble  on  a  low  Reynolds 
number  airfoil.  The  analysis  predicted  a  separation  bubble  over  a  NACA  65-213 
airfoil  from  60  to  84%  chord,  which  was  similar  to  the  experiments  of  Hoheisel  et 
al.  [Ref.  67]  when  the  value  was  reduced  to  40.  The  Chen-Thyson  transition 
model  was  developed  for  high  Reynolds  number  flows  and  specified  a  Gy    value 

of  1200.  The  analysis  of  the  NACA  65-213  airfoil  was  done  for  an  angle  of  attack 
of  0°  and  a  Reynolds  number  of  240,000. 

The  second  critical  factor  for  obtaining  numerical  results  comparable  with 
the  measurements  is  the  proper  location  for  transition  onset.  Michel's  criterion, 
which  was  developed  for  attached  flows  at  higher  Reynolds  numbers,  is  used 
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initially  to  estimate  the  transition  onset  location.  Ekaterinaris  et  al.  [Ref.  62]  used 
Michel's  criterion  and  Chen-Thyson's  model  with  some  success  for  a  flow  at  a 
higher  Reynolds  number  (4xl06)  with  deep  dynamic  stall  and  a  lightly  stalled, 
lower  Reynolds  number  flow  [Ref.  58].  Excessive  surface  pressures  were 
reported  for  potential  and  fully  turbulent  flow  analysis,  but  flow  transition 
consideration  led  to  results  more  in  line  with  the  measurements.  For  flows  at 
lower  Reynolds  numbers,  transition  onset  locations  were  expected  to  be  aft  of 
Michel's  prediction.  Thus,  transition  onset  locations  were  varied  to  give 
computed  results  that  agreed  best  with  the  measurements. 

The  premise  for  this  work  is  that  if  appropriate  transition  lengths  and 
transition  onset  locations  are  found,  then  the  physics  of  a  leading-edge  separation 
bubble  could  be  described  sufficiently  and  adequately.  In  this  work,  the 
sensitivity  of  the  transition  onset  location  and  transition  length  parameters  are 
tested  by  comparing  computed  results  with  the  experimental  measurements  at 
angles  of  attack  near  the  stall  with  substantial  pressure  gradients.  Pressure 
distributions,  skin-friction  coefficients,  and  velocity  profiles  over  the  upper 
leading  edge  of  the  airfoil  are  shown  to  change  substantially  when  slight  changes 
are  specified  in  location  of  transition  onset  and  transition  length. 

B.        STEADY  FLOW  RESULTS 

1.         Transition  Onset  and  Length  Effects 

As  an  example  of  the  effect  of  the  Gy     value  (in  Equation  7.1)  on  the 

transition  length,  results  are  presented  for  a  =  8°  with  transition  onset  specified  at 
x/c  =  0.04  in  Figure  7.1.  The  a  =  8°  case  was  chosen  because  it  was  the  starting 
point  for  the  low-amplitude  forced-oscillation  calculations  and  it  is  also  a 
moderate  angle  of  attack  for  the  steady  flow  analysis.  Values  were  varied  from 
460,  as  computed  by  Cebeci's  correlation  formula  (Equation  4.27),  to  200  and 
then  50.  The  intermittency  showed  extreme  sensitivity  to  changes  in  the  value, 
reducing  the  transition  length  from  12%  to  less  than  4%  chord  in  the  value  range. 
The  final  value  used  in  the  analysis  was  chosen  by  observing  resultant  pressure 
and  skin  friction  predictions  that  agreed  best  with  the  measurements,  as  shown  in 
Figures  7.2  and  7.3. 
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Figure  7.1.  Effect  of  (GYff)  Transition  Value  on  Turbulence  Intermittency 
Distribution.   M^  =  0.3,  a  =  8°,  Re  =  540,000,  (x  /  c)tr  =0.04. 

Pressure  distributions  and  skin-friction  coefficients  are  shown  in  Figure  7.2 
for  the  previously  stated  conditions.  As  can  be  seen  by  the  streamwise  extent  of 
negative  skin-friction  values,  an  order  of  magnitude  change  in  the  value  from  460 
to  50  greatly  reduces  the  reverse  flow  (-Cf)  region  along  the  airfoil  surface  from 
5.4%  to  3.5%  chord.  For  the  smaller  value  of  50,  the  -Cf  region  extends  from 
0.023  <  x/c  <  0.058,  whereas  for  460,  the  -Cf  range  is  from  0.020  <  x/c  <  0.074. 
Note  that  the  Cf  values  are  multiplied  by  100  in  the  plots  for  better  presentation. 
For  all  values  of  Gy  ,  a  laminar  separation  bubble  is  predicted.  In  each  case, 
negative  CfS  are  shown  forward  of  the  specified  transition  onset  location  at  x/c  = 
0.04.   Additionally,  higher  values  of  Gy    give  longer  streamwise  plateaus  in  the 

pressure  distribution  that  are  typical  of  flows  with  separation  bubbles.  However, 
values  greater  than  400  predicted  early  flow  reattachment,  secondary  suction- 
pressure  peaks,  and  skin-friction  coefficients  that  become  positive  in  a  small  range 
before  again  turning  negative.  By  contrast,  the  experimental  results  do  not  show 
what  is  in  effect  a  secondary  separation  bubble,  or  multiple  leading-edge  vortical 
structures.  Thus,  for  the  predictions,  the  Gy  value  of  200  most  closely  matched 
the  measurements. 
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Figure  7.2.  Effect  of  (Gy  )  Transition  Value  on  Pressure  and  Skin-Friction 
Distributions.   M^  =  0.3,  a  =  8°,  Re  =  540,000,  (x  /  c)tr  =0.04. 

The  effect  of  different  transition  onset  locations  on  the  flow  are  shown  in 

Figure  7.3,  using  a  Gy    value  of  200  for  the  same  flow  at  a  =  8°.   Choice  of  the 

forward  transition  onset  location  (x/c  =  0.03)  predicts  a  reverse  flow  region  of 
2.9%  chord  at  the  surface  and  a  very  slight  pressure  plateau,  whereas  a  slightly 
aft  transition  onset  location  (x/c  =  0.04)  predicts  a  reverse  flow  region  of  4.1% 
chord  and  a  slight  secondary  pressure  peak  at  the  end  of  the  plateau.  The 
calculations  at  both  transition  onset  locations  predict  a  laminar-flow  separation. 
On  the  other  hand,  Michel's  criterion  predicts  transition  onset  at  2.5%  chord  (not 
shown)  and  a  transitional-flow  separation  with  a  reverse-flow  region  from  0.027 
<  x/c  <  0.051.  The  predicted  separation  bubble  in  this  case  is  very  minute  when 
compared  with  the  experimental  results.   For  transition  onset  locations  that  were 
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specified  downstream  on  the  airfoil  surface  (x/c  >  0.05),  the  numerical  scheme 
gave  the  formation  of  multiple  vortical  structures  over  the  airfoil  leading-edge 
region.  In  effect,  laminar  flow  extended  too  far  downstream  such  that  large 
multiple  non-physical  reverse-flow  regions  were  predicted.  From  the  combined 
results  shown  in  Figures  7.2  and  7.3,  transition  onset  was  specified  at  x/c  =  0.035 
and  Gy    was  set  equal  to  200  for  a  =  8°. 


U 

x 

o 

o 

a 

U 


4.0 

3.2 

2.4 

1.6" 

0.8" 

o.o- 

-0.8 


-1.6-H 


0.00 


N-S,  8deg  AOA,  Gconst=200 

Cp,  (x/c)tr  onset  =  0.04 

Cp,  (x/c)tr  onset  =  0.03 

Cf,  (x/c)tr  onset  =  0.04 

Cf,  (x7c)tr  onset  =  0.03 


0.02        0.04        0.06        0.08        0.10 
Chordwise  Location,  x/c 


0.12 


Figure  7.3.    Effect  of  Transition  Onset  Location  on  Pressure  and  Skin-Friction 

Distributions.  M^  =  0.3,  a  =  8°,  Re  =  540,000,  Gy    =  200. 
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2.         Density  Contour  Comparisons 

The  next  series  of  plots  overlie  computed  and  measured  PDI  density 
contours.  The  computed  density  ratio  is  ranged  to  give  the  same  number  of 
contours  in  the  flow  field  and  the  same  delta  (p/p0  =  0.0085)  for  each  dark  fringe 

to  provide  direct  comparisons  with  the  PDI  measurements.  Density  ratios  are 
noted  at  the  intersections  of  selected  equivalent  measured  and  computed 
contours.  The  comparisons  are  shown  in  Figure  7.4  at  a  =  0°  with  the  stagnation 
point  at  the  leading  edge  and  a  symmetry  of  contours  about  the  airfoil.    Eight 
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contours  are  shown  in  the  flow  field,  with  the  lowest  density  ratio  of  0.944 
corresponding  to  a  Mach  number  of  0.34  both  over  and  under  the  airfoil  leading 
edge  surface.  Even  though  complete  agreement  is  not  seen,  all  contours  are 
within  the  one  fringe  uncertainty  of  the  techniques. 
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Figure  7.4.    Computed  and  Measured  Density  Contours  (PDI)  Near  the  Airfoil 
Leading  Edge.  MM  =  0.3,  a  =  0°,  Re  =  540,000. 

At  a  =  6°,  Figure  7.5  shows  that  the  stagnation  point  has  moved  under  the 
leading  edge  of  the  airfoil  and  the  number  of  density  contours  has  doubled  to  16 
from  the  eight  contours  at  a  =  0° .  The  number  of  contours  over  the  upper  airfoil 
leading-edge  surface  has  increased  dramatically.  The  density  contour  previously 
lying  on  the  upper  airfoil  surface  at  a  =  0°  now  extends  vertically  from  the 
leading  edge  of  the  airfoil.  The  minimum  flow  field  density  ratio  of  0.876 
corresponds  to  a  local  Mach  number  of  0.52  over  the  airfoil  leading  edge  surface. 
No  separation  bubble  is  indicated  by  either  the  measurements  or  computations,  as 
all  contours  have  a  continuous  and  gradual  curvature  with  the  upper  airfoil 
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surface.  Computed  contours  converged  forward  of  the  measurements  on  the 
upper  leading  edge  of  the  airfoil  surface.  Part  of  the  reason  for  the  differences 
was  reported  by  Cho  et  al.  [Ref.  68],  resulting  from  the  large  density  gradients  at 
the  airfoil  leading-edge  and  optical  distortions  causing  a  slight  downstream  shift 
in  the  measured  separation  bubble. 
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Figure  7.5.    Computed  and  Measured  Density  Contours  (PDI)  Near  the  Airfoil 
Leading  Edge.  MM  =  0.3,  a  =  6°,  Re  =  540,000 . 

Figure  7.6  shows  a  comparison  of  density  contours  near  the  leading  edge 
of  the  airfoil  at  a  =  10°.  A  total  of  27  density-flow-field  contour  levels 
corresponds  to  a  highly  accelerated  flow  over  the  leading  edge  to  2.39  times  the 
free-stream  Mach  of  0.3  (M=  0.72).  The  formation  of  a  separation  bubble  is 
clearly  seen  in  the  region  where  the  density  contours  flatten  (plateau)  over  the 
airfoil  leading  edge  and  then  turn  perpendicular  towards  the  airfoil  surface  as  was 
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described  when  discussing  the  PDI  images.  The  flat  plateaus  in  the  measured 
contours  extend  to  greater  heights  above  the  airfoil  and  extend  further 
downstream  over  the  airfoil  surface  than  the  predicted  separation  bubble.  The 
measured  separation  bubble  is  estimated  to  extend  from  0.01  <  x/c  <  0.06  along 
the  surface  of  the  airfoil  by  including  all  contours  that  turn  perpendicular  to  the 
airfoil  surface  near  5%  chord.  By  using  the  same  criteria,  the  height  of  the 
separation  bubble  is  estimated  to  be  y/c  =  0.02.  The  overall  lateral  extent  of  both 
measured  and  computed  separation  bubbles  is  similar,  remaining  within  the  one 
fringe  uncertainty.  However,  the  predicted  separation  bubble  is  again  slightly 
forward  on  the  airfoil  surface. 
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Figure  7.6.    Computed  and  Measured  Density  Contours  (PDI)  Near  the  Airfoil 
Leading  Edge.  M^  =  0.3,  a  =  10°,  Re  =  540,000 . 

As  the  angle  of  attack  increases  from  10°  to  11.9°,  the  density  contour 
picture  does  not  change  significantly,  as  shown  in  Figure  7.7.  However,  the 
boundary  layer  thickens  slightly  aft  of  the  separation  bubble  and  both  measured 
and  computed  separation  bubbles  extend  to  greater  heights  above  the  airfoil 


143 


surface.  The  PDI  measurements  now  show  six  contours  as  opposed  to  five 
turning  perpendicular  to  the  airfoil  surface  at  their  downstream  edges.  Now  a 
two-fringe  disparity  exists  between  the  measurements  and  computations  at  some 
locations  in  the  flow  field. 
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Figure  7.7.    Computed  and  Measured  Density  Contours  (PDI)  Near  the  Airfoil 
Leading  Edge.  MM  =  0.3,  a  =  1 1.9°,  Re  =  540,000. 

3.         Pressure  Distribution  and  Skin-Friction  Comparisons 

Figure  7.8  shows  a  comparison  between  PDI  derived  and  computed 
pressure  coefficient  distributions  at  a  =  6°  in  steady  flow.  For  this  angle  of 
attack,  as  previously  stated,  no  separation  bubble  was  found  in  the  measurements. 
Likewise,  the  computed  skin  friction  remains  positive  and  thus,  predicts  no 
separation  bubble.  The  pressure  distributions  agree  very  well  with  each  other. 
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Differences  are  well  within  the  experimental  and  computational  uncertainties  with 
transition  onset  specified  at  x/c  =  0.04  and  a  transition  length  value  of  200. 
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Figure  7.8.    Computed  and  PDI  Derived  Pressure  Distributions  Near  the  Airfoil 
Leading  Edge.  M^  =  0.3,  a  =  6°,  Re  =  540,000,  (x  /  c)tr  =  0.04,  GYtr  =  200. 


In  Figure  7.9,  a  steady-flow  solution  is  given  for  a  =  8°  with  a  specified 
transition  onset  location  of  x/c  =  0.035.  The  measured  suction-pressure  peak  is 
slightly  lower  than  the  computed  value  and  the  extent  of  the  predicted  pressure 
plateau  is  less  than  the  measurements.  However,  the  end  of  the  separation  bubble 
cannot  be  precisely  determined  from  the  pressure  distributions  alone,  as 
previously  reported  by  Gaster  [Ref.  15]  and  Tani  [Ref.  5].  The  start  of  predicted 
flow  separation  as  determined  by  the  region  of  negative  skin-friction  coefficient 
agrees  well  with  the  start  of  the  measured  pressure  plateau.  Negative  CfS  are 
indicated  at  x/c  =  0.028,  which  is  forward  of  transition  onset,  and  therefore,  a 
laminar  flow  separation  is  predicted.    In  addition,  CfS  become  positive  at  5.9% 
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chord,  indicating  the  downstream  end  of  the  reverse  flow  region  that  nearly 
coincides  with  the  measured  pressure  recovery  location  at  6%  chord,  confirming 
agreement  between  the  two  data  sets. 
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Figure  7.9.    Computed  and  PDI  Derived  Pressure  Distributions  Near  the  Airfoil 
Leading  Edge.  M^  =  0.3,  a  =  8°,  Re  =  540,000,  (x  /  c)tr  =  0.035,  GYtr  =  200. 

Similar  results  can  be  seen  in  Figure  7.10  for  a  =  10°  with  transition  onset 
set  at  x/c  =  0.023.  Here,  the  computed  pressure  distributions  (Cp)  show  a  second 
rise  and  a  subsequent  fall  (secondary-separation  bubble),  and  the  extent  of  the 
pressure  plateau  again  is  smaller  than  that  seen  in  the  measurements.  The 
measured  suction-pressure  peak  is  downstream  of  the  predictions,  partially 
because  of  the  optical  distortions  stated  earlier,  however,  the  pressure  values  are 
within  the  one-fringe  uncertainty.  The  extent  of  the  reverse-flow  region 
indicated  by  negative  skin-friction  coefficients  gives  good  agreement  for  the 
separation-bubble  length,  and  again,  a  laminar  separation  is  predicted  by  the 
computations. 
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Figure  7.10.  Computed  and  PDI  Derived  Pressure  Distributions  Near  the  Airfoil 
Leading  Edge.  M^  =  0.3,  a  =  10°,  Re  =  540,000,  (x  / c)tr  =  0.023,  GYtf  =  200. 


Figure  7.11  gives  an  interesting  comparison  of  computed  pressures  with 
the  stalled  airfoil  PDI  results  at  a  =  12°  by  specifying  transition  onset  between 
0.012  <  x/c  <  0.017.  The  experiment  showed  that  a  =  12°  was  the  borderline  stall 
angle  of  attack  for  a  0.3  Mach  number  flow.  The  unsteady  nature  of  flow 
separation  at  stall  onset  gave  flows  that  were  intermittently  stalled  or  attached 
with  large  levels  of  suction  pressures.  For  the  case  of  transition  onset  at  x/c  = 
0.012  and  a  transition  length  value  of  200,  the  computations  reproduced  the 
stalled  experimental  results.  Slight  changes  in  the  transition  onset  location  and  a 
lower  transition  length  value  of  100  gave  a  completely  attached  flow  with  high 
levels  of  suction.  Thus,  the  numerical  modeling  with  a  simple  transition  model 
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appears  to  adequately  represent  the  physical  flow  in  the  wind  tunnel  for  steady 
angle  of  attack  analysis.  However,  the  oscillations  in  the  pressure  distribution 
near  the  leading  edge  are  incorrect  predictions  of  multiple  vortical  structures. 
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Figure  7.11.   Computed  and  PDI  Derived  Pressure  Distributions  Near  the  Airfoil 
Leading  Edge.  Mx  =  0.3,  a  =  12°,  Re  =  540,000,  (x  /  c)tr  =0.012,  GY(r  =  200. 

As  steady  angles  of  attack  are  increased  from  6  to  12°  in  the  previous 
solutions,  the  streamwise  transition  length  is  greatly  reduced.  Figure  7.12  shows 
that  at  a  =  6°  the  transitional  flow  intermittency  range  is  predicted  to  extend  from 
4%  to  12%  chord,  but  at  a  =  12°  the  intermittency  range  extends  only  from  1.7% 
to  4%  chord.  For  calculations  giving  the  best  agreement  with  the  measurements, 
Figure  7.12  also  shows  the  upstream  progression  of  the  specified  transition  onset 
locations  with  increasing  angles  of  attack.  In  reviewing  Equation  7.1,  the 
predominant  term  causing  a  reduced  transition  length  at  higher  angles  of  attack  is 
the  larger  boundary-layer  edge  velocity.    The  other  factors  in  the  equation  are 
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either  constant  or  shrinking  in  size  with  increasing  angles  of  incidence,  which 
would  extend  the  transition  length. 
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Figure  7.12  Angle-of- Attack  Effect  on  Intermittency,  Transition  Onset  Specified. 
MM  =  0.3,  Re  =  540,000,  Gv    =  200. 
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4.         Velocity  Profile  Comparisons 

Computed  and  measured  U-component  velocities  are  compared  for  angles 
of  attack  of  8,  10,  and  12°.  In  each  case,  velocity  profile  comparisons  begin  at  x/c 
=  -0.167,  which  is  0.5  inch  forward  of  the  airfoil  leading  edge,  and  then  proceed 
over  the  leading  edge  of  the  airfoil  from  0.0  <  x/c  <  0.167  in  0.033  chord  length 
intervals.  In  the  plots,  both  the  measured  and  computed  U  velocities  are 
normalized  by  the  velocity  (U0)  at  y/c  =  0.167.    An  interpolation  scheme  was 

used  to  predict  computed  velocities  at  the  nodes  of  the  measurement  grid.  Only 
the  positive  computed  velocities  are  shown  close  to  the  airfoil  surface,  whereas 
the  measured  velocities  are  restricted  to  some  distance  above  the  airfoil  because 
of  laser-beam  blockage.     For  all  angles  of  attack,  measured  and  computed 
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velocities  forward  of  the  airfoil  at  x/c  =  -0.167  are  in  full  agreement  and  give  the 
same  trend  with  slight  velocity  reductions  at  lower  hefghts. 

Normalized  velocities  at  8°  angle-of-attack  range  from  0.0  to  1.42,  as 
shown  in  Figure  7.13.  The  maximum  computed  velocity  in  the  flow  field  is  at  x/c 
=  0.033  near  the  airfoil  surface.  The  maximum  measured  velocity  is  aft  at  x/c  = 
0.067,  most  likely  because  of  the  measurement  height  restriction.  All  computed 
and  measured  profiles  for  x/c  >  0.067  are  essentially  in  full  agreement.  Measured 
values  at  y/c  <  0.0117  above  the  leading  edge  of  the  airfoil  are  less  than  the 
computed  values,  but  are  within  the  accuracy  of  the  LDV  system. 
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Figure  7.13.  Computed  and  Measured  Velocity  Profiles  Near  the  Airfoil  Leading 
Edge.  Moo=0.3,  a  =  8°,  Re  =  540,000. 

Comparisons  for  10°  angle  of  attack  are  shown  in  Figure  7.14. 
Maximum  velocities  have  increased  slightly  to  1.54U0,  occurring  at  locations 

similar  to  those  at  8  °  angle  of  attack.    Agreement  between  the  data  sets  is  even 
better  than  at  the  lower  angle  of  attack,  but  the  measurements  in  this  case  were 


150 


restricted  to  a  higher  location  (y/c  =  0.100)  above  the  airfoil  surface.  Maximum 
computed  and  measured  differences  are  <2%  over  the  airfoil  surface.  Again, 
slight  differences  are  seen  at  the  airfoil  leading  edge  and  the  lowest  measurement 
height,  where  measurements  tended  to  be  lower  than  the  computations. 
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Figure  7.14.  Computed  and  Measured  Velocity  Profiles  Near  the  Airfoil  Leading 
Edge.  Moo=0.3,  a  =  10°,  Re  =  540,000. 

For  11.95°  angle  of  attack  in  Figure  7.15,  the  data  comparison  region 
above  the  airfoil  surface  is  essentially  the  same  as  for  the  previous  case  at  8° 
angle  of  attack.  Velocities  over  the  airfoil  have  increased  slightly  to  1 .62  U0 ,  and 

agreement  between  the  data  sets  is  excellent  (<1%).  Again,  at  the  lowest 
measurement  height  near  the  airfoil  leading  edge,  slight  differences  are  seen. 
Because  of  LDV  measurement  system  limitations,  the  PDI  method  was  also  used 
to  make  a  better  assessment  of  the  extent  of  the  leading-edge  separation  bubble. 
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Figure  7.15.  Computed  and  Measured  Velocity  Profiles  Near  the  Airfoil  Leading 
Edge.  M^  =0.3,  a  =  11.95°,  Re  =  540,000. 

C.        OSCILLATING  AIRFOIL  FLOW  (a  =  10°  -  2°sincot,  k  =  0.1) 

The  PDI  image  in  Figure  5.20d  at  11.5°  angle  of  attack  (down)  verified  a 
light  dynamic  stall  and  the  vorticity  shedding  process  over  the  airfoil  surface 
during  the  downstroke  of  the  low-amplitude  oscillating  airfoil  flow.  Rapid  flow 
reattachment  occurs  as  evidenced  by  the  flow  picture  in  Figure  5.20e  at  10.5° 
angle  of  attack  on  airfoil  downward  motion.  Only  the  N-S  analysis  method 
adequately  modeled  separation  bubble  formation  and  development  (the  focus  of 
this  work)  at  the  higher  angles  of  attack  over  the  airfoil  leading  edge.  Thus,  only 
N-S  and  experimental  data  comparisons  and  summary  analysis  are  given  for  the 
low-amplitude  oscillations  (MM  =  0.3,  a  =  10°  -  2°  sin  cot,  k  =  0.1,  Re  =  540,000).  In 
addition,  oscillating  airfoil  results  are  given  only  for  airfoil  upward  motion  at 
angles  of  attack  less  than  11.5°  because  of  the  limited  number  of  validated  LDV 
samples  at  higher  angles  of  incidence  and  an  inability  of  the  analysis  to  predict 
light  dynamic  stall  during  the  airfoil  downstroke. 
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1.         Density  Contour  Comparisons 

As  a  result  of  optical  path  length  differences,  the  dark  and  light  fringes  that 
were  recorded  on  film  represent  density  contours  in  the  flowfield.  Again,  the 
computed  density  ratio  is  ranged  to  give  the  same  number  of  contours  in  the  flow 
field  and  the  same  delta  for  each  dark  fringe  (0.0085)  as  for  the  PDI  measurements 
to  make  direct  data  comparisons.  Thus,  the  figures  give  a  direct  comparison  of 
computed  and  measured  fringe  locations  in  the  vicinity  of  the  separation  bubble. 
Density  ratios  are  noted  in  the  flow  field  at  intersections  of  selected  equivalent 
measured  and  computed  contours. 

Computed  and  measured  density  contours  are  overlaid  for  a  =  8°  angle  of 
attack  (bottom  of  the  oscillatory  cycle)  in  Figure  7.16.  There  are  23  density 
contour  levels  shown  (one  less  than  for  a  =  8°  steady)  with  a  peak  Mach  number 
of  0.65  over  the  upper  surface  of  the  airfoil.  The  predicted  or  computed  contours 
in  the  separation  bubble  region  are  shifted  forward  of  the  measurements  on  the 
airfoil  surface.  The  trailing  edge  of  the  measured  contours  turn  abruptly  toward 
the  airfoil  surface  between  6%  and  8%  chord  downstream,  converging  with  the 
airfoil  surface  0.7%  to  1.7%  chord  downstream  of  equivalent  computed  contours. 
The  computed  contours  converge  more  gradually  with  the  airfoil  surface  at  the 
trailing  ends,  predicting  a  less  defined  separation  bubble.  Computed  contours 
begin  0.9%  to  1.3%  chord  forward  of  the  measured  contours  in  an  estimated 
separation  bubble  region  from  1%  to  8%  chord  downstream.  Measured  density 
contours  extend  0.8%  chord  higher  above  the  airfoil  surface  than  the  predicted 
contours. 

Computed  and  measured  density  contours  are  presented  for  a  =  8.8°  (up) 
in  Figure  7.17.  Contour  levels  (25)  are  shown  with  a  minimum  density  ratio  in  the 
flow  field  of  0.799,  which  corresponds  to  a  peak  local  Mach  number  of  0.68  over 
the  airfoil  surface.  The  measured  contours  parallel  the  airfoil  surface  for  some 
distance  downstream  of  the  suction  peak  and  then  turn  even  more  abruptly 
toward  the  airfoil  surface  at  the  trailing  ends  between  6%  and  8%  chord 
downstream.   In  the  separation  bubble  region,  computed  contours  begin  0.7%  to 
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Figure  7.16.   Computed  and  Measured  Density  Contours  (PDI)  Near  the  Airfoil 
Leading  Edge.     Moo=0.3,  a  =  10°-  20sincot  =  80,  k  =0.1,  Re  =  540,000,  (x/c)tr  = 
0.025,  Gv    =200. 

'tr 

1.0%  chord  forward  of  equivalent  measured  contours,  and  the  trailing  ends 
converge  with  the  airfoil  surface  as  much  as  2%  chord  forward  of  the  measured 
contours.  Also,  measured  contours  are  as  much  as  2%  chord  greater  in  height 
above  the  airfoil  surface  than  equivalent  predicted  contours. 

Figure  7.18  compares  density  contours  at  a  =  10°  (up).  A  total  of  24 
contours  indicates  that  the  flow  reaches  2.23  times  (M=  0.67)  the  free-stream 
Mach  number.  The  formation  of  a  separation  bubble  is  shown  in  the  region  along 
the  airfoil  surface  between  0.01  <  x/c  <  0.07,  where  both  measured  and  computed 
density  contours  flatten  over  the  airfoil  surface  and  then  turn  almost 
perpendicular  to  the  airfoil  surface  at  the  downstream  ends.  At  this  angle  of 
attack,  the  flat  plateaus  in  the  measured  contours  extend  to  only  slightly  greater 
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Figure  7.17.   Computed  and  Measured  Density  Contours  (PDI)  Near  the  Airfoil 

Leading    Edge.     M^  =  0.3,  a  =  10° -2° sin  cot  =  8.8°  (up),  k  =  0.1,  Re  =  540, 000, 
(x/c)tr  =0.025,  GYr  =  200. 

heights  (0.5%  chord)  above  the  airfoil.  But,  major  differences  are  evident  at  the 
trailing  ends  of  the  density  contours,  where  measured  contours  are  up  to  2.5% 
chord  downstream  of  the  equivalent  predicted  contours.  Again,  computed 
contours  at  the  forward  edge  of  the  separation  bubble  begin  on  the  order  of  1 .0% 
chord  forward  of  equivalent  measured  contours.  The  numerical  predictions  for 
the  oscillating  airfoil  at  a  =  10°  have  fewer  density  contours  in  the  flow  field  than 
in  the  corresponding  flow  for  steady  analysis  at  the  same  angle  of  attack.  This  is 
attributable  to  the  lag  effects  that  are  inherent  to  oscillatory/unsteady  flows. 
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Figure  7.18.   Computed  and  Measured  Density  Contours  (PDI)  Near  the  Airfoil 
Leading  Edge.  MM  =  0.3,  a  =  10°  -  2°sincot  =  10°(up),  k  =  0.1, Re=  540,000,  (x/c)tr  = 
0.025,  Gy    =200. 
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As  the  airfoil  pitches  up  from  oc=  10°  to  11.2°,  the  number  of  density 
contours  in  the  flow  field  has  increased  to  29  and  the  local  flow  attains  a 
maximum  Mach  number  of  0.75.  Figure  7.19  shows  that  measured  and  predicted 
contours  are  more  densely  packed  over  the  leading-edge  airfoil  surface.  Some 
inner  measured  contours  are  lower  in  height  over  the  airfoil  surface  than  the 
predictions  and  are  only  slightly  higher  than  the  predictions  in  the  flow  field  at  a 
density  ratio  of  0.833,  for  example.  Both  measured  and  computed  density 
contours  downstream  of  the  suction  peak  extend  along  the  airfoil  surface  much 
the  same  as  those  at  10°  angle  of  attack.  Again,  major  differences  are  seen  near 
the  downstream  end  of  the  separation  bubble,  where  the  measured  contours 
converge  abruptly  with  the  airfoil  surface  and  are  as  much  as  3.0%  chord 
downstream  of  equivalent  predicted  contours,  which  merge  more  gradually  with 
the  airfoil  surface.  The  picture  at  the  leading  edge  of  the  separation  bubble 
remains,  however,  the  same.  Recall,  as  reported  by  Cho  et  al.  [Ref.  58],  the  aft 
shift  in  the  measured  suction  pressure  peak  area  is  in  part  attributable  to  optical 
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distortions  as  a  result  of  the  large  density  gradients  over  the  airfoil  leading  edge. 
Overall,  the  density  comparisons  give  a  measured  separation  bubble  that  is  up  to 
2%  chord  greater  in  length  and  height  than  the  predicted  separation  bubble. 
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Figure  7.19.   Computed  and  Measured  Density  Contours  (PDI)  Near  the  Airfoil 
Leading  Edge.    M^  =  0.3,    a=  10°  -2°sincot  =  11.2° (up),    k  =  0.1,    Re  =540,000, 
(x/c)tr  =0.017,  GYtr=200. 

2.         Pressure  Distributions  and  Computed  Skin-Friction  Coefficients 

Measured  and  computed  pressure  distributions  at  the  extremes  of  the 
oscillation  cycle,  a  =  8°  and  12°,  are  compared  in  Figure  7.20.  There  is  good 
agreement  between  a  majority  of  the  predicted  and  measured  pressure 
distributions  at  a  =  8°  angle  of  attack  with  transition  onset  specified  at  x/c  = 
0.025.  However,  the  measured  suction-pressure  peak  is  17%  lower  than  the 
computed  value  and  the  extent  of  the  pressure  plateau  is  greater  than  in  the 
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computations.  Recall,  as  stated  before,  that  Gaster  [Ref.  15]  and  Tani  [Ref.  5] 
reported  that  the  end  of  the  separation  bubble  cannot  be  precisely  determined 
from  the  pressure  distributions  alone.  The  computed  suction  pressure  peak  is 
slightly  greater  than  that  previously  given  for  8  °  steady  angle  of  attack  in  Figure 
7.10  attributable  to  airfoil  oscillation  in  the  range  between  8  and  12°  angle  of 
attack.  The  computed  skin  friction  distribution  gives  negative  values  in  the 
separation  bubble  region  from  2.3%  to  5.8%  chord  downstream. 

At  a  =  12°  with  transition  onset  specified  at  x/c  =  0.012,  on  the  other  hand, 
the  computed  suction  pressures  significantly  overshoot  the  measured  ones.  But, 
the  overall  negative  computed  skin  friction  area  (0.0116  <  x/c  <  0.08)  still 
indicates  a  separation  bubble  of  comparable  size  to  the  measured  one  along  the 
airfoil  surface.  As  expected,  the  onset  of  stall  is  delayed  during  the  upstroke 
yielding  a  greater  suction  pressure  peak  and  a  stable  solution  at  a  =  12°  in 
contrast  to  the  steady  case  result.  For  both  angles  of  attack,  negative  CfS  are 
predicted  less  than  0.3%  chord  forward  of  transition  onset,  indicating  a  laminar 
flow  separation. 

Computed  and  measured  pressure  distributions  in  Figure  7.21  are  for 
a  =  10°,  during  the  airfoil  upstroke  and  downstroke.  Computed  skin-friction 
distributions  are  included  in  the  plot.  It  is  apparent  from  the  pressure  distribution 
differences  at  10°  angle  of  attack  during  airfoil  downward  motion  that  the 
computations  did  not  predict  the  shedding  of  the  dynamic  stall  vorticity  that  was 
described  earlier.  The  transition  model  was  developed  for  steady  flows  void  of 
the  changing  pressure  gradients  that  are  part  of  the  oscillating  airfoil  flow  at  high 
angles  of  attack.  The  unsteady  airfoil  motion  poses  many  problems  for  the 
numerical  modeling  as  shown  by  the  abrupt  changes  in  the  skin-friction 
coefficient  in  the  separation  bubble  region.  This  aspect  of  the  problem  requires 
much  further  study  and  therefore  is  omitted  from  further  discussion  in  this  work. 
Skin-friction  coefficients  become  negative  at  x/c  =  0.0135  for  10°  angle  of  attack 
(up)  and  at  x/c  =  0.0116  for  10°  angle  of  attack  (down),  indicating  a  laminar 
separation  bubble. 
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Figure  7.20.  Computed  and  PDI  Derived  Pressure  Distributions  Near  the  Airfoil 
Leading  Edge.  M00=  0.3,  a  =10°-  2°sincot=  8  and  12°,  k=  0.1,  Re=  540,000, 
(x/c)tr=  0.025  and  0.0 12,  Gy    =200. 
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Figure  7.21.  Computed  and  PDI  Derived  Pressure  Distributions  Near  the  Airfoil 
Leading  Edge.    M^  =  0.3,  a  =  10°  -  2°sincot  =  10°(up  /  down),  k  =  0.1,  Re  =  540,000, 
(x  /  c)tr  =  0.025  and  0.020,  GTtr  =  200. 
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3.         Transition  Onset,  Flow  Separation  and  Reattachment 

Figure  6.27  showed  the  specified  transition  onset  locations  for  both  the 
steady  and  unsteady  N-S  analysis,  while  holding  the   Gy     value   at   200. 

Transition  onset  clearly  moved  forward  as  the  angle  of  attack  increases  in  both 
cases.  Unsteady  transition  onset  locations  were  specified  aft  of  steady  locations 
on  the  airfoil  upstroke  and  forward  of  the  same  during  the  airfoil  downstroke. 
Only  a  small  "window"  of  solution  convergence  was  possible  at  a  =  12°, 
suggesting  a  flow  on  the  verge  of  stall  and  presenting  some  of  the  difficulties  in 
numerically  modeling  the  flow  physics. 

Flow  separation  and  reattachment  locations  along  the  airfoil  surface  as 
predicted  by  the  N-S  computations  are  given  in  Figure  7.22.  As  expected,  flow 
separation  is  predicted  to  move  forward  as  the  angle  of  attack  increases.  Using 
negative  skin-friction  coefficients  to  indicate  flow  separation  for  the  low- 
frequency  oscillating  airfoil  flow  as  in  the  steady  flow  analysis,  flow- separation 
location  differences  are  shown  to  be  relatively  minor.  Prediction  of  similar  flow- 
separation  locations  is  consistent  with  the  low  frequency  airfoil  oscillation.  Flow 
reattachment  moves  downstream  as  the  angle  of  attack  increases  and  upstream  as 
the  angle  of  attack  decreases.  Also  note  that,  during  the  upstroke  at  angles  of 
attack  exceeding  10°,  reattachment  becomes  difficult  to  define  because  of 
multiple  separation  bubbles  that  are  predicted  by  the  numerical  solution.  Steady 
angle  of  attack  computations  predict  a  separation  bubble  length  along  the  airfoil 
surface  from  2.7  to  5.9%  chord  at  a  =  8°.  At  a  =  12°,  the  separation  bubble 
length  is  predicted  to  increase  in  size,  extending  from  1  to  7%  chord  downstream. 
The  unsteady  airfoil  flow  separation  bubbles  were  similar  in  size  to  slightly  greater 
than  the  steady  results  for  angles  of  attack  less  than  10°.  Thereafter,  multiple 
areas  of  reverse  flow  are  incorrectly  predicted  by  the  N-S  calculations. 

A  summary  of  specified  transition  onset  locations  and  those  predicted  by 
Michel's  criterion  for  both  the  steady  flow  and  the  low-amplitude  oscillations  are 
given  in  Table  7.1.  Likewise,  the  transition  length  values  predicted  by  Cebeci's 
correlation  formula  are  given  in  Table  7.1.  The  N-S  solutions  for  transition  onset 
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Figure  7.22.  Flow  Separation  and  Reattachment  Estimates  for  Steady  and 
Oscillatory  Computations.   M^  =  0.3,  a  =  10°  -  2°  sin  cot,  k  =  0  and  0.1. 

that  most  closely  matched  the  measurements  were  1%  to  1.5%  chord  downstream 
of  those  predicted  by  Michel's  criterion.  This  distance  downstream  along  the 
airfoil  surface  is  appreciable  since  the  predictions  for  transition  onset  were 
generally  very  close  to  the  airfoil  leading  edge  from  1%  to  2%  chord.  Cebeci's 
correlation  formula  for  predicting  transition  length  values  did  not  vary  widely 
throughout  the  angle  of  attack  range.  The  specified  transition  length  value  of 
200  was  approximately  0.50  of  the  values  predicted  by  the  formulation  that  gave 
secondary  separation  bubbles,  unlike  the  measurements. 


161 


Table  7. 1  N-S  Code  Transition  Onset  and  Transition  Length  Criteria. 


Flow  Over  Steady  Airfoil 

Oscillatory  Airfoil  Flow 

Angle  of 

Attack,  deg 

(Oscill.  - 

up/down) 

(x/c)tr, 
(Michel) 

GY„' 

(Cebeci) 

(x/c)tr 
Specified, 
GY  =200 

la 

(x/c)tr, 
(Michel) 

Or.- 

(Cebeci) 

(x/c)tr 

Specified, 

Gv  =200 
it 

6 

8  (bottom; 

9  (up) 

10  (up) 

11  (up) 

12  (top) 
11  (down) 
10  (down) 

9  (down) 
8  (bottom) 

0.025 
0.023 
0.017 
0.0155 
0.0135 
0.009 

425 
417 
400 
390 
381 
360 

0.040 
0.035 
0.025 
0.0225 
0.0155 
0.0117 

0.015 

0.014 

0.013 

0.012 

0.002 

0.013 

0.014 

0.0145 

0.015 

410 
410 
400 
380 
370 
380 
390 
400 
410 

0.025 
0.025 
0.025 
0.0183 
0.0117 
0.0145 
0.020 
0.0245 
0.025 

4.         Velocity  Profile  Comparisons 

Computed  and  measured  U-component  velocities  are  compared  for  angles 
of  attack  of  8°,  9°,  10°,  and  11.2°  during  the  airfoil  upstroke.  As  in  the  steady 
angle  of  attack  presentations,  velocity  profile  comparisons  begin  at  x/c  =  -0.167, 
which  is  0.50  inch  forward  of  the  airfoil  leading  edge,  and  then  proceed  over  the 
leading  edge  of  the  airfoil  from  0.00  <  x/c  <  0.167  in  normalized  chord  intervals  of 
0.033.  Both  measured  and  computed  U  velocities  are  again  normalized  by  the 
outer  velocity  (U0)  at  y/c  =  0.0167.  An  interpolation  scheme  was  used  to  predict 

computed  velocities  at  the  nodes  of  the  measurement  grid.  Again,  as  in  the  steady 
flow  results,  measured  and  computed  velocities  forward  of  the  airfoil  at  x/c  = 
-0.167  are  in  full  agreement. 
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Comparisons  of  velocity  profiles  for  a  =  8°  (bottom  of  the  oscillatory 
cycle)  are  given  in  Figure  7.23.  Solid  lines  in  the  figure  represent  the  computed 
solution  and  the  measurements  are  depicted  by  symbols.  Normalized  velocities 
range  from  0.0  to  1.44  with  the  maximum  computed  velocity  at  x/c  =  0.033  being 
slightly  greater  than  the  steady  maximum  velocity  at  the  same  location  given  in 
Figure  7.13.  Maximum  measured  velocities  in  the  flow  field  (1.1  U0)  are  at  both 

the  airfoil  leading  edge  and  at  3.3%  chord  downstream.  Generally,  the  data 
agreement  is  within  5%  throughout  the  flow  field.  Exceptions  are  the  measured 
velocities  at  the  lowest  heights  and  3.3  and  6.7%  chord  downstream  that  are 
38%  lower  than  the  predictions  and  up  to  40%  lower  (0.82  U0  vice  1.25  U0)  than 

the  steady  velocities  at  the  same  locations  in  Figure  7.13.  The  computations 
predict  a  slowed  flow  region  only  very  close  to  the  airfoil  surface.  As  previously 
discussed,  the  reduced  velocity  region  points  toward  the  separation  bubble 
region  for  the  measurements,  where  the  surface  of  the  airfoil  is  from  0.058  <  y/c  < 
0.064.  Thus,  an  estimated  separation  bubble  height  for  the  measurements  is  from 
2  to  2.5%  chord  above  the  airfoil  surface. 

Figure  7.24  compares  velocity  profiles  at  a  =  8.8°  (up).  Normalized 
velocities  range  from  0.0  to  1.46  with  the  maximum  velocity  predicted  at  x/c  = 
0.033  and  y/c=  0.063,  which  is  only  0.3%  chord  above  the  airfoil  surface.  The 
maximum  measured  velocity  in  the  flow  field  (1.096Uo)  is  above  the  airfoil 

leading  edge  at  y/c  =  0.083  or  approximately  4.5%  chord  above  the  airfoil 
surface.  Generally,  data  agreement  is  within  5%  throughout  the  flow  field,  except 
in  the  region  where  the  reduced  measured  velocities  point  toward  the  separation 
bubble.  Measured  velocities  for  the  forced  oscillation  at  the  lowest  heights  and 
3.3  and  6.7%  chord  downstream  are  as  much  as  45%  less  than  the  predictions. 
Compared  with  the  previous  results  at  the  lower  angle  of  attack,  measured 
velocities  are  lower  by  5  to  10%  in  the  region  from  3.3  to  10%  chord  downstream 
at  the  lower  heights.  Because  of  the  increased  angle  of  attack  and  the  larger 
region  of  reduced  velocities,  the  estimated  separation  bubble  height  for  the 
measurements  at  a  =  8.8°  (up)  remains  from  2  to  2.5%  chord  above  the  airfoil 
surface. 
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Figure  7.23.  Computed  and  Measured  Velocity  Profiles  Near  the  Airfoil  Leading 
Edge.  M^O.3,  a  =  10°  -2°sincot  =  8°,  k  =  0.1,  Re  =  540,000. 
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Figure  7.24.   Computed  and  Measured  Velocity  Profiles  Near  the  Airfoil  Leading 
Edge.    MTO=0.3,  a  =  10°  -2°sincot  =  8.8°(up),  k  =  0.1,   Re  =  540,000. 
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Comparisons  for  a  =  10°  (up)  are  shown  in  Figure  7.25.  Maximum 
predicted  velocities  have  increased  to  1.52  U0,  which  is  slightly  less  than  the 
steady  results  and  occurs  at  3.3%  chord  downstream  and  0.37%  chord  above  the 
airfoil  surface.  Measured  velocities  are  reduced  to  0.67  and  0.69  U0  at  the  lowest 
heights  and  3.3  and  6.7%  chord  locations  downstream.  Increased  differences 
between  the  data  points  are  seen  in  the  region  of  the  reduced  velocities  as 
recorded  by  the  measurements  that  extend  to  greater  heights  in  the  flow  field. 
Appreciable  reductions  in  the  measured  velocities  are  noted  at  y/c  =  0.100  or 
approximately  3%  chord  above  the  airfoil  surface. 
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Figure  7.25.  Computed  and  Measured  Velocity  Profiles  Near  the  Airfoil  Leading 
Edge,   MM  =  0.3,  a  =  10°  -  2°  sin  cot  =  10  degrees  (up),  k  =  0.1,  Re  =  540,000. 


For  a  =  11.2°  (up)  in  Figure  7.26,  the  data  comparisons  show  increased 
differences  with  even  slower  velocities  measured  by  the  LDV  system,  extending 
to  greater  heights  in  the  flow  field.  Predicted  velocities  over  the  airfoil  increased 
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appreciably  to  1.61  U0,  which  occurred  at  3.3%  chord  downstream  and  0.35% 
chord  above  the  airfoil  surface.  Measured  velocities  are  reduced  to  0.64  U0  at  the 
lowest  heights  and  3.3  and  6.7%  chord  locations  downstream.  Appreciable 
reductions  in  the  measured  velocities  are  noted  at  y/c  =  0.117,  or  approximately 
4%  chord  above  the  airfoil  surface. 


0.181 
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12         3         4         5         6         7 
Stations,  U/Uo  (values  0.0-1.61) 


Figure  7.26.  Computed  and  Measured  Velocity  Profiles  Near  the  Airfoil  Leading 
Edge.  Moo=0.3,  oc  =  10°  -2°sincot  =  11.2° (up),  k  =  0.1,  Re  =  540,000. 


In  summary,  this  study  is  among  the  first  to  critically  evaluate  the  effect  of 
transition  onset  locations  and  transition  length  on  computed  results  by  making 
direct  comparisons  with  experimental  measurements  under  identical  conditions. 
The  results  of  the  work  were  first  reported  in  Reference  69.  PDI  derived  and 
computed  pressure  distributions  are  in  excellent  agreement  for  the  steady  flow 
analysis,  but  for  the  unsteady  comparisons  the  predicted  separation  bubble  starts 
forward  of  and  is  much  smaller  in  height  than  the  measured  bubble. 
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VHI.   CONCLUSION  AND  RECOMMENDATIONS 

The  formation,  development  and  effect  of  leading  edge  separation  bubbles 
on  the  stalling  characteristics  of  a  thin  airfoil  was  the  focus  of  study  in  this  thesis. 
The  LDV  velocities  and  PDI  density  contours  with  resultant  pressure  distributions 
were  used  to  evaluate  computed  results  from  boundary-layer/panel  and  N-S  (thin 
layer,  compressible,  Reynolds-averaged)  codes.  The  flow  considered  for  these 
comparisons  was  (1)  predominated  by  leading-edge  separation  bubbles  at  the 
higher  angles  of  incidence,  and  (2)  compressible  in  a  transitional  Reynolds  number 
range. 

Both  the  PDI  and  LDV  measurements  showed  that  the  laminar  separation 
bubble  alters  the  flow  characteristics  considerably.  The  experiments  indicated  that 
the  separation  bubble  formed  at  7  to  8  °  angle  of  attack  in  both  the  steady  and 
oscillating  airfoil  tests.  The  PDI-derived  pressure  plateau,  corresponding  to  the 
surface  extent  of  the  bubble,  ranged  from  2  to  6%  chord,  decreased  slightly  in  size, 
and  moved  slightly  forward  on  the  airfoil  surface  with  increasing  angles  of 
incidence  as  a  result  of  the  increased  adverse  pressure  gradients.  LDV  velocities 
were  as  high  as  1.68LL  at  12°  angle  of  attack  in  a  region  from  5  to  9%  chord, 
enveloping  the  separation  bubble. 

Light  dynamic  stall  was  observed  briefly  on  the  downstroke  (below  the 
static  stall  angle  of  attack)  for  the  low-amplitude  oscillating  airfoil  flow  in  both  the 
PDI  and  LDV  tests.  This  surprising  and  new  result  was  explained  by  vorticity 
balance  arguments.  The  vorticity  was  partially  shed  on  the  airfoil  downstroke  as 
the  flow  adapted  to  the  rapidly  decreasing  angle  of  attack.  Both  steady  and 
dynamic  stall  was  initiated  from  the  bursting  of  the  leading-edge  bubble  as 
opposed  to  trailing-edge  flow  reversal,  pointing  to  the  need  for  inclusion  of  flow 
transition  modeling  in  the  computations.  For  the  oscillating  airfoil,  the  LDV  and 
PDI  measurements  were  among  the  first  attempted  and  successfully  made  for  the 
flow  conditions  specified. 

The  highly  separated  and  complex  flow  physics  associated  with  the  high- 
amplitude,  deep  dynamic  stall  made  it  intractable  for  detailed  analysis  by 
computational  methods.   However,  the  viscous  effects  were  much  smaller  for  the 
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low-amplitude  tests,  providing  an  opportunity  to  apply  different  computational 
analysis  methods.  Thus,  the  computational  studies  were  limited  to  steady  and  low- 
amplitude,  light  dynamic  stall  flow  cases. 

The  computations  were  the  first  to  incorporate  the  flow  physics  of  leading 
edge  separation  bubbles  and  to  make  comparisons  with  measured  results  in  flows 
with  substantial  pressure  gradients  near  the  stalling  angle  of  attack.  Both  the  direct 
boundary  layer  and  thin-layer  N-S  codes  accounted  for  boundary  layer  flow 
transition  to  analyze  the  development  of  separation  bubbles  for  identical 
experimental  conditions.  The  Chen-Thyson  transition  model  was  used  to  describe 
the  extent  of  transition  flow  region  and  Michel's  criterion  was  used  to  make  initial 
predictions  of  transition  onset  locations.  Pressure  distributions,  skin-friction 
coefficients,  and  velocity  profiles  over  the  upper  leading  edge  of  the  airfoil  were 
shown  to  change  substantially  with  small  changes  in  (1)  the  transition  onset 
location,  and  (2)  the  flow  transition  length  over  the  airfoil  surface. 

Computations  that  incorporate  the  flow  transition  process  were  found  to  be 
critical  in  predicting  separation  bubbles  that  were  similar  to  the  measured  ones. 
For  both  the  steady  and  low-amplitude  oscillating  airfoil  (8°  <a<12°),  the  N-S 
code  predicted  a  negative  skin-friction  region  that  agreed  well  with  the  beginning 
of  and  extent  of  the  pressure  plateau  derived  from  the  PDI  data.  Overall,  the 
agreement  between  the  computations  and  the  measurements  was  good  with  slight 
differences  observed  in  the  separation  bubble  sizes  during  the  airfoil  upstroke. 
However,  the  N-S  analysis  did  not  predict  the  vortex  shedding  process  that  was 
evident  in  the  measurements  during  the  airfoil  downstroke. 

Transition  onset  predictions  moved  forward  with  increasing  angles  of 
incidence,  and  the  range  of  specified  transition  onset  locations  (0.040  >  x/c  > 
0.012)  was  downstream  of  Michel's  criterion  (0.025  >  x/c  >  0.002).  This  was 
expected  because  Michel's  criterion  was  developed  for  higher  Reynolds  number 
flows.  The  range  of  flow  transition  length  values  and  transition  onset  locations  are 
summarized  in  Figure  8.1  for  NACA  0012  airfoil  flows  in  the  transitional 
Reynolds  number  range  (104<Rec<106).  The  use  of  Equation  4.22  and  Michel's 
criterion  gave  separation  bubbles  that  were  much  smaller  than  the  measured  ones 
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(dashed  curve).    The  characteristics  of  measured  separation  bubbles  were  most 
closely  matched  by  the  predictions  when  transition  onset  locations  were  specified 
downstream  of  Michel's  criterion  and  a  transition  length  value  (Gv    =200)  less 

I  ir 

than  those  predicted  by  Equation  4.22  were  specified  (solid  line). 
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Figure  8. 1  Transition  Onset  and  Transition  Length  Values  for  NACA  0012  Airfoil 
Flow,  Mo,  =0.3,  0°<a<12°,  104<Re<106. 

For  future  work,  more  definitive  experiments  (such  as  skin-friction 
measurements)  are  recommended  to  determine  when  and  where  transition  from 
laminar  to  turbulent  flow  takes  place,  in  either  the  attached  or  separated  boundary 
layer  for  flows  with  airfoil-type  pressure  gradients.  Efforts  are  needed  to 
understand  more  fully  the  vorticity  shedding  and  flow  reattachment  processes, 
which  require  advanced  turbulence  modeling.  Transition  length  models  should 
account  for  free-stream  turbulence  levels,  flow  curvature  effects,  and  rapidly 
changing  and  steep  pressure  gradients.  For  higher  speed  flows,  an  understanding 
of  shock  effects  on  the  behavior  of  separation  bubbles  is  required.  Systematic 
experimental  and  analytical  studies  of  the  unsteady  flow  effects  on  laminar 
separation  and  transition  should  complement  each  other  as  in  this  work. 
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This  research  is  the  first  work  to  critically  evaluate  the  effect  of  transition 
onset  locations  and  transition  length  on  computed  results  by  making  direct 
comparisons  with  measurements  for  identical  conditions.  For  the  oscillating 
airfoil  comparisons,  the  predicted  separation  bubble  started  forward  of  and  was 
smaller  in  height  than  the  measured  bubble  during  the  airfoil  upstroke.  The 
computational  schemes  need  more  work  to  correctly  predict  the  vortex  shedding 
and  flow  reattachment  processes  that  were  evident  during  the  downstroke  of  the 
oscillating  airfoil.  The  PDI-derived  and  computed  pressure  distributions  as  well  as 
the  LDV  and  computed  velocities  were  in  excellent  agreement  for  the  steady  angle 
of  attack  flow  analysis. 
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